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Wu  Wenquau  and  Liu  Cule 
(Academia  Sinica) 


ABSTRACT 


Based  on  the  basic  three-dimensional  flow  thermo  aerodynamic 
equations  of  turbine  machinery  proposed  by  Wu  Chung-Kua  employing 
non-orthogonal  curvilinear  coordinates  and  their  corresponding 
non-orthogonal  velocity  components,  the  main  function  to  be  solved 
i.e.,  the  second  order  quasillnear  partial  differentiation  equation 
of  the  stream  function,  was  obtained  in  this  paper  by  Introducing 
a  stream  function.  With  the  exception  of  density  related  terms, 
the  various  order  of  derivatives  of  the  stream  function  were 
placed  on  the  left  side  of  the  equation.  By  doing  so,  the 
convergence  rate  was  accelerated.  After  discretization  of  the 
differential  equation  by  using  a  center  nine-point  difference 
format,  the  linear  algebraic  equations  obtained  were  solved  directly 
by  the  decomposition  of  the  matrices  [L]  and  [u].  The  convergence 
rate  of  this  method  is  comparatively  faster  than  other  methods. 

The  coefficient  matrix  is  a  diagonal  band  shaped  scattered  matrix. 

We  adopted  the  two  measures  that:  1)  non-zero  elements  were  numbered 
according  to  the  diagonal  line,  and  2)  fictitious  points  were 
added,  to  greatly  reduce  the  internal  storage  of  the  computer. 

When  calculating  the  density  from  the  stream  function,  an  internally 
stored  density  function  table  interpolation  method  was  adopted. 

The  selection  of  relaxation  factor  was  briefly  discussed.  This 
program  was  used  in  the  calculation  of  the  blade  patterns  of 

f '  - - 

This  paper  was  presented  in  the  1976  National  Turbine  Machinery 

Aerodynamic  Conference. 

*•  Numbers  in  the  margin  indicate  pagination  of  foreign  text. 
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some  compressors  and  turbines.  In  comparison  of  the  experimental 
results  and  the  theoretical  analytical  solutions,  they  were  In 
mutual  agreement. 


TABLE  OF  SYMBOLS 


e^^  base  vector  of  a  coordinate  in  the  coordinate  system 
e^  reciprocal  vector  of  e^ 

a^^j  forward  transformation  tensor  element  of  the  x^  coordinate  system 


a 


the  Jacobrau  determinant 


composed  of  a 


ij 


I  the  enthalpy  of  stagnant  gas  rotor 


x^  an  arbitrary  non-orthogonal  curvilinear  coordinate  system 

w^  the  reverse  transformation  component  of  the  relative  velocity 
vector  of  the  gas 


w^  the  forward  transformation  component  of  the  relative  velocity 
vector  of  the  gas 

the  physical  component  on  the  x^  coordinate  of  the  relative 
velocity  vector  of  the  gas 

u^  unit  vector  of  the  x^  coordinate  system 

0^  the  angle  between  the  x^  coordinate  line  and  the  revolving 
parent  line 

2 

02  the  angle  between  the  x  coordinate  line  and  the  revolving  line 
Z  the  density  function 


a  function  of  the  partial  derivatives  of  the  stream  function 


V  the  relative  motion  stream  function 
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Superscript : 
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o  stagnant  parameter 
'  relative  stagnant  parameter 
*  dimensionless  parameter 

Subscript : 

i  inlet  station  parameter 

J  calculation  station  along  the  direction 

2 

k  calculation  station  along  the  x  direction 

INTRODUCTION 

The  analysis  and  calculation  of  the  cascade  deflected  flow, 
i.e.,  the  computation  of  the  S^^  relative  stream  surface,  is  an 
important  aspect  in  the  three-dimensional  flow  computation  of 
turbine  machinery.  It  not  only  is  a  direct  solution  to  the  blade 
pattern  problem,  but  also  is  an  indispensible  part  in  the  alternating 
iteration  of  the  S^  and  S2  types  of  relative  stream  surface  in 
a  three-dimensional  flow  computation.  By  doing  so,  not  only 
the  flow  conditions  upstream  and  downstream  from  the  blade  are 
satisfied,  but  also  the  flow  conditions  of  the  gas  flow  in  the 
blade  channel  are  ensured. 

In  References C 3-51,  the  basic  equations  of  three-dimensional 
flow  in  turbine  machinery  were  presented  by  employing  arbitrary 
non-orthogonal  curvilinear  coordinates  and  their  corresponding 
non-orthogonal  velocity  components.  Its  characteristics  are  that 
it  can  adapt  to  any  complex  shape  and  type  of  turbine  machinery, 
and  it  is  capable  of  rigorously  satisfying  the  tangential  condition 
between  the  flow  and  the  boundary  along  the  boundary  of  the  object. 
Based  on  this  set  of  basic  equations,  this  paper  presented  the  actual 
computation  method.  Because  of  the  non-orthogonal  curvilinear 
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coordinates,  it  also  Involved  the  computation  of  the  geometric 
characteristics  of  the  coordinates.  These  computation  methods 
belonged  to  the  original  type. 

We  compiled  the  computer  program  on  the  6912  machine  at 
Beijing  University.  Its  internal  storage  is  64K.  The  allowable 
computational  nodal  points  reached  900  points  (approximately  15 
points  between  cascades),  which  is  sufficient  for  the  usual 
engineering  requirements.  If  the  capacity  of  the  computer 
is  less,  it  is  possible  to  reduce  it  further.  Finally,  some 
typical  computational  results  are  given.  A  comparison  of  the  ex¬ 
perimental  results  and  the  theoretical  analytical  solutions  was 
carried  out;  they  were  in  mutual  agreement. 


I.  The  Basic  Set  of  Equations 

Assuming  the  motion  relative  the  blade  is  stable  and  adiabatic, 

n  21 

then  the  basic  equations  are  *  : 

Continuity  equation 

V  .  (^,W)  =  0  ( 1 ) 

Equation  of  motion 

W  X  (V  X  V)  —  TVs)  (2 ) 

First  Law  of  Thermodynamics  — -  energy  equation 

dl/dt  =  W  •  V/  0  (3) 

Where 


/  =  A  +  iW|V2  +  1U1V2  ('<) 


Second  law  of  thermodynamics 


Tdsldt^fs  (5) 


Ideal  gas  state  equation 


p  =  pRT  (6) 


In  the  derivative  of  the  original  equations,  the  equation  of 
motion  took  the  vlscolsty  of  the  gas  Into  account  and  it  contained 
a  viscosity  term.  Due  to  the  viscosity  effect  and  the  shock  wave 
loss,  an  enthalpy  increase  was  produced.  Therefore,  the  equation 
set  is  consistent.  However,  in  the  actual  computation,  the 
viscosity  term  was  neglected  in  the  equation  of  motion.  It  is 
due  not  only  to  the  complexity  of  its  computation,  but  also  to 
the  fact  that  it  is  a  small  term  comparatively  speaking.  Thus, 
the  approximate  relations  in  equations  (2)  and  (3)  were  obtained. 
The  gradient  of  enthalpy  in  the  equation  was  calculated  through 
the  selection  of  empirical  data  such  as  the  efficiency  or  pressume 
loss  coefficient,  and  the  loss  during  passing  through  shock  waves, 
etc.  The  above  equations  were  expressed  in  an  invariant  form  with 
respect  to  the  coordinates.  They  are  applicable  to  all  coordinate 
systems.  Furthermore,  they  are  simple  and  easy  to  understand. 


II.  The  Coordinate  System  and  the  Velocity  Component  Equation 

By  choosing  an  arbitrary  non-orthogonal  curvilinear  coordinate 
12  9 

system  (x  ,  x  ,  x"^)  and  substituting  the  correlations  of  the  degree 
of  revolution  into  the  basic  equations  in  the  previous  section, 
it  is  possible  to  obtain  the  expression  for  each  component. 
Specially  with  regard  to  the  direct  problem  of  the  S^  relative 
stream  surface  calculation,  because  the  shape  of  the  stream  surface 
is  knr  rn,  then  It  is  possible  to  choose  this  known  stream  surface 
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as  a  coordinate  plane  x  -constant.  Furthermore,  the  x  coordinate 

was  made  perpendicular  to  this  plane  as  shown  in  Figure  1.  The 

1  2 

coordinate  lines  (x  ,  x  )  are  shown  in  Figure  2.  Because  the  iso- 
x”^  planes  are  the  stream  surfaces,  the  velocity  vector  must  be 
tangent  to  this  stream  surface.  This 


(7) 


For  the  usual  engineering  requirement ,  it  is  possible  to 
assume  that  this  stream  surface  is  an  arbitrary  revolving  surface. 
Prom  equations  (2)  and  (7),  we  can  obtain  the  component  along  the 
^2  direction  to  be 


—  W 


dx^J 


4“  2(0 


dx^ 


(8) 


In  order  to  more  precisely  show  the  physical  significance  and 
to  facilitate  the  consistent  computation,  the  and  in  equation 
(8)  were  further  expressed  by  the  physical  component  W^; 


tVi 


(9) 


The  distance  between  two  neighboring  iso-x  planes  (l.e.  , 
two  neighboring  stream  surfaces)  is  t.  In  addition,  let  dx  =1 
then 


-v/ajj  = 


(10) 


Due  to  the  x^  coordinate  line  being  chosen  to  be  perpendicular 
to  this  stream  surface,  then 


““  Ob  Of]  =•  0 


(11) 


Prom  the  definitions 


Oil 


dx' 


dx‘ 


9(>'‘P)  dCrip) 
dx'  dx‘ 


V  •hi'/  O]]  cosfl,, 


(12) 
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a  =  s'ft'  ®'2  sill’  ^12 


(13) 


From  equation  (9) 


\2^ 

This  Is  the  form  of  equation  which  is  familiar  to  us.  Equation 

1  2 

(16)  contained  two  relative  velocity  physical  components  W  and  W 
which  were  to  be  determined.  These  could  be  solved  by  using  the 
so-called  streamline  extension  method.  However,  the  more  suitable 
method  was  to  Introduce  a  stream  function  ¥  to  further  transform 
the  equation. 


Figure  1.  The  arbitrary  revolving  Figure  2.  The  arbitrary  non¬ 
plane  of  the  relative  stream  orthogonal  curvilinear  coor- 
surface.  dinate  system  (x"^,  x*^). 
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III.  stream  Function  and  Stream  Function  Equation 


O 

From  the  continuity  equation  (1)  and  Vf  =0 ,  we  get 


Or  Ux 


Define  the  stream  function  T  which  satisfied  the  following 
relation 

sill  =  -- — ;  a  22  sin  0,2  MQ^ 

Or'  Or'  \  -i-y  j 

By  substituting  equation  (19)  Into  (l6),  we  can  obtain  the 
stream  function  equation.  In  order  to  calculate  conveniently, 
the  equation  was  further  rendered  dimensionless .  The  dimensionless 
reference  quantities  are:  'F*  = 'K/g  •  loo;  p*  =  p/p'.-,  =  p*  ^  p/p'-, 

T*  =  r/'r/,  r*  =  1-/V,;  a>*  =  mrJU,  =  1;  =  V  a,,/'',; 

v'iiTi*  =  (»")•  = (»")*  =  II  vt/.;  (20) 

l* I  •  g.Ji'V]-,  H*=-H-sJ/VV, 

where  G  Is  the  rate  of  flow  through  a  blade  duct,  and  U,  Is  the 

r  61  ^ 

revolving  speed  at  the  Inlet.  We  get*-  - 

1  o'y*  2cosfl„  a"F*  ,1  d‘>P*  ^  VI 

a,*.  Or')'  Vafi^/aS  aS  ( Or')'  y/ 7 a*  ^r'  (21) 

where 


^/  flu  ‘ 

\/a5ar' 

sin0i2V^ 

s/ sin  0|j) 

^  COS0I2 d  In  T* 

,  1  Q&n 

V'^Or' 

\J  a,*  Or' 

sin0,2\/ 

/  1  0ln/>*  _ 

COS  0,2  ^lnjo'*‘\ 

ov* 

V^i 

^^idr' 

,  /  1  Olnp* 

cos0,2  dlnp* 

\  a<F* 

a*' 

V^9r' 

-  2r V'r  sin  (T  sin'  6,2  +  (r%*  sin  0,,)'  ^  ^ 
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where 


r*=  Jll  ii  Aip.  cos,?,  (  1  +  - - 


is  a  dimensionless  constant  which  represented  the  various  incoming 
gas  states,  i.e.,  corresponding  to  various  cascade  density  A41  , 
gas  entering  angle  6^,  the  relative  M  number  of  the  entering  gas 
Mwi ,  and  the  revolving  speed  U^,  etc.  The  coefficients  in  equation 
(21)  were  the  only  gecanetrlcal  parameters  of  the  coordinate  system. 
After  the  mesh  was  determined,  these  coefficients  only  needed  to 
be  calculated  once.  The  major  terms  are  already  included  on  the 
left  side  of  the  equation.  On  the  right  hand  side,  there  are  non¬ 
linear  terms  due  to  variation  of  density.  These  require  to  be 
calculated  through  "iteration".  Usually,  the  value  obtained  from 
the  previous  iteration  process  is  chosen.  In  reality,  the  original 
second  order  quasllinear  partial  differential  equation  is  trans¬ 
formed  into  a  second  order  linear  partial  differential  equation  and 
the  solution  is  sought  through  Iteration. 


The  qualitative  determining  formula  of  equation  (21)  is 


A  =“  1/ cos^  ^a/‘>n<tn  ^0  (  22  ) 

We  can  see  that  it  is  an  elliptical  equation  which  is  suitable 
for  subsonic  flow.  The  original  basic  equation  was  not  limited 
by  this  condition.  The  problem  is  located  in  the  treatment  of 
the  density  terms.  As  a  matter  of  fact,  density  is  related  to 
velocity  which  Includes  the  first  order  derivative  of  4'. 

Therefore,  the  right  hand  tenns  actually  contain  a  second  order 
partial  derivative  term  of  'i', which  affects  the  equation  qualitat¬ 
ively.  If  the  right  hand  terms  of  equation  (21)  are  all  moved  to 
the  left  hand  side,  then  the  equation  changes  from  an  elliptical 
type  to  a  hyperbolic  type  as  the  flow  velocity  changes  from 
subsonic  to  supersonic.  However,  our  present  treatment  is  only 
suitable  for  subsonic  flows.  If  we  entered  the  supersonic  region 


( 
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in  the  computation,  the  program  has  an  automatic  alarming 
measure  arranged  in  it. 

In  order  to  solve  equation  (21),  it  is  also  possible  to 

2 

keep  the  partial  derivatives  with  respect  to  x  on  the  left  hand 

side  and  the  remaining  terms  are  moved  to  the  right  hand  side. 

Thus,  the  equation  is  transformed  into  a  normal  differential 
2 

equation  of  x  in  looking  for  the  solution.  This  method  is  the 
so-called  stream  line  extension  method.  The  actual  computation 
showed  that  when  the  computation  mesh  is  denslfled  (mainly 
in  the  streamline  direction)  then  the  convergence  rate  of  the 
flow  field  matrix  method  is  much  faster  than  that  of  the  stream¬ 
line  extension  method.  The  reason  is  that  a  large  number  of 
undetermined  terms  are  moved  to  the  right  hand  side  in  the  streamline 
extension  method  to  be  substituted  by  the  values  obtained  in  the 
previous  iteration  process.  Thus,  this  large  lagging  quantity 
is  affecting  the  convergence  rate.  Therefore,  the  flow  field 
matrix  method  is  used  to  solve  the  equation  in  the  following. 


IV.  The  Computational  Mesh  and  Boundary  Conditions  /22 

For  the  usual  S,  flow  surface  calculation  it  is  possible  to 
1  2  ^ 

choose  the  (x  ,  x  )  coordinate  lines  as  shown  in  Figure  1,  in  which 

2  1 
the  x  line  is  a  circular  arc  and  the  x  line  is  the  line  connecting 

all  the  points  dividing  the  width  of  the  blade  duct  according  to 
the  same  proportion  in  the  blade  duct.  In  the  upstream  and  down¬ 
stream  regions  of  the  blade,  it  can  be  simply  chosen  as  a  straight 
line  on  the  (J)' -  1  plane.  The  angles  of  inclination  with  respect 
to  line  1  are  3^  and  Bg,  respectively.  The  computational  mesh 
formed  this  way  only  will  require  the  input  of  a  few  data;  the 
parent  line  r(z)  of  the  revolving  surface,  the  blade  surface 

coordinates  ({>  (z)  and  ({>  (z)  the  gas  entrance  and  exit  angles 
P  ® 

0-, ,  So.  It  is  very  convenient.  If  the  computational  mesh  is  plotted 
^  ^  \  2 

on  the  X  —  x  plane,  it  is  as  shown  in  Figure  3.  Thus,  it  means 
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that  the  computational  mesh  is  standardized. 


Figure  3.  The  standardized  computational  mesh. 

Key:  1)  inlet  zone;  2)  outlet  zone;  3)  fictitious  point; 

4)  blade. 


In  the  computation,  it  is  also  required  to  obtain  some 
geometric  parameters  such  as  ^22  ’  ®12»  Numerical 

differentiation  was  also  used  here.  To  obtain  ^a^^ ,  there  are 
two  types  of  correlations: 

(23) 


Our  experience  is  that  it  is  more  difficult  based  on  equation 

(24).  Theoretically  speaking,  the  values  of  ya, ,  and  so  on  can 

1 

vary  with  the  arbitrary  selection  of  x  .  However,  when  considering 
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the  error  of  numerical  differentiation,  it  is  better  off  to  have 
the  values  of  v/a^^ ,  \/a^,  etc.  be  close  to  1.  Therefore,  we 
selected  the  line  with  x^=0  on  the  center  line  of  the  blade 
duct.  Furthermore,  we  chose  this  line. 

Because  of  the  use  of  a  curvilinear  coordinate  and  also  because 
of  the  introduction  of  the  stream  function,  the  boundary  con¬ 
ditions  are  very  simple  and  easy  to  understand: 

(1)  along  the  surface  of  the  blade 

The  surface  of  the  blade  is  a  stream  line,  and  T  is  a  constant. 

'^*♦  =  0;  'if*  ^  =  100 

(2)  the  inlet  and  outlet  regions 

\  ^  ^  m  r,  /B]  ^  ^  ^  m 

A  periodical  condition  exists,  i.e. 

For  any  aerodynamic  parameter 

^23 


The  stream  function 


Kk  *  <i>%..k  -  JOO  -  +  100 

(3)  The  Entrance  and  Exit  Gas  Boundaries,  aa’  ,  dd’  . 


Assuming  that  the  entering  and  exiting  gas  boundaries  are 
located  at  an  infinite  distance  away  from  the  blade,  then  the  two 
boundary  flows  are  homogeneous.  Actual  calculations  showed  that 
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it  is  sufficient  If  they  are  In  front  of  and  behind  the  blade  by 
a  cascade  distance.  The  boundary  conditions  here  are: 


(c)=  constant 


9x'  xl,  —  x;  6x‘ 


V.  The  Storage  of  the  Coefficient  Matrix  and  the  Decomposition 
of  [L]  [u] 

The  various  partial  derivatives  in  the  above  equation  (21) 
all  adopted  a  three-point  center  difference  formula  to  form  a  . 
center  nine-point  difference  format  to  discretize  the  equation. 
The  linear  algebraic  equations  can  be  written  as 


=  {p}  (25) 

[M]  is  a  matrix  formed  by  the  coefficients,  {'!'*},  {p}  are  the 
column  vectors.  It  is  possible  to  use  various  methods  to  solve  this 
equation  set  (25)  such  as  point  iteration.  Improved  point  iteration, 
and  linear  relaxation.  This  paper  adopted  the  direct  decomposition 
of  [L]  [u]  method.  The  purpose  is  to  accelerate  the  convergence 
rate.  From  (25) 


=  (p)  (26) 

(27) 

to  solve  and  obtain  <l>*.  The  reality  showed  that  this  method  is 
comparatively  faster. 

The  difficulty  of  this  method  is  that  it  is  necessary  to 
store  the  large  capacity  coefficient  matrix  [M].  Fortunately, 
this  [M]  is  a  diagonally  band  shaped  scattered  matrix.  In  order 
to  reduce  the  internal  storage,  two  measures  were  taken. 
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(1)  Renumbering  the  coefficient  M,  ,  according  to  the  diagonal 

Jj  K 

line. 

If  [M]  is  stored  as  a  two-dimensional  number  group  M[l:mn, 
l:mn],  then  it  is  difficult  to  eliminate  the  large  amount  of  non¬ 
zero  elements. 


i’'  =  i 


Figure  4a.  The  original  matrix.  Figure  4b.  The  matrix  after 

adding  fictitious  points. 

We  renumbered  the  elements  along  the  diagonal  line.  As  shown 
in  Figure  (4a):  J'  represents  the  diagonal  line  number  from  1  to 
4n-l  and  k*  represents  the  element  sequence  along  the  diagonal 
line  direction  from  1  to  mn.  At  this  time,  the  internally  stored 
numbers  are 

N  =  (4n  -  1)  mn 
(2)  -Adding  fictitious  points. 

A  periodical  condition  exists  in  the  inlet  and  outlet  regions 
of  the  blade  which  caused  the  doubling  of  the  non-zero  elements. 

For  this  purpose,  a  row  of  nodal  points  were  added  in  the  inlet  /24 

region  and  outlet  region,  respectively,  as  shown  in  Figure  3. 

These  nodal  points  are  not  the  ones  we  wanted  to  calculate. 

However,  they  are  still  Included  in  the  numbers.  Therefore,  they 
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are  called  fictitious  points.  Thus,  the  band  width  Is  changed 
to  2n+5.  Now,  the  required  Internal  storage  Is 


.V  =  2n‘m  +  (7n  +  5)m 


The  internally  stored  quantity  is  reduced  by 


2n'm  —  (8n  +  5)m 


This  N  increases  with  n  according  to  a  square  relation. 
Thus,  the  internal  storage  is  saved  by  approximately  a  half. 
Along  with  Increasing  n  value,  the  effect  is  especially  obvious. 
If  we. choose  m*50,  n»15,  then  N=28,000.  Similarly,  when 
m,n=10,  then  N=13,750.  Therefore,  this  prograim  can  be  used 
on  a  32K-6ilK  Internal  storage  computer. 


VI.  The  Computation  of  the  Density  p  and  the  Selection  of 
the  Relaxation  Factor 


In  the  density  calculation,  the  variation  of  enthalpy  can 
be  taken  into  account.  However,  in  the  S^^  stream  surface 
calculation,  due  to  lack  of  experimental  data,  usually  it  is 
an  iso-enthalpic  flow.  Therefore,  we  also  treated  it  using 
the  iso-enthalpic  correlation  to  obtain^ 

■  ■  (w  ~  \—!r - )  . 


[l _  .1 _  1 

/*  4-  —  r*' 

f  =  -  "l 

\  h;*  J 

1  2H:*Cr*p*sin0,jr)’ 

1  hC  ^ 

(28) 


//*  +  ^  ^  r,  y  /  d4>^  y 


/24 
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Let 


then 


r  (l  -  <P/I^)fh  (29) 

When  is  solved,  $  can  be  obtained.  Prom  (32)  we  can  find  z, 
and  p*  can  be  further  obtained.  Equation  (29)  is  not  explicit. 

In  order  to  avoid  the  trouble  of  calculation  each  time,  it  is 
possible  to  prefabricate  a  Z— $  table  to  be  stored  in  the  internal 
storage  of  the  computer  in  preparation  for  interpolation  computation. 

'I'  has  a  maximum  0.066  .  This  corresponds  to  M-1  and  the 

dense  flow  maximum  X  o.68.  it  should  also  be  noted  that  there 
are  two  Z  values  (p  values)  corresponding  to  the  same  $  value. 

One  corresponds  to  M<1  and  the  other  corresponds  to  M>1. 

When  solving  equation  (21),  the  right  hand  terras  of  the  equation 

contain  y*  and  p*.  This  requires  us  to  carry  out.  Iterations.  Various 

relaxation  factors  and  must  be  selected,  respectively. 

According  to  our  computational  experience,  the  convergence  of  density 

p  is  better  and  usually  it  is  chosen  that  a  =1.  However,  a,  is 

p  'J' 

chosen  to  be  various  values  under  different  conditions.  When  the 
flow  M  number  is  comparatively  low,  can  be  chosen  to  be  a  larger 
value.  For  example,  in  the  approximate  incompressible  situation, 
then  at|(  *1.  When  the  M  number  increases,  then  aii>  should  be  relatively 
lower,  usually  around  0.5.  When  M  approaches  1,  then  x'P  should  be 
even  lower.  This  is  consistent  with  the  physical  concept,  l.e., 
when  the  flow  is  close  to  the  speed  cf  sound, the  density  flow  value 


/25 
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does  not  converge  easily.  For  the  usual  engineering  requirements 


\  p*io+»  _  <  l()-5 

_  ^•<•'>1  <  10"‘ 

then  10-15  iterations  will  be  sufficient.  The  number  of  iterations 
required  is  even  less  when  the  M  number  is  low. 


VII.  Computational  Examples 

Trial  computations  were  made  for  the  blade  shapes  of  compressors 
and  turbines.  In  order  to  compare  to  the  experimental  results 
or  theoretical  analytical  solutions,  these  are  all  plane  cascades. 

A  number  of  cascade  computations  were  carried  out  and  the  results 
were  all  satisfactory.  The  following  gives  two  typical  results: 

1)  Reference  [7]  provided  a  compressor  cascade  with  an 
incompressible  analytical  solution.  The  relative  thickness  is 
10%  and  the  density  is  relatively  small.  The  comparison  to  the 
calculated  results  is  as  shown  in  Figure  5- 

2)  The  maximum  thickness  of  the  blade  of  the  T^  compressor 
was  shifted  backward  to  S0%  of  the  chord  length  to  thus  reduce  the 
maximum  velocity  peak  of  the  gas  entrance  boundary.  Furthermore, 

a  second  velocity  peak  was  at  6055  of  the  blade  chord.  Thus,  the  allow¬ 
able  inlet  M  number  can  be  improved.  The  comparison  of  the  cal¬ 
culated  results  to  the  experimental  ones  showed  that  they  are 

r  8 1 

consistent  as  shown  in  Figure  S.'- 

Figure  7  shows  the  distributions  of  the  density  flow  on 
the  center  stream  line  and  the  inlet  density  flow  ratio.  It  showed 
that  the  effect  of  blade  thickness  is  not  only  limited  to  the 
blade  duct,  but  also  in  the  blade  inlet  and  outlet  regions.  It  is 
worthwhile  noting  that  these  compressor  blade  results  are  con¬ 
sistent  with  the  turbine  blade  results  obtained  many  years  age. 
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These  must  be  taken  into  account  in  the  design  of  compressors. 
Furthermore,  these  data  are  indispensible  in  the  iterative 
calculation  of  the  and  S2  stream  surfaces. 


Figure  5.  Comparison  between 
the  computed  results  and  the 
theoretical  analytical  solution. 

Key:  1)  numerical  computed 

value;  2)  analytical  solution. 


:.'M  40 


Figure  6.  Comparison  of  computed 
results  to  the  experimental 
results . 

Key:  1)  blade  model;  2)  calculated 
value;  3)  experimental  value. 


U  ‘JO 


:  C 


Figure  7.  The  effect  of  blade  thick¬ 
ness. 

Key:  1)  blade  region;  2)  p  is  cascade 
distance,  t  is  the  blade  thickness. 

(pWT-),.  _  p 

(pWx),  •  p  -  » 
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CONCLUSIONS 


/26 


Summarizing  the  above,  we  believed  that:  1)  the  three-dimen¬ 
sional  flow  basic  equation  set  for  turbine  machinery  which  was 
rigorously  derived  in  Reference  [3]  by  employing  arbitrary 
non-orthogonal  curvilinear  coordinates  and  their  corresponding 
non-orthogonal  velocity  components  can  be  applicable  to  any 
complicated  boundary  conditions.  Actual  computation  showed 
that  it  is  a  feasible  and  effective  method.  The  use  of  this 
coordinate  system  and  its  computational  method  are  original. 
Furthermore,  it  can  be  extended  to  other  engineering  computations. 

2)  By  introducing  the  stream  function,  a  second  order  quasi- 
llnear  partial  differential  equation  of  the  stream  function  'V  was 
established  as  the  major  equation  in  searching  for  the  solution. 
This  not  only  makes  the  boundary  conditions  more  rigorous  and 
precise,  but  also  brings  about  improvements  in  the  numerical 
solution  method. 

3)  In  this  paper,  with  the  exception  of  the  density  terms 
in  the  stream  function  equation  to  be  solved,  the  various  orders 
cf  partial  derivative  terms  of  'V  are  placed  on  the  left  hand  side 
of  the  equation.  This  would  reduce,  to  the  extent  possible, 

the  lagging  quantities  in  the  iteration  computational  process. 

This  is  very  favorable  to  the  acceleration  of  the  convergence 
rate. 


4)  The  direct  decomposition  of  matrices  [L]  [u]  was  used 
to  solve  the  'f  equation.  Due  to  the  fact  that  the  density 
terms  are  shifted  to  the  right  hand  side  of  the  equation,  the 
coefficient  matrix  is  composed  of  geometric  parameters  of  the 
mesh  which  are  Invariant  in  the  Iteration  process.  Therefore, 
the  matrix  decomposition  will  only  have  to  be  carried  out  once. 
Furthermore,  the  density  terms  are  converging  more  easily. 

Actual  experience  showed  that  this  method  converges  comparatively 
faster. 
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5)  In  the  direct  matrix  solution  method,  the  storage  of 
the  coefficient  matrix  is  an  important  problem.  Because  the 
diagonal  line  numbering  and  the  additional  fictitious  points 
methods  were  used,  the  internal  storage  was  greatly  decreased. 

6)  When  finding  density  p  from  the  streaim  function  'i' , 

an  internally  stored  density  function  table  interpolation  com¬ 
putation  was  used.  This  is  favorable  to  accelerating  the  computing 
speed . 

7)  From  the  comparison  of  the  calculated  results  to  the 
analytical  solutions  and  the  experimental  results,  we  can  see 
that  they  are  in  mutual  agreement.  There  is  some  difference 
in  the  gas  entering  boundary.  This  is  mainly  due  to  the 
complicated  variation  of  the  blade  curve  at  the  gas  entering 
boundary,  and  the  simulated  matching  of  the  curve  is  relatively 
more  difficult.  This  will  require  special  treatment  which  will 
be  described  in  another  paper.  In  addition,  for  the  difference 
in  comparison  to  the  experimental  results,  the  local  flow 
separation  is  also  another  reason. 

This  paper  was  carried  out  under  the  direct  guidance  of 
Professor  Wu  Chung-Hua.  We  wish  to  express  our  deepest  gratitude. 
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FLOW  FIELD  LINE  RELAXATION  SOLUTION  FOR  INVERSE  PROBLEiM  OF 
PLOW  ALONG  S2  RELATIVE  STREAM  SURFACE  EMPLOYING 
NON-ORTHOGONAL  CURVILINEAR  COORDINATES  AND  CORRESPONDING 
NON-ORTHOGOMAL  VELOCITY  COMPONENTS* 

Zhurong-guo 
(Academia  Slnica) 


ABSTRACT 


Specifically  with  regard  to  the  calculation  of  the  inverse 
problem  of  the  S2  stream  surface  of  a  turbomachine ,  the  computer 
program  of  the  flow  field  line  relaxation  solution  for  the  inverse 
problem  of  flow  along  the  S2  stream  surface  employing  arbitrary 
non-orthogonal  curvilinear  coordinates  and  non-orthogonal  velo.'ity 
components  was  introduced.  Furthermore,  the  differences 
between  this  program  and  the  presently  commonly  u.-red  vr.l'  -Ity 
extension  method  (streamline  curvature  method,  yrogram  and  the 
direct  matrix  solution  method  program  were  briefly  compared  and 
ccmmented  on.  The  advantages  .  of  the  flow  field  line  relaxation  method 
are  pointed  out.  Especially  when  the  three-dimensional  blade 
design  computation  is  used,  when  it  is  necessary  to  establish 
more  computational  stations  along  the  streamline  direction  of 
the  blade  region,  the  flow  field  line  relaxation  method  is 
worth  recommending. 


I.  Introduction 


Since  Wu  Chung 


^u^l 


established  the  three-dimensional  flow 


basic  equation  of  the  turbomachine  and  the  stream  surface 


theory  for  the  solution  of  three-dimensional  flow  in  the  early 


This  paper  has  been  presented  in  the  Second  National  Engineering 
Thermal  Physics  Conference  in  Haugchow  in  November  1973. 
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fifties,  great  advances  have  been  made  in  the  aerodynamic  numerical 
computations  of  turbomachines.  At  the  present  moment,  it  has 
already  combined  with  experiments  to  become  an  indispensible 
powerful  tool  in  the  design  and  development  of  turbomachines. 

In  the  aerodynamic  design  of  turbomachines,  the  computation 
of  the  S2  stream  surface  occupies  an  extremely  important  position. 

It  determines  the  flow  distribution  along  the  blade  height. 

For  the  calculation  of  the  S2  stream  surface,  much  atten¬ 
tion  was  paid  to  it.  Kany  different  numerical  methods  were  developed. 

In  the  selection  of  the  coordinate  system,  initially  the 
orthogonal  cylindrical  coordinate  system  and  the  corresponding 

r?! 

orthogonal  velocity  components  were  adopted^  .  In  order  to 
adapt  to  the  various  boundary  shapes,  the  three-dimensional 
flow  basic  equation  set^^“^^  of  turbomachlnes  employing  non- 
orthogonal  curvilinear  coordinate  systems  and  non-orthogonal 
velocity  components  was  used.  By  Introducing  a  non-orthogonal 
curvilinear  coordinate  system,  it  can  very  conveniently  transform 
a  physical  surface  with  conplicated  boundary  shape  into  a  rectangular 
mesh  computational  plane.  In  the  meantime,  it  can  naturally 
satisfy  the  physical  condition  that  the  boundary  flow  direction  is 
tangential  to  the  surface. 


In  the  computation  of  the  S2  stream  surface,  from  the 
viewpoint  of  the  equation  form  and  the  computation  method,  there 
are  primarily  two  types  of  methods  widely  used. 


The  first  type  is  the  velocity  extension  method  (streamline 
curvature  method).  It  is  based  on  the  radical  component  equation 
of  the  equation  of  motion  to  write  the  partial  differential  equation 
of  the  velocity  along  the  blade  height.  For  exaunple,  it  is 
written  in  the  following  form^^^: 


dr 


J 


-  y.o  .  ./aw  Os  \  p\.dv. 


(1) 
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The  left  hand  side  of  equation  (1)  is  written  in  the 
velocity  partial  derivatives  form  which  is  characteristic 
of  this  type  of  solution  method.  When  looking  for  the  solution, 
first  a  velocity  value  V  .  is  assumed  at  a  certain  radius 

o  ^  li 

(such  as  the  inner  radius)  and  then  the  velocity  value  along 
the  blade  height  according  to  the  equation.  The  flow  rate  is 
used  to  check  whether  the  selected  initial  velocity  is  proper. 

If  it  is  not  proper,  then  another  Initial  value  is  chosen  until 
the  flow  rate  condition  is  satisfied.  Because  the  right  hand 
of  the  equation  contains  many  unknown  parameters  related  to  the 
streamline  shape,  they  can  only  be  replaced  by  the  approximate  /29 

values.  Thus,  it  requires  repeated  Iteration. 

For  the  quantities  dependent  on  the  streamline  position 
on  the  right  hand  side  of  the  equation,  they  can  be  expressed 
in  various  forms.  Earlier,  such  as  in  Reference  [1],  assuming 
that  the  projection  of  the  streamline  on  the  meridian  plane  is 
a  sinusoidal  curve,  then  the  relevant  term  on  the  right  hand  side 
to  the  streamline  wave  was 


dv,  _  IT-  rj  —  (r,  +  r,)/2  /  AV 
aci/A) 2  A  '  l’ 


(2) 


After  that,  many  authors  more  directly  wrote  the  dv.  dzon  the 
right  hand  term  into  a  form  related  to  the  radius  of  curvature 
of  the  streamline.  Therefore,  this  type  of  method  is  also  commonly 
called  the  streamline  curvature  method,  or  it  is  called  the  velocity 
extension  method. 


The  second  type  is  the  flow  field  matrix  method  or  the 
so-called  stream  function  matrix  method.  It  Involves  the  introduc¬ 
tion  of  the  stream  function  into  the  continuity  equation  and  the 
stream  function  equation  is  substituted  into  the  equation  of 
motion.  Finally,  the  radial  equilibrium  equation  of  motion  expressed 
in  the  stream  function  form  is  obtained.  Then,  the  numerical 
differential  formula  is  used  to  write  the  stream  function 
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differential  equation  into  a  difference  equation  form  with  respect 
to  each  nodal  points  to  obtain  the  linear  algebraic  equation  set 
related  to  the  stream  function  of  the  nodal  points  of  the  entire 
field.  Use  the  numerical  computation  methods  of  the  linear 
algebraic  equation  set  such  as  the  direct  matrix  decomposition 
method  or  the  Iteration  method  to  solve  the  stream  function 
values  of  the  nodal  points  of  the  entire  flow  field.  Subsequently, 
other  physical  quantities  such  as  density,  velocity,  etc.  are  determined. 


II.  The  Line  Relaxation  Solution  Computer  Program  for  S-  Stream 
Surface  Stream  Function  by  Using  Non-orthogonal  Curvilinear 
Coordinates  and  Non-orthogonal  Velocity  Components 


In  the  past,  we  have  compiled  two  S,  computation  programs 

r  7  ^  8 1 

based  on  the  velocity  extension  method^'*  .  Due  to  the  slow 
convergence  inside  the  duct  durlne  computation,  this  more  satisfactory 
computing  program  was  subsequently  compiled.  With  regard  to  the 
details  and  explanation.  References  [7,8]  can  be  referred  to. 

1.  The  Stream  Function  Equation  to  be  Solved. 

The  dimensionless  form  of  the  S.,  stream  surface  radial  equll- 

^  r  4.— 8 1 

Ibriura  equation  under  arbitrary  curvilinear  coordinates  is 

■  v  »M  -0,)  ^ 

dx^dx’- 


where 


Vin,-  dl  RT^og  ds  \ 

-vl  I  -  --  SF-  +  S7‘  -  TTr  ^  a?) 

f  =  j  )  V  dj  I’,,;  )  1  0<p 


(M 

(5) 
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It  Is  briefly  noted  as 


-r-  .i,d-'ir,dx-dx'  +  ./  d->v  o{x  y  +  .■xw,dx-  +  .yd>v  dx‘  -  ^  (g) 

This  is  the  S2  stream  surface  radial  equilibrium  stream 
function  equation  to  be  used  in  seeking  the  solution. 


The  Difference  Equation  and  the  Iteration  Format. 


The  nvimerical  differentiation  is  expressed 
point  differential ■ formula  [9]  as 


using  the  three 


O'P'/ajrO,,*  = 

(a<r/dx‘)„*  -  C2\,V,-,_^  +  C22,V,.^  +  C2l,V,^,,^ 

id^v/dwy),.^  =  Dnyjr,.^.,  +  D\2yy„^  + 

0'y/d(r‘y),.^  =  D21,(P',_,,^  +  D22,'f„^  +  D23,'F,+,,t 
=  C21, Cl  +  C21,C12,'P’,-,.^ 

+  C21,C13*'P', -1.^+1  +  C22,C1 +  C22,C12^'f,.^ 
4-  C22,C\i^V„^+^  +  C23,Cll,<y, +  C23,C  1 2,'^,+,. * 
+  C23,C13*r^,.^^, 


/3O 


(7) 

(8) 

(9) 

(10) 


(11) 


By  substituting  equations  (7)- (11)  into  equation  (6),  the 
difference  equation  related  to  the  j ,k  point  is  obtained.  After 
simplification,  it  is  noted  as 

+  E]2^i-,.^  +  E13V,_,.^+,  +  E21'P‘,.^_,  +  E22>P^,.^  +  E23'P’„,+, 

+  E3Hr,+,.^_,  +  E32(r,+,.^  +  E33'P,+m+.  =  +  F,., 

With  regard  to  equation  (12),  it  is  possible  to  use  various 
methods  to  seek  the  solution.  One  is  the  point  by  point 
iteration  format 


-  EI2'F<'-,'i  -  EnrAl'.Ui 

-  E21(F):iL',  -  E23>^rtV■l^  “  -^3  llF'V.'.i-i  -  /^32r;;-,:i 

-  E33!r)Vi;U,)/F22, 
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Using  matrices  and  vectors,  it  can  be  expressed  as 

.  _  [  ,v/  -  /;][  Q4  ) 

where  [K]  is  the  coefficient  matrix,  [B]  is  the  vector  formed 
by  the  right  hand  terms,  and  [4']  is  the  stream  function  vector 
to  be  solved. 


The  second  method  is  to  list  all  the  difference  equations  with 
respect  to  all  the  j ,k  points  in  equation  (12).  Then,  an  element 
elimination  method  is  used  to  solve  the  linear  algebraic  equation 
set  related  to  t|/.  The  iteration  format  is 


Ell  +  E12'y^V,.^  +  +  £22'^';';' 


+  £2VP-<:U,  +  +  £32r,;',,^ 

= 


(15) 


By  using  matrices  and  vectors,  it  can  be  expressed  as 


With  regard  to  the  solution  of  equation  (16),  it  frequently 
uses  the  decomposition  of  the  coefficient  matrix  [M]  into  the 
lower  triangular  matrix  [L]  and  the  upper  triangular  matrix  [u], 
i.e.,  the  so-called  large  matrix  direct  decomposition  method 
is  used  to  find  the  solution. 


This  program  used  a  piece  by  piece  iteration  format,  and  equation 
(12)  is  briefly  noted  as 


(17) 


For  each  station  j ,  the  column  functions  form  k*2  to  k=n-l 
in  equation  (17)  together  with  the  two  functions  v.  ,  *0  and  y  *0 

J  >  1  J  J  “ 

are  used  to  obtain  using  the  chasing  method.  This  Iteration 

format  can  be  expressed  in  the  following  form  by  using  the  matrix 
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and  vector  form,  i.e.,  [M],  [<))],  and  [B]  in  equation  (16)  are  con¬ 
sidered  to  be  composed  of  some  submatrices  and  component  vectors. 
It  is  written  as 


t  Wu  1 

[m,i]  [m„]  [toj,] 

[ra,.,] 


(18) 


Then,  each  submatrix  block  and  component  vector  is  used  as 
a  unit  to  carry  out  the  iteration  computation.  After  the  component 
vectors  are  obtained,  it  is  relaxed  one  by  one,  i.e.,  by  the  line 
relaxation  or  the  block  iteration  method. 

It  is  actually  written  as 

+  [m, ■,,][!?•,]  +  [«,.,]['!>•;+,]  =  [B, ■]"-■)  (19) 

The  Iteration  format  used  is 

-  Cm, +,]"-•>.  (20) 

With  regard  to  solving  equation  (20),  the  similar  large 
matrix  direct  decomposition  method  is  used  to  find  the  solution. 
Because  the  matrix  [mj  2^  ^  diagonal  band  matrix  with  a  width 

of  three,  the  computation  at  this  time  is  much  simpler.  Its  matrix 
decomposition  and  the  solution  processes  can  be  written  into  an 
iterative  deduction  formula.  This  is  also  the  so-called  chasing 
solution  method.  This  is  a  computation  method  located  between 
large  matrix  direct  decomposition  and  point  by  point  iteration, 
which  accelerates  the  convergence  rate  and  does  not  occupy  a  lot 
of  interal  storage. 
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3.  Finding  the  Solution  for  the  Main  Equation,  Using  the  Chasing 
Method  to  Solve  the  Linear  Algebraic  Equation  Set. 

With  regard  to  equation  (17),  by  listing  all  the  equations 
of  the  station  j ,  the  following  equation  set  Is  obtained: 


'E,  •  =0 

Eh<p-,  4-  £2,<P-,  +  Ei,V,  =  B, 

(21-1)  , 
(21-2)  j 

1 

El^'E*.,  +  E2/E^  +  ES^fE^^x  = 

! 

(21-^)  1 

II 

■ 

1 

(21-/») 

Write  the  equations  (21-2)  to  (21-(n-2))  In  equation  (21)  one 
by  one  In  the  following  standardized  form 


'P'*  +  bnV 


♦+i  ~ 


(22) 


With  respect  to  (21-2),  we  have 

'P'2  +  iElJE2{)¥y  =  I}2/E2t  -  {El,  •  (?■,)/ E2j 


1.  e. 


b,=  El,/E2,  (2^) 

=  B2/E22  -  (£1,  •  (P-,)/E2,  (25) 


Substituting  the  standardized  form  'P’<-t  +  ^♦-i^P'*  =  f*-i  Into  equation 
(21-k),  we  get 


'P*  +  [£3*'(E2*  -  £l»**-i)]'P't  =  —  El^g^-l')/{E2^  —  E\^b,-x') 

Comparing  to  the  standardized  form  (22),  we  have 

b^  =  ( 26  ) 

s*  =  (27) 

co^  “  E2^  —  E\^b^.,  ^28) 


/31 
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According  to  equations  (24)-(28),  all  the  and  are 

obtained  first.  They  can  be  substituted  into  equation  (22)  to 
obtain  the  entire  values. 


4.  The  Computation  of  Other  Parameters. 

The  automatic  formation  of  the  mesh  nodal  points  and  the 
computation  of  the  geometric  parameters: 

In  order  to  reduce  the  input  of  original  data,  the  curvilinear 

mesh  was  formed  automatically  by  the  machine.  It  was  regulated 
2 

that  the  x  coordinate  lines  from  the  root  to  the  top  of  the  blade 

are  second  order  curves.  However,  the  coordinates  of  the  root  /32 

and  the  top  in  the  x^  direction  were  completely  given  in  order 

to  adapt  to  the  arbitrary  shape  of  the  meridian  duct.  The 

computation  of  the  mesh  was  given  in  terms  of  the  maximum  arch 

degree  DZ  and  the  distribution  coefficient  x  of  the  nodal  point 

along  the  blade  height: 

Ri  +  x^(R,-  RJ  (29) 

=  Z,  +  -  Z.)  +  DZ[1  -  4(x^  -  0,5)M  (30) 

The  maximum  arch  degree  DZ  was  given  at  the  front  and  rear 
fringes  of  the  blade.  As  for  the  other  stations,  it  was  calculated 
automatically  according  to  a  linear  distribution. 

After  the  nodal  coordinates  had  been  computed,  then  the  following 
geometric  parameters  were  calculated: 


01  =■  —  arctg  (22/rfr).‘.-,„„„  (  31 ) 

01  —  arctg  (32) 

{dx/dx^y  +  idr/dx'-y  (33) 

-  Vcde/dx')'  +  (dr/dx'Y  (3L) 
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The  dimensionless  computational  equations  for  the  physical 
quantities  are  given  in  the  following: 

The  Vgp  of  the  blade  spacing  stations  along  the  blade  height 
were  given  and  the  values  in  the  blade  duct  were  determined 
using  the  coefficient: 


-  (r,r).„]  (35) 

The  coefficient  is  usually  given  based  on  stream 

surface  computation 


-  A> 

(36> 

H,  +  A[(F,r),*,  -  n- 

1 

(37) 

1  =11  -  AK,r 

(38) 

for  a  moving  blade  A=l,  for  a  static  blade  A=0. 


K,  =  (aifvar-’)/(p/y^,) 

r,+,  -  T,  +  -  1)[(W  -  KV2U,  -  (II  -  V\  l),]/(KRT,^g.)  (^1) 

1'^=  F-:+  vi+  2F,F,sm(fl, -0,)  +  I’J  (42) 

=  p,(t,+,/7-, )■«*-"> .  (43) 

P-Tf,  (44) 

A;  =  ln(/'o/P'o)  (static  blade)  (''5) 

At  =  k/(k  -  1)  •  ln{(H,/H')/t;,(f/„7/'  -  1)  +  HI  (U6) 


for  a  moving  blade  during  compression  process. 

A/  =  k/(>:  -  1)  •  In  {(»,///')/[(//„. 7/',  -  O/jj  +  11)  (47  ) 

for  a  moving  blade  during  expansion  process. 

The  computational  formulas  related  to  the  output  printing 
are  as  follows: 
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total  temperature 


T„  -  li  ■  U-J.K  -  0/V-  •  R  ■  •  g,  ( 48 ) 


total  pressure 


pj  =  (r/7,)‘'<'-') 


•  p 


(49) 


Mach  number 


M  r,V  KTRT,^g,;uu 


(50) 


diffusion  factor 


£>  =  1  -  («'•,„/ •  It'in  •  (r,„  +  r„„,)] 


efficiency 


pressure  coefficient 


Meridian  inclination  angle  of  the  streamline 


<}>  coordinate  of  the  stream  surface 


(Pi  =  rp,_,  +  A/Kiy^./ii'*),-!  + 


where 


(51) 


/33 


^  =  ((/>.ou,/^i„)«‘““''‘  -  ll/(T«.,./T„i„  -  1)  (52) 


op  (/^Uoifr  P\n^/Cpa{n  pin')  (53) 


ff  =  arctg  (/?;)•  (180/»)  (54) 


A/  ==  \/  Ar^  +  As' 


(55) 
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(56) 


thickness  of  the  stream  slice  of  the  stream  surface 

r,  -  •  G  ■  ;P/[200(P  -  /):ryplF„,] 


the  flow  angle 


=  arctg 


(57a) 

(57b) 


5.  Brief  Description  of  the  Program. 

The  computation  program  was  compiled  using  the  ALGOL  language 
on  a  TQ-16  machine.  The  computation  procedures  of  the  program  are 
as  follows: 

1)  Input  the  original  data,  including  Incoming  flow  conditions, 
the  geometric  parameters  of  each  blade  section,  loss  coefficients, 
Vgj,  the  variational  coefficients  of  in  the  duct,  the  blade 
thickness  distribution  and  the  effect  coefficient,  the  flow  rate 
coefficient,  densities,  and  other  parameters. 

2)  Calculate  the  geometric  parameters  for  the  mesh  nodal 
points  according  to  equations  (29)-(34). 

3)  Make  them  dimensionless. 

4)  Form  the  initial  value  fields  of  4*  and  p  in  preparation 
for  iterative  computation. 

5)  Calculate  parameters  such  a.  W,j,>  H,  1,  s,  T,  and  p. 

6)  Find  the  value  of  V  using  the  chase  method  one  after  the 
other  along  the  J  stations.  After  line  relaxation,  the  hew  f 
distribution  of  the  entire  field  is  obtained. 
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7)  Calculate  and  ^ 

8)  Check  the  errors  of  f  and  p  .  If  the  accuracies  are  not 
enough,  then  go  back  to  step  5)  to  repeat  the  next  iterative 
computation. 

9)  If  the  computational  accuracy  is  sufficient,  then  print 
out  the  computed  results.  The  output  of  results  can  be  in  the 
dimensional  or  dimensionless  form  according  to  the  need.  It  Is 
also  possible  to  have  the  output  according  to  the  nodal  point 
position  or  the  streamline  position.  The  computational  accuracy 
can  be  given  according  to  the  need. 


III.  The  Comparison  of  Several  Numerical  Computational  Methods 
and  Discussion 

Although,  whether  a  design  is  successful  or  not  is  not  completely 
determined  by  the  computational  method  adopted  and,  to  a  great  extent, 
is  decided  by  the  design  parameters  and  experimental  data  chosen 
such  as  whether  the  selected  losses,  attack  angle  of  the  blade, 
and  lagging  angle  are  rational.  However,  different  computational 
methods,  due  to  the  different  internal  storage  capacity  requirements 
and  the  iterative  convergence  rates,  have  different  computational 
costs  for  nvimerical  computations.  Sometimes,  the  difference  can 
be  even  very  large. 

1.  The  Velocity  Extension  Method. 

Its  major  characteristic  is  to  alternatively  use  the  equation 
of  motion  and  the  continuity  equation  to  carry  out  the  iteration. 
Therefore,  the  efficiency  of  computation  is,  to  a  great  extent, 
determined  by  whether  this  mutual  iterative  correction  process  can 
converge  rapidly  or  not.  However,  this  type  of  mutual  iteration 
does  not  converge  rapidly  under  all  conditions.  Especially  when 


/3i 
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many  stations  are  established  in  the  blade  duct,  the  computation 
diverges  easily.  In  order  to  prevent  the  computational  process 
from  diverging,  therefore,  it  is  necessary  to  use  a  low  relaxation 
factor  a, which  is  less  than  1  when  calibrating  the  position  of 
the  streamline.  According  to  our  trail  computation  practice  and 
related  analysis  on  the  convergence  of  Iteratlon^"^^^ ,  under  the 
condition  that  other  parameters  are  unchanged,  the  magnitude  of 
the  mesh  Am  in  the  streamline  direction  is  a  major  factor  affecting 
the  convergence.  Furthermore,  there  is  a  quantitative  relationship 
that  the  relaxation  factor  a  is  proportional  to  Am^.  V/ith  x'-espect 
to  the  computation  of  a  two  level  fan,  when  the  aerodynamic  parameters 
are  Identical  and  only  the  mesh  densities  in  the  streamline  direction 
are  different,  there  are  the  following  effects:  a)  when  the 
computational  stations  are  only  taken  in  between  the  blade  gaps, 
the  relaxation  factor  can  be  chosen  as  0.3.  After  12  iterations , the 
accuracy  of  the  streamline  position  can  reach  0.98x10  ;  b)  when 

the  mesh  in  the  streamline  direction  becomes  denser,  each  blade 
zone  has  three  computational  stations  added  to  it.  That  is,  when 
one  computational  station  is  added  at  the  front  arid  rear  fringes 
as  well  as  in  the  middle  of  the  blade,  then  the  convergence 
becomes  apparently  poorer.  The  relaxation  factor  can  only  be 
chosen  as  0.04.  In  order  to  reach  the  same  accuracy  of  0.99x10 
the  number  of  iterations  is  increased  to  80;  c)  if  we  further 
densify  the  mesh  in  the  streamline  direction,  when  the  computational 
stations  in  the  blade  region  are  Increased  to  6  stations,  it 
becomes  divergent  easily  when  the  relaxation  factor  exceeds  0.01. 

After  273  iterations  using  a  0.01  relaxation  factor,  the  accuracy 

—2 

of  the  streamline  position  is  merely  0.58x10  .  Therefore,  we 

can  see  that  the  velocity  extension  method  (streamline  curvature 
method)  is  more  suited  for  the  aerodynamic  computation  in  the 
situation  of  blade  gap  station  of  a  large  mesh  spacing  in  the 
streamline  direction.  When  the  computation  enters  the  blade 
region  and  it  is  required  to  set  up  more  computational  stations, 
its  convergence  becomes  very  poor. 
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2.  Direct  Flow  Field  Matrix  Decomposition  Method. 

Its  greatest  special  feature  is  that  once  the  coefficient 

matrix  [M]  is  decomposed,  when  the  B  value  of  the  right  hand 

term  is  determined,  the  precise  solution  of  T  on  all  the  nodal  points 

can  be  directly  calculated  without  the  need  of  any  iterative  process. 

Therefore,  for  the  calculation  of  incompressible  fluids  and  Iso- 

energetlc  iso-enthalplc  axial  symmetric  flows  along  the  blade 

height,  it  is  possible  to  obtain  the  'P  value  of  the  entire 

field  by  one  computation.  The  computation  is  more  fast  and  accurate. 

However,  there  is  a  very  large  coefficient  matrix  [M]  in  the 

direct  matrix  decomposition  method.  The  storage  capacity  occupied 

by  the  coefficient  matrix  and  the  [L]  and  [u]  after  decomposition 

is  considerable.  It  is  proportional  to  the  number  of  stations 

along  the  streamline  direction  and  it  is  proportional  to  the 

square  of  the  number  of  stations  n  along  the  blade  height  direction, 

2 

i.e.,  «  mn  .  In  addition,  with  regard  to  the  usual  $2  stream  surface 
aerodynamic  calculation,  the  right  hand  term  B  of  equation  (16) 
must  still  be  iterated.  However,  at  this  time  we  can  use  a  very 
large  relaxation  factor.  It  is  possible  to  obtain  a  solution  with 
satisfactory  accuracy  after  a  few  Iterations. 

3.  The  Line  Relaxation  Iteration  Solution  of  the  Flow  Field 
Matrix. 

This  program  is  well  suited  for  the  three-dimensional  flow 
design  blades  when  it  is  necessary  to  set  up  more  computational 
stations  in  the  computation  of  the  S,  stream  surface  for  the 

r  111 

mutual  iteration  with  the  stream  surface*-  .  Under  this 
condition,  the  convergence  rate  of  the  velocity  extension- method 
is  too  slow.  The  typical  relaxation  factor  is  about  0.01.  In 
order  to  satisfy  the  0.1^  accuracy,  it  requires  about  200  iterations. 
The  typical  relaxation  factor  of  the  flow  field  matrix  line  relaxa¬ 
tion  program  is  1-1.65.  To  satisfy  the  0.155  accuracy,  15-20 
iterations  will  do.  It  is  10  times  faster  than  the  velocity 
extension  method. 
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The  flow  field  matrix  direct  decomposition  method,  because 
of  the  high  internal  ’storage  requirement,  such  as  when  m=47 , 
n=ll, the  nodal  point  number  is  517  and  the  internal  storage 
capacity  required  by  the  coefficient  matrix  is  about  28000  units, 
is  more  difficult  to  compute  on  a  Internal  storage 

computer.  However,  the  line  relaxation  method  can  allow  about 
1200  nodal  points  on  a  32K  internal  storage  computer  without 
using  external  storage. 

Therefore,  in  this  situation  the  line  relaxation  computation 
program  avoids  the  shortcomings  of  slow  convergence  and  high 
iteration  numbers  of  the  velocity  extension  method  and  it  also 
overcomes  the  difficulty  of  the  high  internal  storage  capacity 
requirement  of  the  direct  matrix  decomposition  method.  It  maintains 
the  advantages  of  the  low  internal  storage  capacity  requirement 
of  the  velocity  extension  method  and  the  fast  convergence  rate  of 
the  direct  matrix  decomposition  method.  It  is  Indeed  a  good  method 
which  is  worthwhile  recommmending. 

This  work  was  carried  out  under  the  direct  guidance  of  Professor 
Wu  Chung-Hua.  We  wish  to  express  our  deepest  gratitude. 
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COMPUTER  PROGRAM  FOR  DIRECT  PROBLEM  OF  FLOW  ALONG  S„ 
RELATIVE  STREAM  SURFACE  EMPLOYING  NON-ORTHOGONAL 
CURVILINEAR  COORDINATES  FOR  TURBOMACHINE* 


Ding  Yun-Feng  and  Deng  Zl-Tong 
(Computer  Center,  Academia  Sinica) 

ABSTRACT 

With  regard  to  the  direct  problem  of  S2  stream  sur¬ 
face  of  a  turbomachine  with  non-orthogonal  curvilinear 
coordinates  (subsonic  conditions),  a  convergence  algorithm 
was  realized  on  the  computer.  The  verifying  computation 
of  a  highly  loaded  axial-flow  two  stage  fan  showed  that  the 
computed  results  and  experimental  results  are  in  agreement. 

The  principal  deviation  equation  used  was  the  newly 
obtained  non-linear  form — equation  (8).  In  comparison  to 
the  conventionally  used  linear  form — equation  (13)  >  when 
the  high  subsonic  velocity  region  was  comparatively  large, 
the  new  algorithm  improved  the  convergence  of  the  numerical 
solution  process. 


I .  INTRODUCTION 

The  flow  inside  a  turbomachine  is  very  complicated.  The  three- 
dimensional  flow  theory  presented  by  Professor  Wu  Chung-Hua  in  the 
early  50 's  is  still  used  as  the  basis  for  the  computation  of  internal 
flow  motion  of  turbomachines. 

In  [1],  the  principal  equations  of  flow  along  the  stream  sur¬ 
face  and  S2  stream  surface  under  orthogonal  coordinates  were  derived 
in  linear  and  non-linear  forms.  In  the  early  60 's,  the  principal  equa¬ 
tions  of  the  motion  along  the  stream  surface  and  S2  stream  surface 
using  non-orthogonal  curvilinear  coordinates  were  derived  in  [2-J<]  in 
a  linear  form.  [1]  believed  that  it  is  more  convenient  to  use  the 
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linear  form  of  the  principal  equation  for  numerical  solution  for 
subsonic  conditions. 

Along  with  the  development  of  electronic  computers,  in  the  mid- 
60 's  the  three-dimensional  flow  theory  began  its  computer  programming 
realization  stage.  For  over  a  decade  up  to  the  present,  the  princi¬ 
pal  derivation  equation  has  been  the  linear  form  by  tradition. 

This  paper  derived  the  non-linear  form  principal  equation  of  the 
direct  problem  for  S2  stream  surface  using  non-orthogonal  curvilinear 
coordinates.  It  is  used  as  the  starting  equation  in  the  numerical 
solution.  The  difference  equation  of  the  stream  function  4  is  solved 
by  using  the  super  relaxation  iteration  method.  Computations  were 
also  made  with  the  principal  equations  in  a  linear  form.  A  comparison 
of  the  two  methods  showed  that  the  new  method  Improved  the  convergence 
better  in  the  high  subsonic  velocity  region.  In  the  meantime,  with 
regard  to  high  parameter  conditions,  measures  were  taken  to  gradually 
increase  the  parameters  and  to  correct  the  gas  entrance  angle. 

The  program  was  compiled  on  a  TQ-I6  computer  using  ALGOL-60 
language.  It  is  capable  of  analyzing  the  aerodynamic  properties  of 
turbomachines  (under  subsonic  conditions).  It  is  able  to  calculate 
the  entire  flow  field  of  a  two-stage  blade  at  once.  The  computed 
results  agree  well  with  the  experimental  results. 

II.  BASIC  EQUATIONS 

1  2 

Non-orthogonal  curvilinear  coordinates  (x  ,  x  )  were  used  on  the 

meridian  plane.  For  these  coordinates,  the  wheel  shell  line,  the 

wheel  boss  line  and  the  projections  of  the  front  and  rear  fringes  of 

the  blade  on  the  meridian  plane  can  all  be  selected  as  the  x^  line 
2 

or  the  X  line.  Therefore,  it  is  not  only  general  but  is  also  con- 

venlent  to  treat  the  boundary  conditions. 

This  paper  has  been  presented  at  the  2nd  National  Engineering  Thermal 
Physics  Conference  in  Hangchow  in  November  1978. 
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For  a  moving  blade,  let  us  choose  a  moving  coordinate  system 
i  2 

(x  ,  X  ,  which  rotates  at  an  angular  speed  u  with  the  moving  blade. 

1  2 

The  remaining  are  static  coordinates  (x  ,  x  ,  9).  Under  the  assump¬ 
tions  that  it  is  stable  and  adiabatic,  the  flow  on  the  stream  sur¬ 
face  satisfies  the  equations  [2-5]. 

Continuity  equation 


Equation  of  motion 
x^  direction: 


jv 

"J  “n 


Q/iW'  \/ ail)  + 


^  |_L[(}V'>sin(0,  -  di)  +  W^Waii]  -  — 


(1) 


(2) 


B(V„r) 

dx' 


'dx' 


X  direction: 


Energy  equation: 


,  Bcp  d{V„r) 
dx^  9x^  dt 

—  H^^sinC^i  —  0»))V ^\\]  -  [(T^^*sin(0i  -  9i) 

dx* 

^ an  f  ^^9  d  rJl  r)j 

+  U^')y an]  =  -jpr  [-  —  a?  ~  ^  a?  " 

^  =  0 
dt  ■ 


(3) 


(M 


i.e.,  I  is  a  constant  along  the  streamline  (relative  streamline  in 
the  dynamic  blade  region).  We  must  notice  that  I  is  not  continuous 
on  the  boundary  of  the  dynamic  and  static  regions.  The  gap  value  is 
u>rVf/gj,  ^  i.e.,  the  constant  has  a  different  value  in  each  region 
(because  i  ==  I!  —  QrVg/gj), 

Thermodynamic  relations:  j  ^  ^  j- ^  ^  Qi,x]/2g,j  (5) 

dt  =  id\nT  -  (/c  -  \)d\n!,)R/  J(.K  -  0  (6) 

where:  f  is  the  blade  thickness  effect  coefficient,  the  specific  beat 

ratio  k  =  k(T),  the  isobarlc  specific  heat 


f,  =  /?«/(«-  O/.  + 

([j/y  =,  (w-‘y  4-  (If'")'  +  2fFU>'‘sin(0,  -  di)  +  \V‘„ 

d\  =  cosdx  •  tgy, 


J41 


W'd„ 


(dynamic  coordinate  system) 
(static  coordinate  system) 


,  _  <yC V ,ir) 

~  'dx‘  dt 


ai  =  "J a,i  C'  —  1  >  -  )  > 


a.  /^\ 

dx‘^  dx^l 
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The  geometric  parameters  y  and  (}>  of  S2  stream  surface  are  not 
continuous  at  the  boundary  between  the  static  and  dynamic  regions. 


The  principal  derivation  equations:  From  equation  (1),  let  us 
define  the  stream  function  f : 


ay 

dx^ 

d'V 


=  IfdV'VaJK 


- fpW^y/on/l^ 

dx^ 


(7.1) 

(7.2) 


K  is  the  flow  rate  coefficient.  Using  equations  (5),  (6)  and  (7), 
it  is  possible  to  transform  equation  (3)  into 

^  (.-/'  -ha^-iwy-)  -  2[C8>  -  Ovy)  sin  (0,  - 

/7. 


dx'dx^ 


dW 

idxy 


M.-  -  (f^  -  51 1;  1^ 

—  ^  sin  (^1  ^i)  ^  ^  sin  (d\  ^2) 


+  ^  (1  +  1^!  -  —  (-’je  + 

^2  dx^  ax  dx^  \  T  }  dx^ 

il ('■“S-'-'s) 

{|r.  +  f  I? 


2  gsinC^t  -  5a 


a;t^ 


i  4.  ^XlA 

^2  ai  dx^> 


-  II:  (<•>'  +  +  ff:  -  ("'Wfi 

_  a  sin  (9,  -g,)  I£^ 

Ox'  aj  dx^  i(i  ax'i 


(8) 
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where 

A  =  4-  IK' sin  (0,  -di),  B  =  W'sinCd,  -  0j)  +  IK'(I  +  </{). 

Equation  (8)  Is  hyperbolic  form  when  W  >  a,  elliptical 
when  w  <  a. 

Equations  to  calculate  enthalpy  s  and  density  p:  The  viscous 
force  was  neglected  in  the  equation  of  motion.  However,  it  is  poss¬ 
ible  to  realize  part  of  the  viscosity  effect  through  the  calculation 
of  enthalpy  with  the  given  static  blade  total  pressure  loss  coeffi¬ 
cient  a  and  the  dynamic  blade  efficiency  n-  When  passing  through  the 
static  blade  and  dynamic  blade  (compressor),  the  enthalpy  Increments 
along  the  streamline  are 

A/*  =  — intr  static  blade  (9.1) 

=  — - — In - - -  dynamic  blade  (9.2) 

K-l  1  +  7j(H,/W,  -  1) 

The  subscripts,  1,  t,  represent  the  front  and  rear  fringes,  res¬ 
pectively.  S*  =  s  •  J/R.  The  enthalpy  is  calculated  according  to  the 
predetermined-  distribution  in  the  blade  array. 

The  computation  of  p  was  obtained  by  integrating  equation  (6) 
along  the  streamline.  We  obtained 

P  =  (10) 

Boundary  conditions:  The  inlet  gas  flow  parameters  are  given.  The 
wheel  shell  line  and  wheel  boss  line  are  streamlines.  The  outlet  con¬ 
duit  is  properly  extended  when  necessary  to  have  axial  gas  exit  to 
2 

make  W  =0.  Thus,  the  boundary  conditions  are:  inlet  'i'  =  I*  , 

oL 

wheel  boss  f  =  0,  wheel  shell  f  =  G,  outlet 


III.  THE  SUPER  RELAXATION  SOLUTION  METHOD  AND  CONVERGENCE 
ALGORITHM  AND  RESULTS 

If  we  assume  that  the  flow  in  the  gap  is  an  axial  symmetric  flow, 
a  unified  solution  can  be  sought  for  the  dynamic  blade  region  and  the 
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^3 


static  blade  region  simultaneously.  This  program  is  capable  of  cal¬ 
culating  seven  blade  sections,  including  the  upstream  region,  the 
guiding  blade,  the  two-stage  blade  and  the  downstream  region.  In  the 

direction,  a  three-point  unequal  distance  numerical  differential 

2 

equation  was  used.  In  the  x  direction,  it  was  equally  divided  accord 
ing  to  the  duct  height  of  the  loop  using  the  center  difference  equa¬ 
tion.  By  using  a  nine-point  difference  form,  the  difference  equa¬ 
tion  (8)  is  transformed  into  a  difference  equation. 


+  pio,.^  =  o 


(11) 


Using  the  super  relaxation  method  to  solve  the  difference  equa¬ 
tion  (11),  the  tf,.,  solved  from  equation  (11),  is  recorded  as 

(12) 


Then 


where  a  is  the  super  relaxation  factor,  n  is  the  number  of  times  to 
find  coefficients  Ql-QlO  and  p^  is  the  iteration  number  of  'V  under  the 
coefficients  which  are  found  in  the  nth  time.  Iterations  are  carried 
out  until 


The  solution  at  this  time  is  noted  as  which  represents 

the  solution  of  f  with  the  coefficients  Pi)l’— PloJl’  of  the  nth  time. 
This  is  the  iteration  of  f  which  is  called  the  outer  loop. 

Then,  the  and  ‘"Vn  obtained  are  used  to  calculate 

It'*,  tt'*.  If,  fs,  r.  Us*  ,  etc.  Equation  (10)  is  used  to  calculate 
and  Interated  until  <  «»  The  p  obtained  at  this  time 

is  called  .  Therefore,  new  coefficients  pij?**''— pioJl'^*'. 

are  obtained.  This  is  the  interation  of  p  which  is  called  the  inner 
loop . 


Repeat  the  inner  and  outer  loop  until 
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There  is  no  general  proof  on  the  convergence  of  the  afore¬ 
mentioned  f-p  iteration  process.  With  regard  to  the  S2  stream  sur¬ 
face  direct  flow  problem,  the  convergent  algorithm  under  high  para¬ 
meters  condition  is  an  important  problem  because: 

(1)  the  equation  is  non-linear,  as  the  M  number  approaches  1, 
the  degree  of  non-linearity  becomes  larger. 

(2)  the  tgY  term  in  the  coefficient  of  the  equation  has  a  dis¬ 
continuity.  Furthermore,  it  is  not  easy  to  provide  an  accurate  value. 

This  may  cause  the  Vq  obtained  to  be  discontinuous.  Very  frequently 
in  the  inner  loop  to  find  p  the  value  of  p  is  too  small  to  be  ration¬ 
al  and  T  appears  to  be  a  negative  value.  These  would  cause  the  iter¬ 
ations  to  fail. 

Specifically  with  respect  to  these  two  parts,  this  paper  adopted 
the  following  methods: 

1)  the  gradual  raising  of  the  computational  parameters.  That  is 
to  gradually  raise  the  inlet  velocity  or  gradually  increase  the  S2 
stream  surface  parameter  d^. 

2)  the  correction  of  the  S2  stream  surface  geometric  parameters. 

This  is  to  consider  that  the  gas  exit  angle  of  each  blade  array  is 
given  and  unchanged  and  the  gas  entrance  angle  is  determined  by  the 
gas  exit  condition  of  the  preceding  blade  (here  it  is  assumed  that 

r—?  .  When  the  model  is  built  on  a  cylindrical  surface,  6  is 

identical  to  y)-  The  condition  used  is  that  in  the  bladeless  regions. 

VgT  remains  a  constant  along  the  streamline.  The  tgy  in  the  blade 
array  varies  according  to  the  defined  regularity.  This  newly  defined 
tgY  is  used  to  recalculate  the  flow  on  the  corrected  S2  stream  surface 
until  the  difference  in  gas  entrance  angles  before  and  after  the  cor¬ 
rection  is  within  the  error  range  required. 

40 

These  measures  are  taken  according  to  the  need.  Usually,  they 
are  used  only  in  high  parameter  conditions. 
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Figure  1.  The  total  characteristics  of  a  two-stage  fan 
1 — efficiency;  2 — total  compression  ratio;  3 — experimental;  4 — cal¬ 
culated;  5 — revolving  speed,  designed  value  {%)]  6 — experimental; 

7 — calculated 


Figure  2.  The  radial  distribution  of  at  the  front  and 
rear  fringes  of  a  dynamic  blade 
1 — calculated;  2 — experimental;  3 —  (a)dynamic  blade  I; 

4 — calculated;  5 — experimental;  6 —  (b)djl.namic  blade  II 
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For  a  highly  loaded  two-stage  axial  flow  fan  [6],  the  working 
conditions  corresponding  to  velocity  cocJi“ ^  50  and  70  were  comnuted  . 
There  are  four  points  for  the  50%  working  condition  and  five  points 
for  the  70^  working  condition.  They  were  marked  as  50-1,  -2,  -3,  -4, 
and  70-1,  -2,  -3,  -4,  -5,  respectively.  In  the  computation,  48  sta¬ 
tions  were  setup  in  the  direction.  Seven  stations  were  set  up  in 
one  blade  section  in  which  four  stations  were  in  the  blade,  one  sta¬ 
tion  each  at  the  front  and  rear  fringes  and  one  station  in  the  gap. 

2 

11  stations  were  set  up  in  the  x  direction.  There  were  a  total  of 
528  nodal  points.  The  maximum  corresponding  M  number  for  the  70  work¬ 
ing  conditions  was  always  above  0.95.  The  maximum  relative  M  number 
of  70-1  was  0.9852.  When  8^  was  not  corrected,  the  numbers  of  inner 
and  outer  loop  Iterations  for  the  50  and  70  working  conditions  were  7 
and  25,  respectively.  The  iterations  of  'V  were  about  100  and  200 
times,  respectively.  The  computer  time  was  four  minutes  and  nine 
minutes,  respectively.  The  correction  of  6^  required  approximately 
10  minutes  of  computer  time. 


The  calculated  results  and  the  experimental  results  are  very  con¬ 
sistent.  Figure  1  is  the  total  characteristics  of  a  two-stace  fan. 
Figure  2  is  the  radial  distribution  of  the  relative  number  M  at  the 
front  and  rear  fringes  of  the  dynamic  blade  for  the  7-4  working  con¬ 
dition.  Figure  3  is  the  streamline  pattern. 

IV.  DISCUSSION  OF  THE  PRINCIPAL  EQUATION 


Substituting  equations  (7.1),  (7.2)  into  equation  (3)> 
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J  drV 
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}  =  0 

dx^' 

dt 

^  ‘  dx^ 

we  can  get 


(13) 


i.e.,  the  conditions  when  the  rotating  speeds  are  50%  and  70% 
of  the  designed  revolving  speed. 


Figure  3.  The  projections  of  the  streamline  on 
the  meridian  plane 

This  is  the  linear  form  principal  equation.  The  term 
T)pidx‘ it  =  \,l')  is  Included  which  inexplicitly  contains  second 
order  derivative  terms  of  Therefore,  equation  (13)  is  linear  only 
by  form.  In  the  low  subsonic  condition,  the  density  variation  is 
not  too  large.  It  is  actually  also  close  to  linearity.  The  linear 
form  recommended  by  [1]  is  indeed  convenient.  However,  in  high  sub¬ 
sonic  situations,  the  density  variation  is  very  large.  When  using 
equation  (13)  to  perform  Iterations  for  f,  9o/3r‘(<  =  i .  2)  remains 
unchanged.  However,  when  M  «= l  ,  the  variation  of  p  with  respect  to 
f  is  especially  sensitive.  Therefore,  it  is  not  favorable  for  the 
convergence  of  iterations  between  'F  -  p.  If  equation  (6)  is  used,  then 
the  variation  of  5^/9x'(;  =  l,  2)  with  f  is  taken  into  account  which  is 
helpful  for  the  convergence  of  p.  The  number  of  iterations  n  between 
'i'  -  p  is  reduced.  The  number  of  iterations  of  the  4'  field  itself 
could  also  be  decreased. 

In  a  loop  duct  example,  the  ratio  of  the  time  re¬ 
quired  for  computation  by  using  equations  (8)  and  (13)  0.4. 

At  this  time,  the  M  > 0.7  region  occupies  approximately  20%.  The  improve¬ 
ment  in  the  70-2  working  condition  was  not  obvious  because  most  of  the 
region  in  the  flow  field  had  very  low  M  number.  The  M  >  0.7  region 
merely  occupied  5'ii>  although  M  had  already  reached  0.96. 


The  new  solution  method  using  the  non-linearization  as  the  prin¬ 
cipal  equation  is  similarly  applicable  to  the  Inverse  problem  of 
stream  surface  and  S2  stream  surface.  Especially,  the  high  subsonic 
region  at  that  time  would  be  larger.  For  the  transsonic  conditions, 
this  new  method  should  have  special  meaning. 

The  authors  are  very  grateful  for  the  guidance  provided  by  Prof. 
Wu  Chung-Hua. 

APPENDIX  DERIVATION  OF  EQUATION  (8) 


Prom  equation  (6)  we  can  get 


9in/>  ^  Jg,  Sh  _  Ss* 

dx*  a*  Qx*  dx*  U  **  1 5  2 ) 
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Prom  equation  (5),  we  get 
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From  equation  (A. 2)  and  using  (A.l),  we  can  find 
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(  A .  2  ) 


(A. 3) 


(A. 4) 


(A. 5) 


From  (7.1),  we  get 

1  S'V''  _  JJ,,  j  Stnifa,)  ^  gj  Sh  _  St*  \ 

dx'  i  pa ,  dx '  d.v  ^  1  d.f  ‘  *  dx '  dx  *  J 

Substituting  equation  (A. 3)  Into  (A. 4),  we  get 

iifd  »  JlZ-  jl  +  - — - [(1  +  d\)W'  +  W  sm  (6,  -  ^ 

dx'  dx‘dx‘  I  a’  -  (IV)'  J  ).oa, 

W'  +  sin  (S,  -  9.)  3'V  _  J  Sln(Ja,)_  j.  A,  _  Os 
*  JIm]  (a*')’  i  ax'  a‘  »x‘  J 

Similarly,  we  can  get  the  same  kind  of  expressions  for  5H.‘,a'x',  3»"  ax 
In  which  the  expression  of  h^^  is 

*■  («."  *  -  ‘--y  &  -  ■£■) 

~  ^ —  if  'ir*  4.  [(\  +  d\)\w^y 

6x*  dx* 
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By  substituting  the  expressions  of  iw'idx‘,iw'idx‘  into  equation 

(3),  we  get  the  principal  equation  (8). 
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THEORY,  METHOD  AND  APPLICATION  OF  THREE  DIMENSIONAL  FLOW  DESIGN 
OF  TRANSSONIC  AXIAL  FLOW  COMPRESSOR 


Academia  Sinlca  and  Shenyang  Aeroengine  Company 
ABSTRACT 


On  the  basis  of  the  3-dimenslonal  flow  motion  theory 
of  Wu  Chung-Hua,  the  thermoaerodynamic  design  of  a  single 
stage  transsonic  axial  flow  research  purpose  compressor  has 
been  accomplished.  Through  five  cycles  of  iteration  of  the 
two  stream  surfaces,  the  center  streamline  on  the  S^  surface 
of  revolution  basically  coincided  in  shape  with  that  of  the 
corresponding  streamline  on  the  S2  stream  surface.  The  para¬ 
meters  on  it  such  as  VqI’,  As/R  and  stream  surface  thickness 
are  approximately  the  same.  A  solution  convergent  to  the 
three-dimensional  flow  field  is  obtained.  In  the  meantime, 
the  guidelines  to  judge  that  the  design  requirements  are 
satisfied,  the  required  flow  passage  needed  for  the  flow  rate 
and  the  blade  airfoil  are  determined.  The  flow  field  analysis 
with  respect  to  the  supersonic  gas  inlet  cross-section  shows 
that  due  to  the  stagnation  effect  of  the  duct  shockwave, 
quantities  such  as  V^r,  As/R  all  undergo  sudden  jumps. 

I.  INTRODUCTION 

In  the  early  50 *s,  reference  [1]  established  the  theory  for  the 
solution  of  subsonic  and  supersonic  3-dimensional  flow  motion  of  turbo¬ 
machines  using  two  types  of  relative  stream  surfaces.  Now,  it  has 
become  the  theoretical  basis  for  the  turbomachine  3-dlmensional  flow 
design  and  the  solution  of  the  3-dimensional  flow  field. 

For  an  axial  flow  compressor,  it  is  possible  to  neglect  the  rising 
curvature  of  the  stream  surface  and  to  approximately  assume  it  is 
a  surface  of  revolution.  The  design  problem  can  begin  with  the  cal¬ 
culation  of  the  Inverse  problem  of  the  center  S2  stream  surface  (S2jjj). 
After  going  through  shape  building  and  the  direct  problem  calculation 
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of  a  family  of  surfaces  of  revolution  procedures,  a  new  cycle  of 
computation  is  Just  beginning.  In  this  process,  based  on  cer¬ 
tain  Judgment  criteria  determined  by  experiments  or  experience,  cor¬ 
rections  are  made  for  parameters  to  make  the  aerodynamic  plane  more 
rational.  Therefore,  the  entire  iteration  process  is  not  only  the 
solution  process  of  the  S-dlmensional  flow  but  also  the  design  process 
of  the  3-dimenslonal  flow. 

In  a  transsonic  compressor,  the  airfoil  of  the  supersonic  gas 
entering  cross-section  and  the  calculation  of  the  revolving  surface 
flow  field  are  closely  connected. 

This  paper  introduces  the  3-dimensional  flow  design  of  a  single 
stage  transsonic  axial  flow  research  compressor.  The  major  design 
parameters  are:  flow  rate  61  kg/sec,  tangential  velocity  of  the  blade 
tip  of  the  rotor  400  m/sec,  wheel  boss  ratio  0.4,  stage  pressure  ratio 

I. 5  and  adiabatic  efficiency  85?.  Furthermore,  it  is  required  to  have 
sufficient  flutter  margin  and  distortion-resistant  capability.  In 
order  to  improve  the  varying  working  condition  characteristics,  a 
varying  curvature  inlet  guide  blade  is  used. 

II.  COMPUTATION  OF  THE  CENTER  S,  STREAM  SURFACE  AND  THE 
DESIGN  OF  THE  FLOW  PASSAGE 

1.  Simplification  of  radial  direction  computation 

Before  computing  because  the  radial  distributions  of  the 

flow  passage  and  the  work  done  are  not  determined, 

for  simplification,  the  radial  directloft  computation  is  carried  out 
first.  Seven  factors  such  as  the  stage  pressure  ratio,  the  law  of  work 
performed,  the  inner  and  outer  diameters  of  the  rotor  and  the  stator 
are  selected.  Three  levels  are  chosen  for  each  factor  and  an  ortho¬ 
gonal  design  of  flexible  level  is  used  according  to  the  table  beading 
^27^^  )  to  carry  out  27  computations.  The  results  are  as  follows: 


This  paper  has  been  presented  in  the  2nd  National  Engineering  Thermal 
Physics  Conference  in  Hangchow  in  November  1978. 
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(1)  Under  the  condition  of  =  i|00  m/sec,  1.5  pressure  ratio 
and  85%  efficiency  can  be  reached.  Therefore,  they  are  chosen  as 
the  design  pairameters. 

(li)  In  the  C  type,  equal  work  and  S  type  work  performing 
laws,  the  S type  not  only  adds  more  work  on  the  top  to  sufficiently 
utilize  the  work  performing  capability  of  the  supersonic  flow,  but 
also  greatly  reduces  the  stress  at  the  root  which  is  advantageous  to 
the  deflection,  diffusion  factor  and  M  number  of  the  flow  at  that 
place.  Consequently,  it  is  favorable  to  the  assurance  of  the  efficiency 
of  the  root  and  the  flutter  margin  and  distortion-resistance  of  the 
stage.  Therefore,  we  decided  to  adopt  the  S  type  work  law  (Figure  1). 

(ill)  The  optimum  value  ranges  of  the  inner  and  outer  diameters 
of  the  flow  passage  are  obtained  to  enable  the  computation  of  82^^^  to 
be  carried  out  in  a  more  concentrated  manner. 

2.  Characteristics  of  center  S2  stream  surface  computation 

It  completely  abandons  the  axial  symmetry  assumption  and  it  is 
also  different  from  the  principle  of  the  average  S2  stream  surface. 

(1)  Consideration  of  the  effect  of  blade  thickness.  On  the  sur¬ 
face  of  revolution  with  P  as  the  cascade  distance,  the  flow  rate  can 
be  expressed  as 

=  -  /)k 

From  this,  the  stream  surface  thickness  used  as  the  input  of  the  S2jjj 
computation  is 

~  (2»ror'‘VoI^’io)/(rL'>pkH'u)m 
-  (.2^n/0/(.P  -  t/P\ 

where  and  are  the  stream  surface  thicknesses  of  the  and 

$2  stream  surfaces,  respectively.  The  subscripts  0,  k  and  m  are  the 


Figure  1.  Distribution  of  Figure  2.  Distribution  of  (P-t)/ 

Vgr  along  the  blade  height  ?P  along  the  radial  direction 


parameters  in  the  gap,  in  the  passage  and  on  the  center  streamline, 
respectively.  From  the  above  equation,  we  can  see  that  the  physical 
meaning  of  c,  is  the  ratio  of  the  density  flow  on  the  center  stream¬ 
line  (pW^)m  and  the  average  density  flow  on  the  surface  of  revolution 
pW^ .  It  represents  the  effect  of  the  tangential  thickness  of  the 
blade  on  the  center  S2  stream  surface  flow  parameters  (in  the  mean 
stream  surface  computation,  ?  =  1.0). 

In  the  S2JJJ  computation,  ?  directly  takes  the  computational  results 
from  the  surface  of  revolution.  Finally,  the  radial  distribution 
of  (P-t)/?P  is  as  shown  in  Figure  2.  Figure  3  is  the  distribution 
regularity  of  the  M  number  computed  with  various  ^  values  with  respect 
to  S2JJJ  and  the  flow  angle  6.  We  can  see  the  magnitude  of  the  effect 
of  c  (in  the  figure  the  large  ?  is  the  aforementioned  calculated  value, 
the  small  ?  is  value  obtained  by  using  the  blade  front  fringe  radius  as 
the  tangential  thickness  at  the  front  fringe). 

(11)  The  work  law  in  the  passage  along  the  streamline 

-  [(F,r),  -  (F,r)J/((F.r),  -  (tV),] 

is  based  on  the  results  on  the  center  streamline  obtained  in  the 
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computation  of  the  surface  of  revolution  in  the  previous  cycle.  The 
typical  distribution  of  CVgr  on  the  subsonic  and  supersonic  cross- 
section  is  shown  in  Figure  4,  We  can  see  that  the  regularity  of  the 
subsonic  cross-section  can  approximately  be  considered  as  a  linear 
variation.  However,  the  work  performed  on  the  majority  of  the  super¬ 
sonic  cross-section  is  accomplished  by  shockwaves. 
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Figure  3*  The  effect  of  large 
and  small  ?  values  on  the 
results  2m 

1 — M  (large  ?);  2 — M  (small 
C);  5 — 6  (small  c);  4--6  (large 
?);  5 — z(mm) 


Figure  4.  The  variation  of  Vgr 
along  the  flow  passage 
1 — cross  section  II 


(ill)  The  enthalpy  Increase  law  in  the  blade  passage  along  the 


streamline  is 


*=  (sy  —  Sd^/Oj  —  Jo) 


In  the  subsonic  cross-section,  let  us  assume  a  linear  distribu¬ 
tion.  In  the  supersonic  part,  the  shape  of  shockwave  oosition  is  pre¬ 
sumed  and  the  enthalpy  Increase  is  calculated  on  the  middle  streamline 
(Figure  5). 

Because  the  projections  of  the  blade  gas  entrance  and  exit  bound¬ 
aries  on  the  meridian  plane  are  curves,  we  adopted  the  computer  program 
[61  compiled  according  to  C3>^»5l. 

3.  Selection  of  parameters  in  flow  passage  design 


55 


Figure  5.  Variation  of  enthalpy  along  the  flow  passage 
1 — cross-section  IV 

The  finally  determined  flow  passage  is  as  shown  in  Figure  6 
(part  of  the  streamlines  and  the  computation  stations  are  also  plotted 
on  it).  It  is  obtained  in  the  omputation  by  adjusting  the  curva¬ 
ture  of  the  flow  passage.  The  purpose  of  the  adjustment  is  to  change  47 

the  velocity  distribution  of  the  local  regions  in  the  root  and  at  the 
top  to  improve  the  elemental  characteristics  of  the  heavily  loaded 
rotor  top  and  roots  of  the  rotor  and  stator. 


-  ItMJ  6  lUU  MU 

Figure  6.  Flow  passage  pattern 
1 — (mm);  2 — cross-section  II;  3 — computation  station;  4 — z(mm) 

In  the  S2JJJ  computation,  we  combined  a  series  of  experimental  data 
and  directly  established  the  nco^p/la-D  curve  based  on  various  blade 
airfoils,  cross-sectional  positions  and  M  numbers  to  select  the  abso¬ 
lute  loss  value  for  various  positions.  Because  the  numerical  value 
of  the  D  factor  is  not  known  beforehand,  therefore,  it  requires  iter¬ 
ation  in  the  computation  of  the  loss  value.  However,  when  the  gas 
exciting  angle  turns  beyond  the  axial  direction,  the  value  of  the  D 
factor  can  no  longer  reflect  the  load  condition.  This  should  be 
noticed  in  using  the  loss  curve.  The  finally  selected  rotor  efficiency 
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and  stator  total  pressure  restoration  coefficient  are  given  in 
Figure  7. 


With  regard  to  the  flow  rate 
coefficient,  on  the  basis  of  analyz¬ 
ing  large  amounts  of  experimental 
data  according  to  the  condition  of 
this  research  compressor,  corres¬ 
ponding  numerical  values  are  chosen 
for  the  gap  stations.  In  the  blade 
array,  linear  interpolation  is  used. 
The  flow  rate  coefficient  values  of 
the  gap  stations  are  shown  in  Table 
1. 


4.  Judgment  guideline 


Figure  7.  Radial  distribu¬ 
tions  of  rotor  efficiency  and 
stator  total  pressure  restor¬ 
ation  coefficient 


In  this 
trolled: 

rotor 


design,  the  ranges  of  the  following  parameters  are  con- 

D,  <  0.4  A^h<4  5° 

A/h.  <0.7S  CF,,/F.,)  ^0.8S 


Stator 


O),  <  0.5  A«i,  45° 

5s  l.OO  (P,/P,)  >  1.03 


5.  Calculated  results 


The  convergent  results  of  computation  after  five  cycles  of 
iteration  are  shown  in  Figures  8-11  and  Table  2. 

In  order  to  facilitate  the  comparison  to  experiments,  we  also 
carried  out  S2jjj  computations  of  a  single  stage  and  single  rotor  with 
larger  gaps.  The  calculated  results  indicate  that  the  gas  entrance 
and  exit  angles,  the  inlet  M  numbers  and  the  diffusion  factors  of  the 
three  are  comparatively  close.  Thus,  as  long  as  we  carry  out  a  design 
for  a  single  stage  without  enlarging  the  gap,  we  should  be  able  to 
obtain  the  blade  and  wheel  which  satisfy  the  three  requirements  of 
the  experimental  piece. 
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TABLE  1 


far 

!  guiding  1 

f  rotor 

r  stator 

Stator 

upstream 

i  vane  j 

j  blade 

i  blade 

blade 

1  inlet  1 

inlet 

I  inlet 

exit 

1.000 

,  1.010  i 

_i _ : 

1.015 

!  1.025 

i  1.028 

1 _ 

1— far  forward  position;  2— inlet  of  the  guide  blade- 
the  dynamic  blade;  4— outlet  of  the  static  blade; 


3 — inlet 


of 


H? 


TABLE  2 
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inlet 

inlet 

exit 

axial 

static  1 

total 

total 

diff- 

M 

angle 

angle 

speed 

press.  ( 

press . 

temp . 

usion 

! 

1  ... 

1 

ratio 

!  rat  io  : 

ratio 

iratio 

factor 

tip 

1.340 

-64.85 

-58.22 

0.8905 

1.492 

1.159 

1.0396 

rotor 

mid-span 

1.115 

-56.49 

-42.97 

0.9060 

1.545  i 

1.144 

i  0 .432 

root 

0.6826 

-43.67 

0.40 

1.186  i 

1.454 

1.129 

0.290 

tip 

0.5798 

35.13 

0  i 

1 

1.043  ' 

0.9700 

1.000 

0.403 

stator 

mid-span 

1 

0.6529 

36.75  J 

0 

1.067  , 

0.9902 

1.000. 

0.368 

■root  J 

0.8010 

45.35 

0 

1.011 

1.214 
_ 1 

0.967c 

1.000 
i _ 

0.485 

jL; 

I  n>  I  ^  / 

/  / 

'  f 


U.6  .  /  H  ,  liLll 


iW 


Figure  8.  Radial  distributions 
of  the  inlet  and  outlet  M  nos. 

1 — stator  outlet;  2 — rotor  out¬ 
let;  3— rotor  inlet;  4 — stator 
inlet 


Figure  9.  Radial  distributions 
of  the  inlet/outlet  B  angle 
1— static  blade  inlet;  2~dynamlc 
blade  inlet;  3— dynamic  blade  out¬ 
let;  4 — B  (degree) 
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Figure  10.  Radial  distribution 
of  diffusion  factor 

1 —  dynamic  blade; 

2 —  static  blade 


Figure  11.  Total  pressure  dis¬ 
tribution  of  dynamic  blade  outlet 
1 — static  blade;  2 — dynamic 
blade;  3 — P°  (kg/cm^) 


III.  COMPUTATION  OF  SUBSONIC  SURFACES  OF  REVOLUTION 
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1.  The  computation  of  direct  surface  of  revolution  proble'^. 

Different  from  the  computation  of  surface  of  resolution  is 

carried  out  on  a  curved  surface  with  pre-determined  shape.  Its  main 
contents  are: 

(2) 

(1)  to  calculate  to  obtain  the  V^r,  AS/R  and  t'  on  the  middle 
streamline  as  the  input  data  of  S2jjj  computation  to  form  the  next  cycle 
of  computation. 

(ii)  to  compare  the  results  of  different  airfoils  and  different 
attack  angle  shapes,  to  choose,  adjust  and  correct  the  airfoil  to  mini¬ 
mize  flow  loss  to  the  extent  possible. 

The  selection  and  comparison  of  the  airfoil  are  based  on:  the 
small  pressure  gradient  of  the  rear  part  of  the  blade,  the  uniform 
velocity  reduction  of  the  flow  passage  and  the  low  peak  value  of  the 
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Figure  13.  Diagrammatic  example  of  the  blade  model  parameters 
1 — position  of  the  front  fringe;  2 — projection  of  the  blade  section 
on  the  meridian  plane;  3 — position  of  the  rear  fringe;  — vertex  of 
the  conical  surface;  5 — middle  arc  line  of  the  airfoil  on  the 
expansion  surface  of  the  conical  surface;  6 — axial  direction  position 

M  number  in  front  of  the  suction  surface.  Thus,  it  is  possible  to 
reduce  the  blade  surface  boundary  separation  loss  or  to  cause  the 
cascade  to  deviate  from  critical  over  flow  during  the  over  rotating 
work  condition  without  the  occurrence  of  plugging.  The  computation 
of  the  surfaces  of  revolution  adopts  the  computer  program  [7]  compiled 
based  on  [3j^,5].  The  typical  calculated  results  are  shown  in  Figure 
12. 


Through  the  computation  of  the  surfaces  of  revolution,  based  on 
the  parametric  range  in  this  design,  with  the  exception  that  the  top 
of  the  dynamic  blade  is  a  multi-arc  model,  the  other  sections  are  the 
double  circular  arc  blade  airfoil.  The  attack  angle  at  the  root  of  the 
dynamic  blade  is  chosen  to  be  5®. 
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2.  Geometric  model 
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Because  it  is  difficult  to  model  on  the  surface  of  revolution, 
the  geometric  modeling  is  carried  out  on  the  expansion  surface  of 
the  conical  surface  (Figure  13) •  On  the  expansion  surface  of  the 
conical  surface,  the  bending  angle  of  the  airfoil  is 

(ft  “  ^  —  I  —  E 

Where  E  is  the  angle  between  the  two  parent  lines  passing  through  the 
front  and  rear  fringes  on  the  expansion  surface.  The  E  angle  of  the 
positive  cone  is  positive. 


In  order  to  be  consistent  with  the  S2j^  computation,  in  this  design 
the  projection  length  of  the  airfoil  on  the  meridian  plane  is  fixed 
to  be  —  ,  For  a  certain  attack  angle  and  lag  angle,  the 

chord  length  b  is  determined  by  iteration.  In  the  meantime,  consider¬ 
ing  the  methods  to  obtain  the  experimental  cascade  data,  we  use  the  M 
number  and  6  angle  of  the  gap  station  in  the  calculation  as  the 

basis  for  the  model  on  the  conical  surface  formed  by  the  front  and 
rear  fringes.  The  lag  angle  computation  formula  used  is: 


where 


8 


(MsiL 

V  <r/m. 


+  A5 


m,  +  0.0002ft,  a  -  +  /?,) 

(N  is  the  number  of  blades),  A6  is  the  correction  value  of  the 
angle  which  is  chosen  by  experience  (Figure  lU). 


lag 


The  throat  choking  margin  is  an  important  Judgment  guideline  of 
geometric  modeling.  If  the  throat  flow  has  reached  the  critical  stage, 
then  it  is  possible  to  list  the  continuity  equation  from  the  inlet 
to  the  throat: 

Because  of  the  boundary  layer  and  three-dimensional  effects,  the 
actual  throat  area  should  be  greater  than  A^.  In  this  design,  the 
throat  aerodynamic  margin  of  the  various  sections  of  the  rotor  is 

.ft///.*  ^  l.Oi. 


61 


iMi  \ 


.  I 


-■=1  X 


'  \ 


\ 


\ 


3  t.' JJ.IV 


For  stators,  according  to  the  results  of 
cascade  experiments,  in  order  to  obtain 
the  minimum  flow  loss,  the  throat  area 
margin  A^/A^  is  about  1.03  within  the  M 
number  ran-^e  (0.6-0. 8)  in  this  design. 

The  entire  model  design  has  been  compiled 
into  a  computer  program.  The  details  are 
shown  in  [8]. 

3.  The  center  streamline  method  aero¬ 
dynamic  model 


f^igure  14.  Variation  of 
the  correction  amount  of 
the  lag  angles  A5 
1 — static  blade;  2 — 
dynamic  blade;  3 — cor¬ 
rection  of  the  lag  angle 
(degree ) 


For  the  subsonic  flow  on  the  surface 
of  revolution,  the  various  physical  quanti¬ 
ties  q  (1,  (t))  of  the  flow  field  can  be  ob¬ 
tained  from  the  value  <f*)  obtained  from 

computation  on  the  center  streamline 
through  a  Taylor  series  expansion  along  the 


tangential  direction.  For  a  higher  density  airfoil,  usually  suffi¬ 


cient  accuracy  is  reached  by  taking  three  terms: 


in  which  the  first  and  second  order  derivatives  are  obtained  from  the 
continuity  equation  and  the  momentum  equation.  Then,  from  the  flow 
rate  relation,  by  taking  the  +50%  streamlines  of  the  cascade  flow  rate, 
it  is  possible  to  obtain  the  upper  and  lower  airfoil  coordinates  of 
the  flow  passage  /  ^  /  *  U*.  <»»*)  .  Simultaneously,  it  is  possible  to 

obtain  the  pressure  distribution  of  the  surface  of  the  airfoil  and  the 
distribution  of  the  number  M. 


The  sectional  airfoil  of  the  rotor  root  obtained  from  computation 
using  the  computer  program  [9]  compiled  by  this  method  is  shown  in 
Figure  15.  In  the  figure,  the  surface  of  the  geometric  model  is  also 
plotted.  The  two  are  comparatively  close.  The  distribution  of  the  M 
number  of  the  surface  obtained  using  the  center  streamline  method 
model  is  shown  in  Figure  12. 
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_^ai  stre^-’.line  method 


double  arcs 


Figure  15.  Comparison  of  center  streamline  method  model 
to  geometric  model,  the  solid  line  is  the  center  streamline 
method  model,  the  dotted  line  is  the  double  circular  arc 
1 — center  streamline  method  model;  2 — double  circular  arc 

It  must  be  pointed  out  that  whether  the  model  surface  obtained 
by  the  center  streamline  method  model  is  rational  or  not  depends  on 
the  given  center  streamline  shape  and  the  parameter  distribution  on 
it;  the  latter  is  obtained  from  S2jjj  computation.  Therefore,  although 
the  center  streamline  method  model  does  not  explicitly  include  the 
attack  angle  and  the  lag  angle,  yet,  because  of  its  relation  with 
computation,  it  must  be  correlated  with  certain  characteristics  such 
as  the  attack  angle  and  lag  angle.  As  a  matter  of  fact,  only  after 
the  convergence  of  stream  surface  iteration,  can  it  then  be  compared 
to  the  corresponding  geometrical  model.  Owing  to  this  fact,  the  key 
to  this  method  is  S2jjj  computation.  Once  the  two  stream  surfaces  iter¬ 
ation  computation  converges  to  obtain  the  accurate  82^^^  data,  it  will 
be  possible  to  get  a  smoothly  surfaced,  reasonably  thick  airfoil  with¬ 
out  any  difficulty »  This  method  can  also  be  used  to  modify  the  airfoil. 


IV.  COMPUTATION  OF  SUPERSONIC  SURFACES  OF  REVOLUTION 
1 ,  Multi-arc  profile  of  sections  at  rotor  tip 

The  inlet  M  number  of  the  rotor  from  the  root  to  the  top  varies 
from  0.68  to  1.3^.  Therefore,  different  airfoils  are  adopted.  The 
sections  with  M  >  1.15  use  the  multiple  circular  arc  and  the  remaining 
sections  use  the  double  circular  arc. 


The  multiple  circular  arc  airfoil  is  plotted  on  the  conical  sur¬ 
face  expansion  surface.  Because  it  is  possible  to  adjust  the  bending 
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ratio  of  the  front  and  rear  segments,  the  position  of  the  connecting 
point  of  the  two  circular  arc  segments  and  the  maximum  thickness  posi¬ 
tion,  it  is  possible  to  control  the  wave  series  in  the  flow  passage  to 
a  certain  extent  to  reduce  shock  wave  losses. 

In  order  to  smoothly  connect  the  multiple  circular  arc  airfoil 
to  the  double  circular  arc  airfoil,  it  is  possible  to  properly  choose 
the  forward  segment  bending  ratio  and  the  forward  segment  chord  length 
of  the  multiple  circular  arc  airfoil  near  the  connecting  point  of  the 
two  airfoils  to  make  the  two  circular  arc  segments  approach  one  circu¬ 
lar  arc. 

For  a  cascade  airfoil  with  supersonic  gas  inlet,  it  is  necessary 
to  satisfy  the  only  gas  entrance  angle  relation.  Therefore,  it  is 
necessary  to  carry  out  Iteration  in  the  modeling  process. 

2.  Computation  at  the  inlet  region 

Both  experimental  results  and  theoretical  analysis  proved  that 
there  exists  a  stable  periodic  shock  wave  and  expansion  wave  series 
in  the  head  regions  of  a  two-dimensional  supersonic  gas  Intake  cascade. 
They  have  the  function  of  adjusting  the  gas  Intake  condition.  Thus, 
corresponding  to  a  certain  gas  intake  M^,  there  exists  the  only  gas 
Intake  angle  3^ •  The  correlation  can  be  obtained  through  the  solution 
of  the  basic  equation  established  by  the  inlet  flow  field  between  the 
upstream  and  inlet  regions. 

The  actual  method  is:  to  use  the  detached  shock  wave  model  [10] 
of  a  two-dimensional  isolated  wing.  It  is  placed  on  the  conical  sur¬ 
face  expansion  surface  as  the  shock  wave  ahead  of  the  cascade.  With 
regard  to  the  supersonic  flow  field,  the  characteristic  line  equation 
in  [11]  and  the  correlation  are  used  to  find  the  solution.  The  com¬ 
puter  program  is  shown  in  [12].  The  computer  program  starts  from  the 
semi-infinite  cascade  of  a  blade  with  head  and  then  propagates  from 
upstream  to  downstream.  After  passing  through  3-^  cascade  distances, 
the  M  number  of  the  corresponding  position  in  the  flow  field  is  equal 
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Figure  16.  The  only  gas  Intake  Figure  17.  Shock  wave 

angle  relation  of  supersonic  gas  track  of  the  dynamic  blade 
intake  cascade  1 — shock  wave  track 

1 — 6^  (degree) 

to  the  value  of  3;  the  periodic  condition  of  the  cascade  is  satisfied. 
Then,  the  invariance  relation  between  the  incoming  flow  and  the  inlet 
parameters  is  established  for  a  control  system  whose  tangential  dir¬ 
ection  is  one  cascade  distance  in  order  to  obtain  the  M  -3  relation 

o  o 

of  a  uniform  incoming  flow.  The  typical  result  is  shown  in  Figure  16. 

When  the  number  of  the  Incoming  flow  is  not  too  large  and  the 
taper  of  the  outer  wall  of  the  machine  box  is  not  too  high,  and  con¬ 
sequently,  the  variations  of  stream  surface  thickness  and  radius  are 
not  large,  the  supersonic  inlet  region  can  be  approximately  considered 
as  single  wave  flow.  Consequently,  the  Prandtl-Meyer  equation  is 
obtained.  Its  advantage  is  that  the  parametric  relation  between  the 
inlet  region  and  the  upstream  area  can  be  directly  established.  It 
is  only  necessary  to  analyze  one  flow  passage  of  an  Infinite  cascade 
in  order  to  obtain  the  upstream  condition.  Thus,  the  work  load  is 
greatly  reduced. 

Under  these  conditions,  with  the  assltance  of  the  neutral  point 
and  neutral  wave  concepts,  from  the  continuity  equation 

f»PoW' o/  cos^o  —  T,p,W,b  sin  fi. 
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(where  b  is  the  width  of  the  first  passage  swallowing  a  Mach  wave  in 
the  cascade  passage)  and  the  P-M  correlation 

(where  , - ; -  , - - — -  ^ _ 

i'  =  V  (^  +  i)/(^  —  O  tg"‘v  —  i)/(^  +  i)(.v/-  —  1)  —  tg"'\/.v/'  —  0 

it  is  very  easy  to  obtain  the  two  unknowns  M  and  6  . 

o  o 

With  regard  to  the  effect  of  three-dimensional  effect  on  the 
computation  of  the  inlet  region,  typical  computed  data  shows  that  when 
the  stream  surface  thickness  differs  by  1%,  the  calculated  attack  angle 
deviates  by  1°  and  the  throat  choking  margin  will  be  off  by  1% . 

3-  Estimation  of  the  Vgr,  As/R  and  c  values  on  the  center  streamline 

Due  to  the  difficulty  in  accurately  determining  the  position  and 
shape  of  the  shock  wave  in  the  passage,  a  model  [13]  with  a  certain 
accuracy  as  compared  to  the  actual  experimental  data  is  used  to  esti¬ 
mate  shock  wave  loss.  Thus,  it  is  approximately  considered  that:  in 
front  of  passage  shock  wave,  it  is  iso-enthalplc  expansion.  Through 
the  shock  wave  enthalpy  Increment  jump,  the  shock  wave  rear  streamline 
airfoil  loss  Increases  linearly  until  the  outlet  (Figure  5).  The 
corresponding  meridian  plane  shock  wave  track  line  is  shown  in  Figure 
17. 


From  the  calculated  results  of  V^r,  we  can  see  that  on  the  super¬ 
sonic  gas  Intake  cross-section,  shock  wave  work  occupies  a  great  pro¬ 
portion.  Along  with  the  increases  in  the  gas  Intake  M  number,  this 
proportion  also  Increases.  Especially  for  the  sections  near  the  top, 
the  work  done  is  almost  completely  accomplished  by  shock  wave.  There 
is  almost  no  increase  in  the  passage  behind  the  wave.  Even  accelera¬ 
tion  would  occur.  In  view  of  the  Important  effect  of  the  design  of 
the  sections  near  the  top  on  the  entire  stage,  it  is  necessary  to 
carefully  select  the  blade  number,  chord  length,  flow  turning  angle  and 
other  geometric  parameters  of  the  airfoil  to  ensure  the  design  require¬ 
ments  . 
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Figure  18.  Center  streamlines  on  the  surface  of 
revolution  during  iteration  convergence  and  on  82^^^ 

1 — l(m);  2 — 5th  cycle  calculated  result  of  32^^^;  3 — ^th  cycle  cal¬ 
culated  result  of  S2JJJ;  ^ calculated  result;  5 — 4>(radian) 

V.  THE  FINAL  RESULTS  OP  STREAM  SURFACE  ITERATION  AND  CUMULATIVE 
ITERATION 


After  five  cycles  of  Iteration,  convergent  result  is  obtained. 
The  major  signs  are  that:  The  center  streamline  on  the  surface  of 
revolution  basically  coincides  with  the  corresponding  streamline  on 
S2jjj  in  shape  (Figure  I8).  The 
parameters,  such  as  Vgr,  AS/R 
and  thickness  distribution  are 
approximately  the  same  on  them. 

At  this  time,  the  Judgment  cri¬ 
teria  for  each  term  have  been 
met . 


The  airfoils  of  all 
the  stream  surfaces  place  their 
centers  of  gravity  on  the  same 
radial  line  to  carry  out  cumulative 


Figure  19.  Schematic  diagram 
of  airfoil  cumulative  iteration 

iteration  [1^].  In  addition,  the' 


meridian  projection  of  the  chord  length  and  the  density  value  are 
adjusted  to  obtain  the  entire  blade.  The  special  feature  of  the 
cumulated  Iteration  in  this  design  is  that  the  weight  (represented  by 
the  thickness)  of  each  stream  surface  is  considered  when  the  position 
of  center  of  gravity  is  calculated.  Consequently,  the  cumulative 
iteration  method  is  consistent  with  the  stream  surface  theory  (Figure 
19).  67 


After  airfoil  cumulated  iteration,  we  must  also  calculate 
elastic  distortion  of  the  blade  and  carry  out  correction  for  the 
blade  coordinate. 

VI.  CONCLUSIONS 

The  two-stream  surfaces  theory  was  originally  proposed  for 
entirely  subsonic  or  entirely  supersonic  flow  motion.  In  recent  years, 
it  has  been  widely  applied  and  verified.  Now,  we  used  it  to  design 
a  transsonic  compressor,  l.e.,  through  the  iterative  computation 
between  a  center  S2  stream  surface  and  a  family  of  S^  surfaces  of 
revolution  to  carry  out  the  thermoaerodynamic  computation  as  well  as 
to  select  and  modify  the  airfoil.  It  should  be  said  that  this  was 
an  attempt.  We  hope  that  through  the  repeated  process  of  experiments 
and  modifications  later  to  continually  perfect  and  develop  the  three- 
dimensional  flow  motion  theory  to  gradually  establish  a  more  perfect 
three-dimensional  flow  design  system  of  a  transsonic  turbomachlne . 

The  supersonic  inlet  section,  due  to  the  presence  of  shock  waves, 
makes  the  work  done  and  enthalpy  increase  rapidly.  Especially  for 
the  sections  near  the  top,  acceleration  may  occur  even  in  the  passage 
behind  the  shock  wave.  Therefore,  to  accurately  control  the  wave 
series  is  an  Important  problem  in  the  design  of  a  transsonic  compressor. 

The  entire  work  was  carried  out  under  the  guidance  and  participa¬ 
tion  of  Prof.  Wu  Chung-Hua.  Those  who  participated  at  various  stages 
are:  Xu  Jian-zhong,  Lu  Wen-qian,  Jiang  Zuo-ren,  Zhao  Llan-xlng,  Wu 

Pang-xlm,  Tu  Zhuan-She,  Yang  Jin-sheng,  Bai  Yln-ming,  Zhang  Yin  and  Lu 
Gen-shou  of  Academia  Slnlca;  Zhou  Wei-guo,  Zhang  Bao-xiang,  Gao  Yao- 
lin,  Pau  Yun-hu,  Xu  Bao-mlng,  He  Shi-hui,  Wu  Ke,  Zhang  Zhi-fang  and 
Yu  Jln-shl  of  the  Third  Mechanical  Department  of  Shenyang  Aeroengine 
Company . 
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EXPERIMENT  OF  TWO-DIMENSIONAL  TRANSONIC  TURBINE  CASCADES* 


Ai  Kiao-Yi** 


(Beijing  Heavy  Electric  Machinery  Plant) 


ABSTRACT 


Based  on  the  experimental  results  of  three  cascades,  this 
paper  analyzes  the  generation  and  development  of  the  shock 
flow  patterns  in  the  transonic  turbine  cascade  channel  and 
their  influence  on  the  cascade  aerodynamic  performance. 

The  important  characteristics  of  the  limiting  regime, 
limiting  loading  and  limiting  inlet  angle  of  the  cascades 
are  discussed.  It  is  apparent  that  the  transonic  turbine 
cascade  has  the  advantage  of  high  loading  capacity  and  high 
efficiency . 


1.  INTRODUCTION 


Two-dimensional  cascade  experiments  are  an  essential  means  to 
study  the  aerodynamic  performance  of  transonic  turbines.  Wind 
tunnel  tests  have  been  conducted  for  the  850  mm  rotor  blades  on  our 
10,000  kilowatt  transonic  turbine  with  three  different  blade  pro¬ 
files  at  various  angles  of  attack.  Figure  1  shows  the  three 
airfoil  profiles  used  in  the  tests. 

The  experiments  are  conducted  on  the  transonic  cascade  wind 
tunnel  at  the  Chlng-he  Laboratory  of  our  plant.  The  cross  section 
of  the  test  section  is  400  x  100.  The  Mach  number  M2^  at  the  exit 
of  the  cascade  can  reach  about  1.7.  The  duration  for  steady  operation 


*  This  paper  was  presented  at  the  Second  National  Convention  of 
Engineering  Thermophysics,  November  1978. 

**  Other  contributors  to  this  study:  Zhang  Ru-muo,  Min  Da-fu, 
Kong  Qlng-Xi  and  Yang  De . 
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Figure  1.  Airfoil  profiles  of  the  blades. 
1 — Cascade 


is  50  -  100  seconds.  The  wind  tunnel  system  satisfies  the  require¬ 
ments  of  cascade  testing  after  examination.  The  error  of  pressure 
measurements  is  0.4%  and  the  error  in  loss  coefficient  measurements 
is  ±0.  37o  .  The  exit  angle  is  measured  by  self-alignment  probes. 

Its  error  is  ±15'.  During  the  experiment,  the  total  pressure  and 
the  static  pressure  both  upstream  and  downstream  of  the  cascade  and 
the  surface  pressure  along  the  blades  are  all  recorded  electronically. 
The  wave  spectrum  in  the  cascade  channel  is  also  recorded  by  means 
of  a  505-1  shadowgraph  camera  with  a  200  mm-diameter  lens. 


The  condition  of  supersonic  flow  at  the  exit  of  the  cascade 
must  be  considered  when  taking  pressure  data.  The  total  pressure  Pg^ 
measured  by  the  probe  is  actually  the  total  pressure  behind  the 
normal  shock  at  the  tip  of  the  probe,  while  the  static  pressure 
P2  on  the  side  wall  downstream  of  the  cascade  is  the  static  pressure 
upstream  of  the  shock.  Hence,  upon  introducing  the  normal  shock 
equations : 

1 

K-l 


61  _  (£+1  (-2^ M!  - 

P.  \  2  /  XJC+l  K+l/ 


and 


Pjj 

Pm 


K  - 
K  +  1/ 
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the  total  pressure  ?q2  downstream  of  the  cascade  without  probe 
disturbance  can  be  found  readily.  The  resulting  loss  coefficient 
calculated  eliminates  the  disturbance  due  to  the  probe. 

All  the  experimental  data  are  arranged  in  the  following 
parameters : 

Isentropic  Mach  number  downstream  of  the  cascade 

Mu  =  J—^  -  n 

y  K  -  \ 

Approaching  Mach  number  upstream  of  the  cascade 

w,  =■  [(Poi/Pic,)^  -  n 

V  ^  1 

Mean  loss  coefficient 

^  (P.i//’«)  *  -  1 

Mean  exit  angle 

?icf  “  "Z"  ( 

S  Jo 

Isentropic  Mach  number  on  the  blade  surface 

M.,  -  {iPjPif^  -  1] 

Loading  coefficient 

C,  ”  M\cot§\  +  M*  cot^i 

Flow  rate  coefficient 
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where  ~  sm (//5)  effective  exit  angle  of  the 

cascade.  This  paper  attempts  to  perform  a  preliminary  analysis 
based  on  the  experimental  data. 

2.  GENERATION  AND  DEVELOPMENT  OP  SHOCK  WAVE  PATTERN  IN  THE 

CASCADE  CHANNEL 

Shadowgraphs  of  the  II#  cascade  channel  at  the  design  inlet 
angle  for  different  regimes  are  taken  and  are  illustrated  on  the 
next  page.  Figure  2  shows  the  isentropic  Mach  number  distribution 
along  the  blade  surface  of  the  cascade.  For  convenience,  the  wave  /57 
system  in  the  channel  is  reconstructed  in  Figure  3  according  to  the 
shadowgraph . 


Figure  3.  Wave  pattern  in 
channel . 

1 —  Shock  wave;  2 —  Expansion  wave; 
3 —  Wake 


Figure  2.  Mach  number  distribu¬ 
tion  along  blade  of  II*  cascade 
at  design  inlet  angle. 


From  Figure  2,  it  can  be  noticed  that,  as  the  regime  Mach 
number  M^^  >  0.7,  a  local  supersonic  region  appears  on  the  portion 
of  the  blade  surface  where  the  curvature  is  the  largest.  A  fan  of 
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relatively  weak  A  shock  waves  from  the  local  supersonic  region  can 
be  seen  from  shadowgraph  a.  They  are  represented  by  AB  in  Figure  3. 
The  first  peak  value  of  the  velocity  distribution  curve  in  Figure  2 
gives  the  position  of  the  generation  of  wave  AB.  As  exceeds  a 

certain  value,  the  pressure  at  the  exit  of  the  cascade  drops  drasti¬ 
cally  due  to  the  sudden  circulating  flow  around  the  trailing  edge 
of  the  blade.  Thus  supersonic  region  may  also  occur  near  the 
trailing  edge,  especially  on  the  upper  surface.  The  flows  along 
the  upper  and  the  lower  surfaces  of  the  blade  meet  at  a  certain 
distance  away  from  the  trailing  edge  and  the  two  fans  of  the  com¬ 
pression  waves  reinforce  into  a  pair  of  oblique  shocks,  that  is,  the 
left  running  wave  FH  and  the  right  running  wave  FG  (Figure  3). 
Together,  they  are  called  the  sparrow  tall  waves  of  the  blade. 

As  ^2t  Increased,  the  local  supersonic  region  enlarges 
accordingly.  When  sonic  line  spans  across  the  entire 

cascade  channel,  causing  the  flow  to  be  choked.  For  such  a  condi¬ 
tion,  the  right-running  wave  PG  will  meet  wave  AB  (shadowgraph  b). 
Since  the  flow  behind  wave  AB  is  subsonic  due  to  a  sudden  increase 
of  pressure  accompanied  by  a  sudden  drop  in  velocity  across  the 
wave  (Figure  2),  wave  PG  ceases  to  propagate  after  crossing  wave 
AB.  On  the  other  hand,  as  the  flow  along  the  lower  surface  of 
the  blade  passes  the  trailing  edge  C,  it  experiences  a  sudden  drop 
of  pressure  and  accelerates  around  the  edge  rapidly.  An  expansion 
wave  fan  CDE  centered  at  C  is  then  generated  between  the  sonic 
line  and  the  tall  wave,  and  reflected  expansion  waves  are 
generated  from  the  upper  surface  of  the  adjacent  blade.  A 
further  Increase  of  the  flow  regime  would  result  in  a  corresponding 
increase  of  the  strength  and  the  range  of  the  expansion  wave  fan. 

At  the  same  time,  the  wave  AB  generated  from  the  upper  surface  of 
the  adjacent  blade  passes  further  downstream,  and  its  strength  is 
reduced  under  the  interference  of  the  expansion  wave  fan  CDE.  Prom 
shadowgraph  c,  the  disturbed  wave  AB  in  the  channel  is  not  very 
clear  due  to  the  poor  transmissivity  of  organic  glass.  But  it  is 
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still  very  sharp  in  the  region  downstream  of  the  cascade.  The 
velocity  behind  wave  AB  is  increased  correspondingly,  such  that  the 
entire  rear  section  of  the  cascade  channel  is  supersonic  (Figure  2). 
The  right-running  tail  wave  PG  will  now  pierce  the  entire  channel 
span  to  reach  the  upper  surface  of  the  adjacent  blade.  Under  the 
Influence  of  the  boundary  layer,  a  shock  wave  GI  with  complicated 
structure  and  large  dimensions  is  reflected.  The  wake  behind  the 
cascade  is  a  parallel  flow  of  varying  velocity.  Hence  when  wave 
AB  and  the  reflected  wave  GI  come  across  the  wave,  they  are 
deflected  and  partially  reflected  (shadowgraph  c). 

With  the  increase  of  intensities  of  all  the  wave 

systems  increase  and  the  Mach  angles  decrease  correspondingly, 
causing  the  incident  point  G  of  the  right-running  tail  wave  on  the 
upper  blade  surface  to  displace  downstream  (shadowgraphs  c  and  d). 
The  second  peak  value  of  the  velocity  distribution  curve  in 
Figure  2  indicates  the  location  of  the  incident  point  G.  If  is 
apparent  that  it  displaces  downstream  gradually  as  M2^is  increased. 
When  >  1-^,  the  incident  point  reaches  the  exit  trailing  edge, 

and  the  reflected  wave  aligns  with  the  left-running  tail  wave 
(shadowgraph  e).  A  further  increase  of  M  would  cause  the  right- 
running  tall  wave  PG  to  depart  from  the  cascade  channel.  Due  to 
the  non-uniform  flow  behind  the  cascade,  it  is  deformed  into  a 
curve  with  its  wave  surface  curving  downstream  (shadowgraph  f). 

The  mutual  interference  of  the  shock  wave  and  the  wake  will  become 
the  decisive  factor  of  the  flow  regime. 

3-  CRITICAL  MACH  NUMBER,  CHOKING  MACH  NUMBER  AMD  FLOW  RATS 

COEFFICIENT 

The  critical  Mach  number  is  the  operating  Mach  number  at 

which  locj.l  supersonic  regions  appear  in  the  cascade  channel.  Based 
on  the  tested  data  of  the  velocity  distribution  along  the  blade 
surface  for  various  operating  conditions,  the  relation  between  the 
maximum  Mach  number  on  the  blade  surface  and  the  regime  Mach  number 
can  be  determined  (Figure  4).  The  critical  Mach  number  can  then  be 
obtained  from  such  curves.  It  is  found  that  for  I*  cascade:  = 

*  0,76,  and  for  II*  cascade:  W2tKP  * 


/ 
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V/hen  the  cascade  is  operating  in  the  transonic  range  and  the 
sonic  line  spans  across  the  entire  channel,  the  flow  rate  will  no 
longer  increase  with  that  is,  the  channel  is  choked.  The 

approaching  Mach  number  upstream  of  the  cascade  also  does  not 
increase  and  is  termed  the  choking  Mach  number  depends 

on  the  geometric  dimensions  of  the  airfoils  and  the  cascade,  the 
inlet  angle,  and  the  viscosity  of  the  working  fluid.  Figure  5 
shows  the  variation  of  the  approaching  Mach  number  M^  with  the 
regime  Mach  number  measured  for  the  I*  cascade.  The  choking  Mach 
number  can  be  read  off  directly.  Figure  6  shows  the  measured 
variation  of  the  choking  Mach  number  versus  different  inlet  angles 
for  the  III  cascade.  The  dotted  curve  is  the  calculated  curve 
based  on  one-dimensional  isentropic  theory  , assuming  that  the  sonic 
line  passes  through  the  geometric  throat  of  the  channel.  The  two 
curves  are  close  to  each  other  at  low  inlet  angles. 


Figure  4.  Relation  between  maximum  Mach  number  M 
on  the  blade  surface  and  the  regime  Mach  number  M 


itmax 
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The  flow  rate  coefficient  y  is  defined  as  the  ratio  of  the 
actual  choked  flow  rate  to  the  theoretical  design  flow  rate.  In 
general,  the  theoretical  flow  rate  is  calculated  by  using  the 
geometric  throat  area  of  the  channel,  while  the  actual  flow  rate 
depends  on  the  aerodynamic  throat  area.  In  the  experiment,  it 
can  be  computed  according  to  one-dimensional  isentropic  theory 
by  employing  the  measured  approaching  Hach  number  and  the  inlet 
cross  section.  Figure  6  also  shows  the  variation  or  the  flow-rate 
coefficient  with  the  inlet  angle  for  the  III*  cascade.  Although 
''^Imax  coi  i’espondingly  as  S-j^  increases,  this  does  not 

mean  that  the  throughflow  capacity  of  the  cascade  would  increase. 
It  is  because  the  throughflow  capacity  of  the  transonic  cascade 
mainly  depends  on  the  actual  dimensions  of  the  aerodynamic  throat. 
Either  an  extremely  small  or  an  extremely  large  inlet  angle  would 
make  the  flow  characteristics  of  the  cascade  worse,  causing  the 
distribution  and  thickness  of  the  boundary  layer  along  the  blade 
surface  to  change.  Especially  when  separation  of  the  boundary 
layer  occurs,  the  size  of  the  actual  throat  is  greatly  reduced. 
Hence  for  a  certain  flow  rate  coefficient  u,  there  exists  an 
optimum  range  of  inlet  angle  3^.  For  the  III*  cascade,  y  always 
exceeds  0.99  for  the  range  lAl°  <  3^<  149°. 


Figure  5.  Variation  of 
approaching  Mach  number 
versus  Regime  Mach  number 
for  I*  Cascade. 


Figure  6.  Variation  of  Loss 
Coefficient,  Choking  Mach 
number  and  flow  rate  coeffi¬ 
cient  versus  inlet  angle  for 
III*  Cascade. 
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1—  Calculated;  2 —  Experiment 


4. 


COEFFICIENT  OF  ENERGY  LOSS  AND  EXIT  FLOW  ANGLE 
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Figure  7  shows  the  experimental  variation  curve  of  energy  loss 
coefficient  versus  operating  Mach  number  ^2t  three  sets 

of  cascades  at  the  design  inlet  flow  angle.  Although  the  order  of 
magnitude  of  the  loss  coefficients  of  the  three  different  cascades 
are  different  from  each  other,  the  trends  of  variation  are  similar. 
These  are  typical  loss  characteristics  of  transonic  turbine  cascades. 
The  following  analysis  of  the  relation  between  the  variation  of  the 
loss  coefficient  and  the  development  of  wave  systems  are  performed 
with  the  II*  cascade  as  an  example. 


Figure  7.  Experimental  relation  between  energy 
loss  coefficient  and  operating  Mach  number. 

1 —  Cascade 


In  the  subcrltical  region  0.69),  the  loss 'coefficient 

due  to  friction  and  wake  mixing  is  of  the  order  2%.  After  reaching 
the  supercritical  condition,  local  supersonic  regions  appear  on 
the  top  of  the  blade,  producing  A  waves.  The  addition  of  shock  wave 
loss  gradually  increases  the  total  loss.  V/hen  M2^  =  1,  the  choking 
condition  is  reached  and  the  tail  wave  meets  the  X  waves  from  the 
blade  top.  The  increase  of  shock  wave  loss  causes  the  early  separa¬ 
tion  of  the  boundary  layer  along  the  blade  top,  increasing  the 
energy  loss.  The  order  of  magnitude  of  ^he  loss  coefficient 
reaches  4%,  which  is  about  twice  that  for  the  subsonic  operating 


79 


condition.  When  M2^,  is  further  increased,  a  system  of  expansion 
waves  occurs,  weakening  the  waves  on  the  blade  top.  At  the  same 
time,  the  incident  point-  of  the  right-running  tail  wave  on  the 
blade  top  is  displaced  further  downstream,  causing  the  boundary 
layer  separation  region  to  move  toward  the  exit  and  the  energy 
loss  to  drop.  For  the  range  ^2x,  ~  ^'07  -  1.35,  the  measured  loss 
coefficient  of  the  cascade  is  2.5  -  3%,  which  is  slightly  higher 
than  that  of  the  subsonic  operation  condition.  This  shows  that  the 
efficiency  is  not  low  at  all.  For  ^2^  ^  above,  the  inci¬ 

dent  point  of  the  tail  wave  on  the  blade  top  still  displaces  down¬ 
stream,  but  the  loss  increases  slightly  due  to  the  increase  of  the 
shock  strength.  At  about  =  1.^5,  the  tail  wave  propagates 

beyond  the  channel  and  hits  the  subsonic  region  of  the  wake  of  the 
adjacent  blade.  The  signal  of  sudden  pressure  increase  behind  the 
wave  .would  propagate  upstream  through  the  subsonic  region  of  the 
wake,  causing  a  pressure  increase  along  with  a  velocity  drop  at  the 
exit  of  the  cascade  and  a  thickness  increase  and  separation  of  the 
boundary  layer.  At  the  same  time,  the  left-running  tall  wave 
occurring  there  also  displaces  forward  along  the  separation  region 
of  the  boundary  layer,  causing  the  separation  region  to  enlarge  and 
mix  with  the  wake  and  hence  producing  a  great  loss.  From  Figure  7, 
it  can  be  seen  that  the  loss  coefficient  at  such  operating  conditions 
can  Jump  to  over  8^.  As  is  further  Increased,  the  downstream 

displacement  of  the  tail  wave  may  reach  the  supersonic  region  of 
the  wake.  The  signal  of  pressure  increase  behind  the  wave  can  no 
longer  be  propagated  upstream.  The  degree  of  separation  of  the 
boundary  layer  near  the  trailing  edge  is  lowered  gradually  and  the 
flow  condition  is  Improved,  causing  the  loss  coefficient  to  drop 
gradually.  The  operating  Mach  number  at  which  the  loss  coefficient 

starts  to  increase  sharply  is  called  the  limit  operating  Mach  number 

* 

M'2t*  Prom  experiments:  for  I*  cascade  M'.,^  =1.1  ,  for  II  cascade 
M'  =  1.45,  f'or  m  cascade  M'2^=  1.55.  the  above  analysis ,  It  can 

be  concluded" that.  If  the  operating  condition  of  the  cascade  Is  kept 
under  the  limit  operation  Mach  number,  then  the:ioss  coefficient  Is 
by  no  means  significantly  large  and  relatively  high  efficiency  Is  attained 
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Factors  influencing  the  loss  coefficient  of  transonic  turbine 
cascades  are  very  complicated.  The  variation  of  the  strength  and 
the  position  of  the  wave  systems,  their  reflections  and  mutual 
interferences,  especially  the  interference  between  the  shock  wave 
and  the  boundary  layer  along  the  blade,  and  the  mixing  process  of 
the  wake  behind  the  cascade,  etc.,  all  have  a  great  effect  on  the 
characteristics  of  the  energy  loss.  These  factors  make  the  calcu¬ 
lation  of  loss  extremely  complicated  and  difficult.  Hence  nowadays 
people  still  rely  on  experiments  to  determine  the  loss  coefficient 
and  the  limit  operating  Mach  number  ^'2^  transonic  cascades  and 
to  define  the  optimum  operating  range. 

Figure  8  shows  the  variation  of  the  mean  exit  flow  angle  B2Qp 
with  the  operating  Mach  number  for  the  three  sets  of  cascades  at 
design  inlet  flow  angle.  Their  trend  of  variation  is  basically  the 
same.  In  the  subsonic  range,  the  exit  flow  angle  is  basically  con¬ 
stant.  In  the  supercritical  range,  shock  waves  are  generated  in 
the  channel.  Refraction  of  the  streamlines  occurs  when  passing 
through  shock  waves,  producing  an  obvious  change  of  exit  flow  angle. 
Upon  taking  II*  cascade  as  an  example,  it  can  be  observed  that 
when  M2^  >  1,  the  exit  flow  angle  drops  slightly  due  to  the  effect 
of  the  right-running  tall  wave.  But  due  to  the  Interaction  of  the 
shock  wave  generated  from  the  blade  top,  the  reflected  tall  wave 
and  the  expansion  wave  fan,  the  exit  flow  angle  increases  with  the 
increase  of  M2^.  At  about  M2^  =  1.36,  the  exit  flow  angle  reaches 
the  first  peak  value  of  about  35°.  As  M2^  is  further  increased,  the 
incident  point  of  the  right-running  wave  on  the  blade  top  approaches 
the  trailing  edge,  causing  the  boundary  layer  separation  region 
to  displace  downstream.  Hence,  the  flow  remains  closer  to  the 
blade  once  again  and  the  exit  flow  angle  drops  slightly.  Beyond  the 
limit  operating  Mach  number,  the  exit  flow  angle  Increases  sharply 
due  to  the  large  separation  region  at  the  trailing  edge.  When  M^^.  = 
1.53,  the  angle  P'2CP  high  as  4^°,  reaching  a  second  peak 

value . 
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Figure  8.  Tested  curve  of  the  relation  between 
mean  exit  flow  angle  S^pp  and  operating  Mach  number 

"2f 

1 —  Cascade 


5.  OPERATING  CAPABILITY  AND  LIMIT  LOAD 

Operating  -capability  is  an  important  gasdynamic  characteristic 
of  the  cascade.  It  is  expressed  by  the  loading  coefficient 


C,  “  M*  cot  fit  +  A/?  cos/?, 

For  transonic  cascades  operating  a+supercritical  pressure  ratio,  the 
choking  condition  is  reached  and  the  upstream  Mach  number  is 
kept  constant.  Hence  the  inlet  flow  angle  is  constant, 

the  operating  capability  of  the  cascade  can  only  be  expressed  by 
M|  cos  which  is  the  tangential  exit  Mach  number. 

When  the  cascade  is  operating  in  the  subsonic  range,  the 
expansion  process  of  the  flow  is  completed  mainly  in  the  mainstream 
passage  upstream  of  the  throat  of  the  cascade  channel.  After  the 
flow  passes  the  throat,  a  portion  of  the  expansion  is  eliminated 
due  to  a  pressure  increase  and  velocity  decrease  along  the  blade  top. 
With  the  increase  of  the  operating  Mach  number  and  when  the  super¬ 
critical  pressure  ratio  is  reached,  the  cascade  is  operating  at 
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transonic  speed  and  the  flow  at  the  throat  is  sonic.  The  super¬ 
sonic  flow  in  the  divergent  section  continues  to  accelerate  and  ex¬ 
pand,  thus  improving  the  operating  capability  of  the  cascade.  In 
Figure  9,  the  variation  curves  of  with  from  experiments  are 
shown  for  the  three  sets  of  cascades.  It  can  be  observed  that,  when 
>  1,  continues  to  increase,  verifying  that  transonic  cascades 
have  higher  operating  capability.  However,  it  does  not  increase 
indefinitely  with  .  For  each  set  of  cascades,  there  exists  a 
maximum  value  for  C^,  which  is  termed  "the  limit  load".  The  operating 
Mach  number  corresponding  to  the , limit  load  is  called  the  limit 
load  Mach  number.  It  is  clear  from  the  analysis  of  the  wave  systems 
that  the  supersonic  expansion  process  in  the  divergent  section 
relies  mainly  upon  the  expansion  fan  CDE  (Figure  3)  and  their  re¬ 
flected  expansion  waves.  With  the  increase  of  Mach  number,  the 
strength  and  the  Mach  angle  of  the  expansion  wave  fan  also  increase. 

Their  Incident  points  on  the  blade  top  are  thus  displaced  downstream, 
enlarging  the  expansion  ratio  and  hence  the  operating  capability 
of  the  cascade.  But  when  the  last  expansion  wave  CE  hits  the  trail¬ 
ing  edge  of  the  blade,  the  expansion  in  the  divergent  section  ends. 

Prom  Figure  2,  it  is  clear  that  velocity  increases  along  the  blade 
top  until  the  trailing  edge  at  the  exit  is  reached,  making  full  use 
of  the  expansion  function  of  the  divergent  section.  As  the  operating 
Mach  number  is  further  increased,  the  flow  behind  the  cascade  expands 
with  constant  tangential  velocity  due  to  the  enlarged  exit  flow 
angle.  At  this  point,  the  operating  capability  not  only  ceases  to 
Increase,  it  may  in  fact  decrease  slightly  due  to  the  separation 
of  the  boundary  layer  at  the  trailing  edge.  Hence  the  problem  of 
maximum  operating  capability  of  transonic  cascade;  that  is, the 
limit  load  problem,  does  exist. 

The  analysis  and  study  of  the  physical  process  of  the  formation  /  62 
of  limit  loads  in  transonic  cascades  and  the  determination  of  the 
magnitude  of  the  limit  load  and  the  corresponding  limit  load  Mach 
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number  is  very  important .  A  comparison  between  Figure  7  and 
Figure  9  shows  that  the  limit  load  Mach  number  of  the  cascade  and  the 
limit  operating  Mach  number  at  which  loss  starts  to  increase 
sharply  are  very  close  to  each  other.  This  may  be  due  to  the  fact 
that  for  common  cascades,  especially  those  with  blades  of  a  very 
thin  exit  edge,  the  last  expansion  wave  is  very  close  to  the  right- 
running  tall  wave,  and  they  hit  the  trailing  edge  of  the  blade 
at  almost  the  same  operating  conditions.  This  result  enables  us 
to  make  good  use  of  the  advantages  of  transonic  cascades,  namely, 
large  operating  capability  while  performing  with  acceptable 
efficiency.  Currently  approximate  methods  for  the  calculation  of 
limit  load  of  cascades  already  exist.  But  wind  tunnel  tests  are 
still  very  practical  in  the  determination  of  limit  loads  of  tran¬ 
sonic  cascades. 


Figure  9.  Tested  curve  of  relation  between  cascade 
load  coefficient  and  operating  Mach  number  M2^ . 

1 —  Cascade 

6.  LIMIT  INLET  PLOW  ANGLE 

Figure  10  shows  the  =  f(M2|.)  curves  of  the  III*  cascade 

at  different  inlet  flow  angles  obtained  experimentally.  It  is 
observed  that  the  effect  of  the  inlet  flow  angle  becomes  more 
explicit  only  after  >  1.4.  For  >  149°,  the  minimum  value 
of  loss  rises  slightly,  but  the  corresponding  optimum  operating 
Mach  number  Increases  significantly.  At  the  same  time,  the 
limit  operating  Mach  number  also  increases  with  the  inlet  flow 
angle.  This  implies  that  the  practical  operating  range  of  the 
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cascade  is  extended.  This  may  be  due  to  the  extension  of  the 
supersonic  expansion  region  when  the  sonic  line  moves  upstream 
from  the  throat  at  higher  inlet  flow  angle.  But  there  exists  a 
certain  limitation  to  the  inlet  flow  angle  Increase  of  cascades 
with  3^  >90°  and  the  inlet  flow  angle  decrease  of  cascades  with 
<  90°.  Analysis  of  the  III*  cascade  >90°)  shows  that  when 
Bj  increases  to  a  certain  value,  such  that  the  width  of  the 
approaching  stream  is  less  than  the  throat  width  of  the  channel, 
then  the  sonic  line  at  which  subsonic  flow  transits  into  supersonic 
flow  would  displace  forward  up  to  the  guiding  stream  in  front  of 
the  cascade  channel.  ThuS;the  cascade  inlet  flow  is  already 
supersonic.  Two  oblique  shocks  are  generated  from  the  leading  edge 
of  the  blade.  As  illustrated  in  Figure  11,  the  right-running  shock 
wave  destroys  the  uniform  flow  field  at  the  inlet  region  of  the 
cascade,  while  the  left-running  shock  wave  propagates  across  the 
channel  and  produces  a  reflection  wave  after  hitting  the  concave 
surface  of  the  blade.  The  cascade  channel  becomes  entirely  super¬ 
sonic  and  more  complicated  wave  systems  appear.  Their  superposition 
and  their  interference  with  the  boundary  layer  along  the  blade 
cause  the  energy  loss  of  the  cascade  to  increase  sharply.  Figure  12 
shows  the  variation  of  M  along  the  cascade  pitch  of  the  III*  cascade 
operating  at  ^2x,  ~  various  inlet  flow  angles.  When  3^  =  l6l°, 

the  distribution  becomes  non-uniform,  and  when  3^^  *  l62°, 
exceeds  1  at  a  certain  section  of  the  cascade  pitch  and  supersonic 
region  appears.  At  3^  =  l62.6°,  a  sudden  Jump  of  distribution 
occurs,  under  the  interference  of  the  leading  edge  shock  wave.  It 
is  certain  that  energy  loss  is  very  large  at  this  stage.  The  experi¬ 
mental  curve  in  Figure  6  verifies  this  Judgment:  For  3j_  >159°,  the 
loss  starts  to  Increase  sharply,  and  when  3^  *  l6l°  is  reached,  its 
value  has  already  exceeded  10^. 

To  restrict  the  inlet  flow  angle,  the  inlet  flow  angle  at  which 
the  sonic  line  occurs  Just  in  front  of  the  channel,  such  that  the 
approaching  flow  is  supersonic,  is  defined  as  the  limit  inlet  flow 
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Figure  10.  Experimental  variation 
of  loss  coefficient  of  III*  Cascade 
at  var'^ous  inlet  flow  angles. 


Figure  11.  Waves  in  the  Cascade 
when  the  inlet  flow  angle  ex¬ 
ceeds  the  limit. 
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Figure  12.  Distribution  of  approach¬ 
ing  Mach  number  of  III*  Cascade. 

1 —  Cascade  pitch 


angle.  If  the  width  of  the  guiding  stream  at  such  a  condition  is 
assumed  to  be  equal ’to  the  geometric  throat  of  the  channel,  then 
the  limit  inlet  flow  angle  of  the  III  cascade  can  be  calculated 
to  be  “  159°10’.  Actually,  due  to  the  two-dimensional  and 
viscous  effects  of  the  flow,  is  slightly  larger  than  the  com¬ 

puted  value.  Using  the  same  method,  the  limit  inlet  flow  angle  of 
the  I*  cascade  is  also  calculated  to  be  =  38°,  which  is 
close  to  the  measured  value. 
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EXPERIMENTAL  INVESTIGATION  AND  ANALYSIS  OF  THE  AXIAL  FLOW 
LOW  PRESSURE  COMPRESSOR  OF  THE  B,  ENGINE  WITH  INLET  FLOW 
DISTORTIONS 


Gao  Yao-lln,  Cai  Yun-jin,  Wang  Hui-Iln  and  Zhong  Yu-llng 
(Shenyang  Aeroengine  Company) 

ABSTRACT 

This  paper  is  a  summary  of  progress  made  in  exper¬ 
iments  on  inlet  flow  distortions  of  the  low  pressure 
compressor  of  the  B^  engine. 

In  the  experiment,  the  equal  rotational  speed  char¬ 
acteristic  lines  and  the  surge  lines  of  compressors  with 
blades  of  two  different  aspect  ratios  were  recorded  under 
different  inlet  flow  distortions.  Furthermore,  the  cri¬ 
tical  sector  angle  of  distortion  resistance  and  attenuation 
rate  of  flow  distortion  for  each  case  were  determined. 

In  the  paper,  the  experimental  results  were  analyzed.  The 
experimental  results  Indicate  that  the  response  charact¬ 
eristics  of  inlet  flow  distortions  for  compressors  with 
blades  of  different  aspect  ratios  were  not  the  same.  Com¬ 
pared  to  a  large  aspect  ratio,  compressor  blades  of 
smaller  aspect  ratios  have  better  distortion  resistant  capa¬ 
bility  and  larger  surge  margin.  This  improved  the  opera¬ 
ting  characteristics  of  the  engine  which  was  more  suitable 
for  flight  maneuverability  of  airplanes. 

I.  EXPERIMENTAL  OBJECT  AND  DISTORTION  APPARATUS 

The  experiments  were  carried  out  with  a  modified  experimental 
J211-950  unit.  The  modification  consists  of  removing  the  fourth 
stage  and  to  restoring  the  flow  regulating  blades  in  the  third  stage 
to  the  original  position.  Then,  the  corresponding  measurement  sec¬ 
tion  and  the  disturbance  device  are  Installed. 

*Thls  paper  was  read  at  the  Second  Engineering  Thermal  Physics 
Conference  in  Hangzhou  in  November  ]978. 
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The  engine  low  pressure  com¬ 
pressor  is  a  three  stage  transonic  axial 
flow  compressor.  The  only  difference 
in  the  two  experimental  compressors 
mentioned  in  this  paper  is  that  the 
aspect  ratios  ^  of  the  first  stage  work¬ 
ing  blades  were  not  the  same.  One  of 
h  ^ 

them  has  r-  -2.0  and  the  other  one  has 
b  ^ 

^  2.6.  For  convenience,  they  are 

called  Compressor  I  and  Compressor  II 
for  short  from  here  on. 

The  steady  state  inlet  radial  press¬ 
ure  flow  distortion  of  the  compressors  used  in  this  experiment  was 
realized  by  relying  on  a  disturbance  device  upstream.  The  simple 
structures  of  this  device  were  most  commonly  holes  in  a  plate  and  dis¬ 
turbance  net  types.  In  this  experiment,  the  disturbance  net  structure 
was  adopted.  Figure  1  is  a  schematic  diagram  of  the  experimental  unit 
and  the  disturbance  device.  The  disturbance  net  was  installed  at  0.75 
Dj^  in  the  upstream.  The  support  frame  was  composed  of  12  radial  bars. 

A  support  net  woven  of  1  mm  metal  wires  with  20x20  eyelets  was  placed 
on  the  support  frame.  Finally,  the  required  disturbance  net  was  secured 
at  the  support  net  and  the  frame.  Thus,  the  distortion  apparatus  of 
this  experiment  was  created. 

The  quasi-uniform  flow  in  the  paper  referred  to  the  flow  field  of 
the  inlet  flow  after  passing  through  the  support  frame  and  support 
net.  Various  flow  distortions  were  accomplished  by  the  sector  shaped 
disturbance  net  of  different  blockage  ratios.  Simultaneously,  in 
order  to  find  the  critical  disturbance  angle  (sector  angle)  of  the 
compressor,  various  sector  shaped  disturbance  nets  were  prepared  before 
hand.  The  passage  way  blockage  ratios  (%)  for  various  kinds  of  sec¬ 
tor  shaped  disturbance  nets  were  20,  25,  30,  35,  ^0,  64  and  69. 


support  net;  3 — support 
frame 
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In  order  to  satisfy  the  required  flow  distortions,  in  the  exper-  /65 
iment  we  frequently  overlapped  various  sector  disturbance  nets  on 
the  support  frame  and  net.  Actual  results  proved  that  the  effect 
was  good.  Using  only  a  few  disturbance  nets,  we  could  obtain  a  lot 
of  distortions.  In  addition,  in  this  experiment,  we  also  used  struc¬ 
tural  measures  for  adjusting  the  axial  distance  from  the  disturbance 
net  to  the  compressor  inlet  in  order  to  obtain  more  distortions. 

II.  PARAMETRIC  MEASUREMENT  AND  DATA  PROCESSING 

The  parametric  measurement  cross-section  I-I  of  the  compressor 
inlet  was  located  30  mm  in  front  of  the  front  fringe  of  the  root  of 
the  first  stage  blades.  The  disturbance  net  was  installed  at  the 
0-0  cross-section  which  was  ^80  mm  from  the  I-I  cross-section  (see 
Figure  1). 

The  inlet  flow  distortion  experiments  had  placed  greater  require¬ 
ments  on  the  testing  methods.  Not  only  the  flow  parameters  must  be 
accurately  measured  at  the  measuring  point,  but  also  the  distributions 
of  the  measured  flow  parameters  across  the  entire  measuring  cross- 
section  must  be  determined  because  they  could  directly  influence  the 
accuracy  of  the  experimental  results. 

In  the  experiment,  several  fixed  pressure  transducer  probes  were 
installed  on  the  I-I  cross-section  to  facilitate  the  determination 
of  the  total  pressure  distribution  of  the  compressor  inlet  cross- 
section.  The  pressure  probes  were  of  two  types:  one  was  the  ten 
point  synthesized  probes  which  were  placed  radially  along  the  total 
static  pressure  measuring  point  l8  mm  apart;  the  other  was  the  6p 
pressure  probes  distributed  radially  along  the  measuring  points  20  mm 
apart. 

From  Figure  7,  we  can  see  that  in  order  to  accurately  calculate 
the  average  total  pressure  value  P]^  at  the  I-I  cross-section  of 
the  inlet,  it  was  necessary  to  go  through  the  lengthy  circumferential 
and  radial  averaging  process.  For  convenience,  we  selected  the  aver¬ 
age  value  of  the  pressures  measured  by  a  synthesized  probe  and  a  6p 
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probe  behind  the  net  as  the  low  pressure  value  of  the  inlet  cross- 
section  disturbance  region  and  chose  the  average  value  of  total 
pressure  measured  by  the  synthesized  probed  inserted  in  areas  not 
affected  by  the  disturbance  as  the  high  pressure  value  of  the  inlet 
I-I  cross-section.  Then,  the  average  value  of  the  total  pressure  of 
the  flow  at  the  inlet  could  be  obtained  according  to  the  following 
formula : 

«.*  =  ei-si-i?!  +  />•»..(  360®  ~  0°) 

- - - 

360® 

where  8°  — disturbance  sector  angle 

^Imin — flow  average  total  pressure  in  the  disturbance 
neb  region 

* 

?!  — inlet  flow  average  total  pressure  in  the  area  without 

disturbance 


In  the  III-IV  cross-section  located  far  away  from  the  outlet  of 
the  compressor,  a  total  pressure  target  with  nine  measuring  points 
along  the  blade  spacing  was  installed.  It  could  move  along  the  direct¬ 
ion  of  blade  height.  The  total  flow  precsure  was  measured  at  5  radial 
positions  and  their  average  value  was  taken  as  the  total  pressure  of 
the  compressor  outlet.  The  measuring  plane  of  the  compressor  outlet 
flow  distortions  was  II-II  which  was  close  to  the  outlet  of  the  third 
stage  flow  regulating  blades.  Several  pressure  and  temperature  probes 
were  installed  in  that  plane. 


In  order  to  evaluate  the  magnitude  of  the  inlet  flow  distortion, 
the  following  two  distortion  indicators  were  used 


DC* 


D  -  -  f:,.) 

Pu, 


(1) 


(2) 


where  "T  “  (p*  ~  ~  ri  Kf*  ”  ?<)••«  ~  ~  (^r  fOw 

*  360 

and  (/-r-p.)...  were  related  to  the  flow  velocities  of  the  disturbance 
free  flow  and  disturbance  flow,  respectively. 
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The  two  indicators  above  are  both  distortion  indicators  based 
on  the  calculation  of  the  total  pressure  loss.  Their  magnitudes 
demonstrated  the  degrees  of  flow  distortions.  In  this  experiment, 
the  experimenvax  ranges  of  DCq  and  D  were  24-5555  and  8-l8!5  respect¬ 
ively. 

IIT.  EXPERIMENTAL  RESULTS  AND  DISCUSSION 


We  recorded  the  characteristic  lines  and  surge  lines  of  the 
quasi-uniform  flow  of  Compressor  I  and  Compressor  II  under  rotational 
speeds  at  n  =  1.0,  O.98,  0.95,  0.90,  0.85  and  0.80  as  the  base  line 
in  comparison  with  the  performance  of  the  compressor  under -flow  dis¬ 
tortion.  In  addition,  in  order  to  evaluate  the  distortion-resistant 
capability  of  Compressor  I  and  Compressor  II,  we  used  a  fixed  working 
line,  obtained  by  calculation,  to  estimate  the  variations  of  the 
surge  margin  of  the  above  two  compressors  under  various  inlet  flow 
distortions.  It  should  be  pointed  out  that  up  until  the  present  mom¬ 
ent  there  has  not  been  a  standard  for  the  definition  of  surge  margin. 
There  are  various  definitions.  Some  of  them  used  the  working  point 
on  the  equal  rotating  speed  line  and  the  converted  flow  at  the  surge 
point  to  carry  out  an  overall  comparison  with  the  pressure  ratio 
paraimeter,  some  only  used  the  pressure  ratio  parameter  for  comparison; 
some  used  the  pressure  ratio  parameter  at  the  working  point  under  con¬ 
stant  flow  and  at  the  surge  point  for  comparison,  etc.  We,  based  on 
the  actual  situation,  adopted  the  following  definition  of  surge  margin 
to  calculate  the  surge  margin  of  the  compressor. 

Surge  margin  (%)  “>  —  i]  x  100% 

where  and  xS  were  the  converted  flow  and  pressure  ratio  of 

the  compressor,  respectively. 


(I)  Critical  sector  angle,  Q|-]-j^tical 


The  critical  sector  angle  is  one  of  the  Indicators  of  the  dis¬ 
tortion-resistant  capability  of  the  compressor.  If  ®cj.itlcal 
large,  this  Indicates  that  the  response  of  the  compressor  to  the  inlet 


flow  distortion  is  not  sensitive.  It  can  stand  a  larger  area  of 
disturbance  of  the  inlet  flow.  In  this  experiment,  two  pieces  of 
sector  disturbance  nets  (overlapping  each  other)  with  passage  blockage 
ratios  at  25%  and  30%,  respectively,  were  used  to  determine  the  cri¬ 
tical  sector  angle  of  the  compressor.  The  sector  angles  of  these 
nets  were  45°,  60°,  75°,  90°  and  120°.  Figures  2  and  3  are  curves 
obtained  experimentally. 


Prom  Figure  2  we  found  that  the  critical  sector  angle  ^d.^^ical 

of  Compressor  I  was  between  75°  and  90°.  In  the  region  less  than 

9  when  9  was  increased,  the  ratio 

critical’  ’  distorted, surge 


/ 


of  the  compressor  dropped  rapidly  almost  in  a  linear 


quasi, surge 

fashion.  After  reaching  ®Qj,j_tical’  ®  increased  further,  that 

ratio  maintained  a  constant  value.  Figure  3  is  the  corresponding 
experimental  curve  for  Compressor  II  whose  critical  sector  angle  was 
60°.  It  was  smaller  than  the  critical  sector  angle  of  Compressor  I. 


Figure  2.  Experimental  curve 
of  the  critical  disturbance 
angle  of  Compressor  I. 

1 —  jrJ  (distortion, surge)/,* 
(quasi, surge)? 

2 —  disturbance  sector  angle  9 


Figure  3.  Experimental  curve  of 
the  critical  disturbance  angle  of 
Compressor  II. 

1 —  wf  (distortion, surge )/ *« 
(quasi, surge)?; 

2 —  disturbance  sector  angle  9 


TABLE  1.  Critical  sector  angles  of  several  compressors  under 
distortion 


rw  14 

/ 

Afilrtll  ^ 

JJE'itt  I 

60* 

Ml' 

Mi'* 

75*-90** 

60* 

1 — spade;  2 — high  pressure  A  engine;  3 — Compressor  I;  4 — Compressor  II 
^  ^critical 
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Table  1  presents  the  distortion  critical  sector  angles  of  sev¬ 
eral  types  of  compressors.  The  first  three,  due  to  the  difference  in 
experimental  conditions,  will  not  be  compared.  The  latter  two  had 
exactly  the  same  experimental  conditions.  The  only  difference  was 
that  the  aspect  ratios  of  the  first  stage  rotor  blades  were  not  the 
same.  From  the  comparison  of  ^Qp^tical’  showed  that  Compressor  I 
was  not  sensitive  to  the  inlet  flow  distortion  which  was  in  agreement 
with  the  conclusion  in  [4].  This  was  also  consistent  with  the  fact 
that  Compressor  I  had  a  larger  surge  margin  under  flow  distortion 
which  will  be  discussed  later. 


(II)  The  effect  of  inlet  flow  distortion  on  compressor 
characteristics  and  surge  lines 

From  the  experimental  results,  we  found  that  at  higher  rotation¬ 
al  speeds  (n  ^  0.90)  the  characteristic  lines  were  most  steep  while 
under  quasl-homogeneous  flow  and  most  flat  while  under  flow  distortion. 
Along  with  Increasing  inlet  flow  distortion  D,  this  trend  became  more 
obvious  (refer  to  Figures  4-6).  Therefore,  based  on  the  experimental 
results,  there  is  such  a  regularity  in  existence.  This  Indicated 
that  the  distortion  of  inlet  flow  would  cause  the  compressor  perfor¬ 
mance  to  decrease  and  that  the  characteristic  lines  were  greatly  varied. 
For  rotational  speeds  n  =  1.0  and  0.98,  the  flow  region  contained 
was  Increased.  The  medium  rotational  speed  lines  did  not  vary  signi¬ 
ficantly  in  shape,  however,  their  positions  were  shifting  toward  the 
direction  of  lower  performance  but  it  was  not  as  obvious,  as  in  the 
high  speed  case.  This  is  most  probably  because  the  corresponding  dis¬ 
tortions  under  these  rotational  speed  lines  were  small.  It  should  be 
pointed  out  that  in  Figures  4-6  presented  in  this  paper,  for  the 
various  characteristic  lines  in  the  same  figure,  their  corresponding 
inlet  flow  distortion  values  were  not  identical.  With  decreasing  rota¬ 
tional  speed,  the  corresponding  inlet  flow  distortion  was  also  de-  /68 
creased.  Our  experimental  method  could  not  guarantee  to  keep  the 
inlet  flow  distortion  under  various  rotational  speeds  to  be  constant. 

In  fact,  it  was  not  necessary  because  that  would  further  deviate  away 
from  the  actual  working  conditions.'  , 
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<r^  AiijtiSfes:  « =•  w  CsS%(«-i.o}  AjEnttii*JJjJ3i*tt«A 

Figure  4 

1 — (kg/sec);  2 — Compressor_I,  0  quasl-homogeneous  flow, 

A  distorted  flow:  9=90°  D=8/5  (n=1.0)  k  surge  point  of  quasi- 
boraogeneous  flow  of  Compressor  II. 


1 —  Gtiir  (kg/sec);  2 — Compressor  I,  0  quasi-homogeneous  flow, 

A  distorted  flow:  0=90°  D=125J  (n=1.0) 


Figure  4  showed  that  in  a  quasl-homogeneous  flow,  the  shapes  and 
positions  of  the  surge  lines  for  Compressors  I  and  II  were  basically 
the  same.  This  indicated  that  the  effect  of  increasing  blade  chord 
length  to  reduce  the  aspect  ratio  on  the  surge  lines  of  the  compressor 
under  homogeneous  inlet  flow  conditions  is  not  significant. 

In  some  papers  fl-3],  it  was  reported  that  the 

conclusion  "the  stable  working  range  of  transonic  compressors  could 
be  expanded  by  reducing  the  aspect  ratio  of  the  blade"  This  was  veri¬ 
fied  based  on  experimental  results  of  single  stage  transonic  compress¬ 
ors.  However,  from  the  experimental  results  of  the  two  different 
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^  snoii  oUkiSi^JiMiS  fl»90*  iJ  -  I8%(s  -  i.O) 

Figure  6 

1 —  c,„  (kg/sec);  2 — Compressor  I,  0  quasl-homogeneous  flow, 
A  distorted  flow,  0=90°,  D=l8%  (n=1.0) 


TABLE  2.  The  surge  margins  of  Compressors  I  and  II  when  D-=0.98= 
15-16!?.  " 


1.0 

0.98 

0.9U 

0.85 

0.8U 

H.2 

7.2 

.  7.9 

9 

10 

6.1 

9 

10 

7,2 

£  I  «  -  0.98  ti  5-16% 

4 

3 

3 

4.7 

-  o.‘#8  ft  5-15% 

0.5 

-O.-l 

1 

1.8 

1 — condition;  2 — surge  margin  !?;  3 — quasi-homogeneous  flow  (Com- 
pressor_I);  h — £uasi-homogeneous  flow  (Compressor  II)i.  5 — Compressor 
I  when  n=0.98,  D=l6!?;  6 — Compressor  II,  when  n=0.98,  D=15/? 


multi-stage  compressors  with  different  aspect  ratios  in  the  first 
stage  blades,  this  conclusion  was  not  completely  verified. 


Nevertheless,  this  experiment  also  verified  another  important 
phenomenon:  Compressor  I  with  a  smaller  aspect  ratio  blades  compared 

to  Compressor  II  with  larger  aspect  ratio  blades  was  able  to  improve 
the  surge  lines  under  inlet  flow  distortions  more.  The  greater  the 

ft 

distortion  was,  the  more  obvious  the  effect  became  (see  Figures  4-6)  . 
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Unfortunately,  we  did  not  further  use  even  small  aspect  ratio 
blades  in  our  experiments. 

TABLE  3-  The  surge  margins  of  Compressor  I  under  various  distortions 


* 

^  k; — 

1.0 

0.V8 

0.V5 

0.9U 

0.8S 

1 

0.10 

3 

8.3 

7.2 

IS 

WSM 

10 

^  * 

=«  l.O  Sj  K% 

3.5 

3.5 

mm 

9.5 

S  • 

=  l.o  \2% 

2.5 

3.1 

iB 

5.3 

7 

* 

«  1 .0  It 

-1.0 

_  1 

0 

— i 

2 

3.S 

6 

7  n 

•  0.98  5  -  16% 

<  i 

3 

3 

4.7 

1 — condition;  2 — surge_margin_?;  k — when  n=l_.0,  D=85S_;_  5 — when 
n=1.0,  D=12%i  6 — when  n=1.0,  D=l8/?;  7 — when  n=0.98,  D=l6%; 

3 — quasl-homogeneous  flow 

In  fact,  it  has  already  been  verified  in  [1]  that  another  major 
reason  to  use  small  aspect  ratio  blades  in  the  first  stage  of  the 
transonic  compressor  was  to  eliminate  the  danger  of  vibrational  loss 
of  speed  of  the  blades. 

Comparing  Figures  4-6,  we  found  that  for  Compressor  I  or  Com¬ 
pressor  II,  with  increasing  inlet  flow  distortion  D,  the  surge  lines 
were  shifted  towards  the  lower  right  and  the  surge  point  flow  was 
decreased.  The  latter  was  even  more  severe  than  the  former.  From 

Table  2  we  found  that  under  the  same  distortion  D-  .  *  15-16/5, 

ufo  .  9 o 

Compressor  I  still  had  3-4.7?  of  surge  margin  while  Compressor  II 
almost  ran  out  of  any  margin.  From  Table  3  we  found  that  the  surge 
margin  of  a  compressor  decreases  with  increasing  D. 

From  the  experimental  curves  of  *5  ^  obtained  at 

1 - 

(from  previous  page) 

Due  to  the  page  limitations,  the  experimental  curves  corresponding 
to  Compressor  II  were  omitted. 
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rotational  speeds  n  =  0.98  and  0.90  for  Compressors  I  and  II,  we 

found  that  the  value  of  w*  decreased  linearly  with  increasing 

_  surge  _ 

D.  Compressor  I  at  a  rotational  speed  of  n  =  0.98,  in  its  peak  effi¬ 
ciency,  also  decreased  linearly  with  Increasing  D. 

III.  DISTORTION  ATTENUATION  RATE 

For  a  twin  rotor  engine,  the  outlet  flow  distortion  of  the  low 
pressure  compressor  could  also  cause  an  adverse  effect  on  the  perfor¬ 
mance  of  the  high  pressure  compressor.  For  a  single  axis  engine, 
excessive  compressor  outlet  flow  distortion  would  have  an  adverse 
effect  on  combustion  chamber  performance.  Therefore,  in  this  experi¬ 
ment,  preliminary  measurements  were  made  on  the  outlet  flow  of  Com¬ 
pressor  I.  The  results  were  not  satisfactory  because  the  individual 
probes  on  the  II-II  cross-section  were  facing  the  rear  fringe  of  the 
third  stage  flow  regulating  blades  and  the  results  were  affected  by 
the  tail  traces.  However,  after  partial  correction,  it  was  still 
possible  to  approximately  observe  the  pressure  distribution  of  the 
outlet  flow.  In  the  experiment,  we  determined  the  total  outlet 
pressure  when  n  =  0.98  (see  Figure  7).  The  corresponding  distortion 
D  value  was  3-5%-  At  this  time,  the  D  value  of  the  flow  distortion 
with  respect  to  the  inlet  and  outlet  was  10%.  This  showed  that  the 
compressor  inlet  flow  distortion  was  strongly  attenuated  after  passing 
through  the  three  stage  compressor.  The  distortion  value  at  the  out¬ 
let  was  about  one-third  of  that  of  the  inlet. 


Figure  7(a) 
1— condition  point 


Figure  7(b) 
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In  the  meantime,  the  outlet  temperature  distributions  of  Com¬ 
pressor  I, when  n  =  0.98  with  sector  angle  9,  60®  and  90®, were  mea¬ 
sured  in  the  experiment.  The  measured  results  indicated  that  the 
pressure  distortion  of  the  inlet  flow  of  the  compressor  causes  a  dis¬ 
tortion  of  the  outlet  temperature  field.  Furthermore,  the  outlet 
temperature  field  relative  to  the  inlet  pressure  field  corresponds 
to  rotation  in  the  direction  of  rotor  rotation. 

/TO 

IV.  PRELIMINARY  CONCLUSIONS 

(I)  The  distortion  critical  disturbance  angle  of  Compressor  I  is 

®Icritical  “  75°-90®,  and  encritical  Compressor  II  =  60®.  The 
former  was  not  sensitive  to  inlet  flow  distortion  and  could  stand 
inlet  flow  distortion  disturbances  over  a  larger  area.  It  also  bad 
more  intense  distortion  attenuation  capability. 

(II)  Under  Identical  conditions,  compressors  with  smaller  aspect 
ratio  blades,  relative  to  compressors  with  larger  aspect  ratio  blades, 
have  a  larger  stable  working  range,  especially  under  inlet  flow  dis¬ 
tortion  conditions.  They  had  larger  surge  margins  which  meant  that 
the  distortion-resistant  capability  was  better.  It  was  more  suitable 
for  airplane  maneuvering. 

(III)  The  pressure  distortion  of  the  inlet  flow  of  the  compressor 
causes  a  temperature  distortion  of  the  outlet  flow  of  the  compressor. 
The  outlet  ten?3erabure  relative  to  the  inlet  pressure  corresponds  to 
a  turn  in  the  direction  of  rotor  rotation. 

V.  CONCLUSIONS 

Experimental  work  on  inlet  flow  distortion  of  compressors  has 
Just  begun.  Due  to  lack  of  experience  and  other  objective  reasons, 
there  were  imperfections  in  all  the  links  which  are  yet  to  be  Improved 
upon.  The  readers  are  welcome  to  present  their  valuable  suggestions. 
This  experiment  was  carried  out  on  the  compressor  testing  device  at 
the  Research  Institute  of  Shenyang  Aeroengine  Company.  The  Comrades 
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of  the  second  group  in  the  third  room  have  performed  a  larpe  am.ount 
of  work.  We  wish  to  express  our  thanks  here. 


With  the  increasing  flight  range  of  modern  airplane  flights, 
more  rigorous  requirements  are  placed  on  the  operating  characteristics 
of  engines.  Hence,  Improvement  of  distortion  tolerance  has  already 
become  one  of  the  major  topics  to  be  studied.  In  order  to  further 
improve  the  compressor  distortion  tolerance  and  to  expand  the  maneuv¬ 
ering  range  of  the  flight  of  airplanes,  we  are  carrying  out  experi¬ 
mental  studies  in  various  casing  treatments.  This  report  only  con¬ 
tained  part  of  the  experimental  plan. 
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EXPERIMENTAL  INVESTIGATIONS  ON  THE  TURBULENT  FILM  COOLING 
EFFECTIVENESS  OVER  AN  ADIABATIC  FLAT  PLATE* 


Ko  Shao-yen,  Liu  Deng-ying,  Huang  Liang-fu,  Zhang  Zbeng-fang  and 
Pan  Wen-xin  (Academia  Slnlca) 

ABSTRACT 

Experimental  methods  were  used  to  determine  the 
adiabatic  temperature  ratio  n  of  turbulent  flow  cool¬ 
ing  over  a  flat  plate.  The  mainstream  velocities  In 
the  experiment  were  5^-101  (m/sec).  The  jet  velocities 
were  15-1^3  (m/sec).  The  mass  flow  rates  M  were  0.2- 
3.2.  The  mainstream  temperature  variation  range  was  40- 
148  (°C).  The  Jet  temperature  variation  range  was  24- 
67  (°C).  The  experimental  results  are  expressed  using 
several  correlations. 

INTRODUCTION 

A  cooling  Jet  was  Injected  In  the  tangential  direction  near  the 
wall  surface  to  separate  the  high  temperature  gas  from  the  wall  sur¬ 
face.  This  type  of  protective  cooling  Is  flow  film  cooling.  Flow 
film  cooling  has  been  widely  used  In  main  combustion  chambers  and/or 
the  acceleration  combustion  chambers  of  modern  gas  burning  turbine 
Jet  engines  and  the  pulse  pressure  engines.  The  combustion  chambers 
and  nozzles  of  rocket  engines  can  also  use  turbulent  film  cooling. 

Recently,  some  authors  carried  out  more  studies  In  the  turbulent 
film  cooling  process.  However,  the  results  obtained  disagreed  widely 
and  could  not  satisfy  the  design  requirements.  For  example.  In  the 
actual  application,  the  effective  cooling  length  x/s  Is  relatively 
short  and  the  flow  rate  M  Is  approximately  around  1.  However,  It  was 


This  paper  was  presented  In  the  Second  Engineering  Thermal  Physics 
Conference  In  Han  zhou  In  ’ovember  1978. 
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reported  in  the  literature  [7-9]  that  the  mass  flow  rate  M  was  low 
and  the  experimental  plate  was  too  long  due  to  the  limitation  of 
experimental  conditions.  The  rear  part  of  the  experimental  plate 
has  already  entered  the  turbulent  boundary  laminar  flow  region  and 
the  flow  film  cooling  effect  was  very  poor.  Under  this  condition, 
the  adiabatic  temperature  ratio  formula  obtained  could  not  agree  with 
the  actual  design  requirements  [5].  The  experimental  parameters  were 
high  and  were  close  to  the  condition  of  actual  engines.  But,  the 
actual  temperature  rise  of  the  jet  flow  in  the  jet  channel  was  not 
considered  when  the  adiabatic  temperature  ratio  was  organized.  There¬ 
fore,  the  results  obtained  were  on  the  low  side  [4,6].  The  more 
practical  adiabatic  temperature  ratio  correlation  formula  were  pre¬ 
sented. 

According  to  the  above  situation  described,  this  paper  studied 
the  calculating  method  of  the  adiabatic  temperature  ratio — emphasiz¬ 
ing  the  organization  of  the  correlation  formula  of  the  adiabatic  tem¬ 
perature  ratio  in  the  jet  mixing  region.  Part  of  the  experimental 
results  were  published  in  [11]. 

BASIC  EQUATIONS  OP  FILM  COOLING 

The  flow  of  flat  plate  turbulent  film  cooling  is  as  shown  in 
Figure  1.  Under  normal  conditions,  y^  is  far  less  than  the  diameter 
of  the  flaming  cylinder.  Therefore,  the  flow  film  cooling  flow  Inside 
the  flame  cylinder  may  be  treated  as  the  flat  plate  flow. 

The  velocity  and  temperture  distributions  in  the  jet  mixing 
region  are: 


(1) 

r 

‘  J 

t  --  Uw  _  1 

(2) 

where  u^  is  the  flow  speed  of  the  boundary  layer  6  of  the  flow. 
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Figure  1.  Simplified  diagram  of  the  flat  plate  turbulent 
film  cooling  flow. 

1 — initial  section;  2 — main  section;  3 — boundary  laminar  flow  region 

The  adiabatic  temperature  ratio  n  is  an  important  dimensionless 
temperature  ratio  correlation  of  an  adiabatic  wall  flow  film  cooling. 
Its  definition  equation  is 

=  —  /..)/(/-  -  /.)  (  3  ) 

where  is  the  adiabatic  wall  temperature  which  is  the  wall  temper¬ 

ature  of  flow  film  cooling  when  the  experimental  plate  is  not  heated; 

».  is  the  average  temperature  of  the  jet  gas  being  injected  out  of 
the  jet  channel  (x  =  0).  In  general,  n  is  a  function  with  the  follow¬ 
ing  correlations: 

1,  -  /(xA,  M .  Rt„  Pr./PrJ)  (  ) 

EXPERIMENTAL  APPARATUS  AND  DETERMINATION  OF  THE  FLOW 

As  for  the  mainstream  flow  entering  the  experimental  section 
(see  Figure  2),  after  being  heated  by  the  electric  heater,  the  flow 
is  homogeneous  and  the  turbulence  has  been  reduced.  The  jet  flow  was 
Injected  into  the  experimental  section  in  the  tangential  direction 
after  pressure  reduction  stabilization.  The  jet  channel  light  s  was 
limited  by  a  jet  spacing  plate  with  a  thickness  A  =  2  and  s  *  2.3  nun. 
The  spacing  plates  used  were  the  removalbe  flat  head  and  sharp  head 
types . 
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The  experimental  flow  plate  was 
398  mm  long,  12  mm  wide  and  the  mat¬ 
erial  was  2xl8  stainless  steel.  In 
order  to  reduce  the  transverse  heat 
flow  along  the  plate  length,  there 
were  ^0  deep  grooves  on  the  back  of 
the  plate.  The  bottom  of  the  groove 
was  0.4-0. 5  mm  away  from  the  surface. 

The  grooves  began  from  the  front  and 
were  distributed  first  densely  and 
then  far  apart.  The  (|)  0.5  mm  nickel 
chromium-nickel  aluminum  thermo¬ 
couples  used  to  measure  the  wall  temp¬ 
erature  were  installed  on  the  bottom  of  22  grooves  on  the  back  of  the 
plate.  Sufficient  heat  Insulating  material  was  filled  under  the  exper¬ 
imental  flat  plate  to  create  a  good  thermal  Insulating  condition  for 
the  experimental  flat  plate. 

This  experiment  used  the  determined  velocity  distribution  and 
temperature  distribution  to  evaluate  the  flow  field  in  the  channel  of 
the  experimental  section.  The  temperature  distribution  was  shown  in 
Figure  3  which  agreed  with  the  temperature  distribution  equation  (2) 
approximately  [3]. 

PHOTOGRAPHY  OF  THE  TRACING  DEVICE 

An  80  mm  diameter  light  beam  tracing  device  was  used  to  carry 
out  observations  of  the  flow  field  near  the  Jet  tank  and  a  series  of 
photographs  had  been  taken  for  the  jet  flow  mixed  boundaries.  Figure 
4  shows  two  typical  photographs . 

Figure  4a  showed  that  when  the  sharp  headed  spacing  plate  was  /75 
used,  there  was  no  displacement  in  the  Jet  extremum  point.  The  flow 
was  more  homogeneous  and  the  inner  and  outer  boundaries  of  the  mixing 
region  were  both  straight  lines.  However,  the  disadvantage  was  that 
the  effective  Jet  groove  height  was  increased  near  the  tip  of  the 
plate.  Correspondingly,  the  Jet  flow  velocity  was  decreased  and  the 
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Figure  2.  The  experimental 
section  and  the  tracing 
device 


Figure  3-  The  dimensionless  temperature  distribution  of 
adiabatic  flat  plate  under  flow  film  cooling. 

1 — symbol;  2 — velocity  (m  /sec);  3 — temperature  (°C);  4 — distance 
(mm);  5 — slot  height  (mm);  6 — remarks;  7 — experimental  sequence 


1/.-  ».0m/i  V,  -  l25.Uiu/> 
A.--I7.2V  /, -25.0X: 

V  -  i.tO 


b 

I/.  »  7l.lm/i  (/,  »  62.0iii/> 

.If  >  o.ys 


Figure  4.  Tracing  photographs  a,  b 
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Figure  5-  Variation  of  the  tangent  of  the  angle  between 
the  mixing  region  boundaries  with  m  measured  by  the  tracing 
method. 

1 — shadow  device,  reference  [l]j  2 — tracing  device,  this  work 


Figure  6.  Variation  of  the 
tangent  of  the  Jet  stream 
mixing  region  boundary  angle 
with  m. 

1 — reference  [5] 


Figure  7.  Variation  of  the 
initial  section  length  in  the 
Jet  flow  mixing  region  with  m. 

1 — thermal  couple;  2--flve;  3 — 
four;  4 — three;  5 — two;  6 — 
experimental  point;  7 — calculated 
point 
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flow  film  cooling  effect  was  reduced.  Figure  is  the  formation 
condition  of  the  Jet  boundary  when  the  flat  top  plate  height  was  M. 
The  mixing  boundary  line  was  approximately  in  the  form  of  a  straight 
line . 


Figure  5  is  the  variation  regularity  of  the  tangent  of  the 
angle  between  the  inner  and  outer  boundaries  (*i  =•  tana,,  b,  =»  tanu,) 
with  m  in  the  Jet  mixing  region  measured  from  a  batch  of  photographs. 
The  results  approximately  were  in  agreement  with  [1].  Figure  6  shows 
the  tangents  of  the  Jet  flow  boundary  angles,  i.e.,  (tan  a),j,  and 
(tan  measured  according  to  the  temperature  and  velocity  distri¬ 

butions.  At  this  time,  the  mixing  boundary  was  selected  at  955?  of 
the  mainstream  temperature  and  velocity.  When  m  >  1,  the  results  of 
the  inverse  of  m  and  thdse  of  m  <  1  were  symmetrical.  The  photographs 
of  this  type  of  flow  film  cooling  and  boundary  formation  were  rarely 
reported  in  the  literature. 

THE  PROBLEM  OF  t^ 

The  average  Jet  temperature  E  of  the  Jet  tank  outlet  is  an 

o 

Important  parameter  in  flow  film  cooling.  The  difference  is  that  the 

different  measurement  methods  of  ^  was  one  of  the  factors  which  caused 

s 

the  inconsistency  of  the  results  in  the  literature. 

In  this  paper,  the  lowest  adiabatic  wall  temperature  measured, 
using  the  second  through  the  fourth  thermocouples  (their  positions 
were  x/s  =  0.217,  5-65,  8.26)  near  the  Jet  tank, was  used  as  t  .  There 

o 

were  two  major  reasons  for  this  selection  method:  (1)  the  insulation 
of  the  experimental  plate  was  really  good.  Experimental  results 
indicated  that  the  temperature  distribution  of  the  flow  near  the  slot 
was  perpendicular  to  the  tested  flat  plate.  Hence,  it  is  possible  to 
use  an  adiabatic  wall  temperature  to  express  the  flow  temperature  on 
the  cross-section  near  the  wall  surface.  (2)  from  the  model  in 
Figure  1  we  found  that  in  the  Initial  section  range  the  unpenetrated 
Jet  flow  covered  the  adiabatic  plate.  According  to  the  half  wall  Jet 
stream  theory,  the  temperature  at  any  point  on  the  adiabatic  plate  at 
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this  time  is  equal  to  the  temperature  of  the  Jet  stream.  As  shown 
in  Figure  7,  the  initial  section  length  of  the  Jet  stream  mixing 
region  lengthens  with  increasing  ra  =  u^/u^. 

ADIABATIC  TEMPERATURE  RATIO 

The  great  amount  of  experimental  results  in  this  work  should 
that  when  the  mainstream  and  Jet  stream  were  air,  the  adiabatic  tem¬ 
perature  ratio  n  mainly  varied  with  x/s  and  M.  Prom  Figure  8j  we 

found  that  the  exponent  of  x/M  varied  with  M.  The  exponent  of 

s 

x/Mg  was  a  constant  because  its  experimental  M  value  was  also  a  con¬ 
stant  [7,10]. 

Figure  9  gives  the  variation  of  n  with  M.  Obviously,  M  <  1 
and  M  >  1  were  two  different  regions.  Furthermore,  we  found  that  in 
the  initial  section  and  the  beginning  of  the  main  section  of  the  Jet 
stream  (about  x/s  <  44.8),  after  M  reached  a  certain  value,  even 
if  when  M  was  increased,  n  no  longer  continued  to  Increase,  but  sta¬ 
bilized  at  a  given  value.  This  indicated  that  at  this  time, even  by 
increasing  Jet  flow, the  cooling  effect  no  longer  would  be  improved. 
This  region  was  the  right  upper  corner  triangular  area  at  M  <  1  in 
Figure  9-  The  reflection  point  of  n  when  reaching  a  stabilized  value 
was  on  a  straight  line. 

From  the  experimental  results,  we  obtained  the  following  adia¬ 
batic  temperature  ratio  correlation  formula:  (1)  when  M  <  1 

II  -  2.73M'-«(rA)-*-*  <  20)  (  5  ) 

,  -  4.55W*-'(x,  A-*-*  (20  <  xls  <  150)  (  6  ) 

,  »  :.04(x//)--'’  (7) 

when  x/s  <  65.2  and  0.45  <  M  <  1,  the  lower  value  of  n  from  Equations 
(16)  and  (7)  was  chosen. 

(2)  when  1  <  M  <  2: 

,1  —  l.96Af"-”(»/>)"^'*  (*A<150)  (8) 

(3)  when  2  <  M  <  3.5: 
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I 

t 

j,  =  2.7I(,/0-"-*  (xA<150)  (9)  ^ 

From  Figure  10,  the  direct  comparison  of  equations  (6)  and  (8) 
in  this  work  and  several  other  equations  in  the  literature  were 
given.  In  the  figure,  we  chose  M  =  0.3,  0.8,  1.3  and  1.5.  x/s  was 
chosen  to  be  between  10-100.  The  formula  used  are  listed  in  the 
following  table. 


Comparison  table  of  adiabatic  temperature  ratio  formula 


^  a 

IB  m 

Ix^xm 

Sees*  (IV54) 

y) »  2.20(*/.W/)-»  ' 

tn 

Harineu  Dirkeback  F.ckerl(l961^ 

ju-o.2a 

x/Mi  wm  60— lOUU 

Pi 

Seban  (1!;60) 

x/Afi*  »>56 

w 

it(») 

M<»,20<jr/,<150 

i«;m<2 

1 — author;  2 — formula;  3 — range;  4 — [reference];  5 — this  work; 
6 — equation 


From  Figure  10,  we  could  see 
that  for  various  M  values  the  for¬ 
mula  given  in  this  paper  gave  more 
reasonable  results.  Especially 
when  M  >  1,  equation  (8)  provided  a 
reasonable  result  while  the  other 
formulas  were  not  applicable.  When 
M  =  0.3,  the  results  by  Hartnett  and 

Seban  are  more  or  less  in  agreanent  with  versus*  M^measured?^ ' 

equation  (6)  in  this  paper.  However,  as  described  earlier  in  the 
paper,  low  M  values  do  not  have  any  practical  value. 

CONCLUSIONS 

1.  In  the  given  experimental  range,  the  adiabatic  temperature 
range  can  be  expressed  by  equations  (5)-(9).  Among  them,  equations 
(6)  and  (8)  are  more  Important.  The  remaining  three  equations  are 


Figure  8.  A  ref.  [8] 
ref.  [10] 
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the  auxiliary  equations  of  the  two  above  equations. 

2.  Prom  Figures  8,  9  and  10,  we  can  find  that  the  x/M  correl- 

s 

at  ion  commonly  used  in  the  literature  and  the  -0.8  exponent  of  this 
correlation  are  mostly  limited  to  M  =  -.01-0.5  and  in  the  range  when 
x/s  is  high.  This  range  does  not  agree  with  the  actual  condition  of 

turbine  Jet  engine  flame  cylinders.  Therefore,  finally  in  this 

paper,  these  types  of  correlation  formulas  were  used. 

TABLE  OF  SYMBOLS 


m  =  u„/u_  flow  velocity  ratio 

S  00 

M  =  u^r„/u  r„  mass  flow  rate 
S  S  00  00 

s  Jet  tank  height,  mm 

X  length  of  the  flat  plate  starting  from  the  Jet 

tank  outlet,  mm 


y,  distance  from  the  outer  boundary  of  the  Jet  stream 

mixing  region  to  the  flat  plate,  mm 

a,.  Op  angle  between  the  inner  and  outer  boundaries  of  the 

Jet  stream  mixing  region 

Subscripts : 


aw 


adiabatic  condition 


s 


Jet  flow  condition 
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EXPERIENCE  ON  DEVELOPING  AN  ANNULAR  COMBUSTOR  OP  "A"  TYPE* 
TURBOJET  ENGINE 


Zhang  Bao-cbeng,  Hou  Yu-tian,  Huang  Dlan-fu  and  Huang  Zhao-xlang 
(Shenyang  Aeroengine  Company)  (Academia  Sinica) 

ABSTRACT 

This  paper  briefly  described  the  technological 
targets  and  the  design  characteristics  of  the  annular 
flame  tube  of  the  type  A  engine.  Emphasis  was  placed 
on  explaining  that  the  adjustment  of  the  flame  tube  on 
an  engine  test  stand  can  control  the  flow  distribution. 

In  the  adjustment,  the  flow  distribution  in  the  main 
combustion  zone  dilution  zone — auxiliary  zone  should  be 
close  to  the  ratio  of  30%-30%-^0% .  Air  used  in  flow 
film  cooling  should  be  controlled  at  under  2555.  In  the 
paper,  it  was  also  shown  that  by  controlling  the  rele¬ 
vant  parameters  of  temperature,  the  homogeneity  of  the 
temperature  field  could  be  improved.  By  choosing  the 
proper  flow  passing  through  the  inner  and  outer  annular 
ducts,  it  was  possible  to  ensure  the  better  distribution 
of  the  combustion  gas  temperature  along  the  blade  height 
direction.  In  the  adjustment  method,  one  adjusts  temp- 
erature  t^^  to  Indirectly  adjust  t^.  The  experiment 
results  realized  the  design  targets  of  raising  R  by  2j5 
and  reducing  Cj^  by  1.5%  • 

INTRODUCTION 

The  annular  combustion  chamber  has  more  advantages  over  the  loop 
combustion  chamber  (Figure  1),  especially  in  the  low  flow  loss  of  the 
combustion  chamber,  homogeneous  outlet  temperature,  better  flame 
connection,  low  breakdown  rate  and  simplicity  in  technology.  If  one 

» - 

This  paper  was  presented  in  the  Second  National  Engineering  Thermal 
Physics  Conference  in  Hang  zbou  in  November  1978. 
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considers  an  annulus  and  a  ring 
of  the  same  Inner  and  outer  dia¬ 
meters,  If  the  same  unit  area 
cooling  air  mass  Is  required, 
the  annulus  can  save  one-third 
of  the  cooling  air  [1].  Thus,  It 
Is  more  advantageous  to  organize 
the  combustion  and  to  use  the 
gas  economically. 


At  present,  the  advanced 
engines  abroad  have  already  wide¬ 
ly  adopted  this  annular  combus¬ 
tion  chamber.  The  change  of  type 
A  turbojet  engine  (from  here  on 
It  Is  called  type  A  engine  for  short)  to  the  annular  flame  tube  was 
based  on  the  Improvement  of  the  flame  tube  structural  design  In  order 
to  raise  the  component  performance,  to  reduce  the  weight,  to  mini¬ 
mize  breakdown  and  to  extend  the  useful  life. 


Figure  1.  Comparison  of  two  types 
of  combustion  chambers. 

1 — annulus;  2 — ring;  3 — external 
shell  of  the  combustion  chamber; 

4 — flame  connecting  tube;  5 — flame 
tube;  6 — vortex  device  and  fuel 
nozzle 


I.  DETERMINATION  OF  TECHNICAL  OBJECTIVES 


The  technological  objectives  of  the  design  adjustment  of  type  A 
annular  flame  tube  (annulus  I  for  short  from  here  on)  are  as  follows: 

1)  The  annulus  I  should  raise  R  by  2!?  and  reduce  C^^  1.5%  for 
type  A  engines. 

2)  The  ground  starting  characteristics  when  using  two  Ignition 
devices  (originally  four  Ignitors  were  used)  must  not  be  lower  than 
those  of  the  original  flame  tube  for  annulus  I.  It  can  guarantee 

the  restarting  of  the  engine  at  a  flight  altitude  of  under  9000  meters. 

3)  Good  lean  and  rock  fuel  stopping  characteristics.  The  eng¬ 
ine  should  not  stop  when  the  throttle  Is  rapidly  cut  and  no  surge  and 
cracking  sound  should  occur  during  acceleration. 
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4)  The  outlet  temperature  t^  dlstrihution  is: 

9,  <  20%;  ill,.//,*,  =»  1.05. 

X 

Based  on  the  experience  on  the  indirect  adjustment  of  t_  by  adjust- 

X  X  J 

ing  t|j  distribution  for  type  A  engine,  it  must  make  t^^  satisfy  the 
following  to  realize  the  above  objectives: 

a. 

b .  > ■lox: .  i •  <  6201:: ;  -  /j:,,  <  1 20<t: ; 

c.  l/jVi  -  /r,i|  <  lO'-C;  AT  =  !/•„  -  /•.!  <  lO'X;; 

4  5)  The  weight  is  lighter  by  3  kilograms  as  compared  to  the 

ring  tube  for  each  unit. 

6)  The  useful  life  is  250-300  hours. 

II.  ADJUSTMENT  METHODS 

The  flime  tube  testing  was  divided  into  a  simulated  single  tube 
test  and  the  complete  engine  test  stand  adjustment.  Single  tube 

X 

tests  can  directly  measure  the  n-,*  <7-3 »  stable  combustion  range  and 
X  •  3  o 

t^  distribution  and  after  the  flame  observation  capability. 

The  complete  engine  test  stand  adjustment  does  not  have  any  sim¬ 
ulated  problem.  It  can  provide  the  actual  flow  field  and  the  flow 
parameters.  It  can  also  directly  access  the  effectiveness  of  the 
improvement  procedures  made  on  the  combustion  chamber.  For  the 
Improvement  of  the  engine  combustion  chamber  already  in  mass  product¬ 
ion,  it  is  more  suitable  for  complete  engine  test  stand  adjustment. 

4  In  general,  it  is  possible  to  determine  the  individual  adjust¬ 

ment  of  the  effectiveness  of  the  flame  tube  of  the  complete  engine 
through  the  following  parametric  evaluation  when  other  components  are 
kept  unchanged: 

1)  The  increase  in  R  can  be  considered  as  the  improvement  in 

X  X 

By  rough  calculation,  to  improve  by  1?  for  type  A  engines, 
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Figure  2.  vs. 

for  type  A  engines 


curves 


Figure  3-  Schematic  diagram  of 
fuel  gas  passing  through  the 
turbine  channel  o 

1 — measuring  position  of  t^  200  point 


R  can  be  improved  by  0.65%. 


2)  The  decrease  in  C-  is  primarily  the  result  of  the  Improvement 

* 

in  The  effect  of  is  small.  By  rough  estimation,  to  raise 

by  1%,  can  approximately  decrease  by  1%. 

.  » 

3)  For  type  A  medium  pressure  ratio  engines,  the  measured  t^^j 

curves  (Figure  2)  when  the  atmospheric  temperature  t,  has  at  -20°- 

* 

30®C  showed  that;  when  n^  remained  unchanged  t^^pj  did  not  vary  sig¬ 
nificantly  with  t^q.  Hence,  it  is  possible  to  use  this  characteris¬ 
tic  to  determine  the  temperature  of  this  type  of  engine.  For  exper¬ 
iments  carried  out  under  different  atmospheric  conditions,  the  abso- 

« 

lute  values  of  t^j  are  comparable. 

«  « 

Between  t^-t^,  there  must  be  some  corresponding  relations  as 
shown  in  Figure  3:  the  high  temperature  fuel  gas  after  passing  through 
the  guide  and  turbine  to  reach  the  measuring  thermocouple  position  for 
must  have  been  shifted  by  an  angle  in  the  direction  of  rotation  of 
the  guide  and  the  turbine.  The  rotation  angle  of  type  A  engine  was 
experimentally  determined  to  be  about  36®.  Hence,  based  on  the 

5,,  and  the  high  temperature  distribution  along  the  blade  height  mea- 

^  » 

sured  on  the  test  stand,  the  t^  distribution  can  be  determined  to  eva¬ 
luate  the  effectiveness  of  Improvement  by  "operating  on"  the  flame 
tube. 
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4)  Through  tho  acceleration,  deceleration,  maneuverability  of 
the  engine  to  test  the  combustion  stability  and  temperature  rise 
rate. 


5)  Through  the  test  stand  starting  test  to  evaluate  the  Ignl-  /97 
tlon  and  flame  connection  characteristics  of  the  flame  tube. 

6)  Through  the  painting  of  temperature  Indicating  paint  on  the 
wall  to  measure  the  temperature  distribution  on  the  flame  tube  wall 
when  the  test  stand  was  under  maximum  operating  conditions  to  expose 
the  local  high  temperature  zone  on  the  wall  and  to  observe  the  effect¬ 
iveness  of  the  Improvement  measures. 

III.  DESIGN  CHARACTERISTICS  AND  ANALYSIS 

1.  The  use  of  elliptical  tips 

Annulus  I  used  10  elliptical  tips  to  satisfy  the  structural  con¬ 
dition  of  the  original  10  full  nozzles  of  type  A  engine.  The  photo¬ 
graph  slows  Its  external  appearance. 

The  use  of  an  elliptical  tip  avoided  the  high  flame  tube  cavity 
caused  by  using  a  round  tip  (Figure  4)  to  make  the  structure  more  com¬ 
pact.  In  order  to  obtain  a  reasonable  structure  and  to  extend  the 
lifetime.  It  Is  necessary  to  keep  the_  heat  capacity  lower  than  that 
of  the  original  flame  tube.  Let  us  take  the  cavity  height  as  127  mm 
(the  original  flame  tube  cavity  height  was  136  iran)  for  annulus  I. 

The  length  to  height  ratio  of  the  flame  tube  was  3.55  which  agreed 
with  the  2. 6-5.0  range  for  an  ordinary  flame  tube  length  to  height 
ratio. 

The  use  of  elliptical  tips  also  could  sufficiently  utilize  the 
radial  space  around  the  tips  to  Increase  the  volume  of  the  main  com¬ 
bustion  zone.  According  to  a  calculation,  the  volume  of  the  ellip¬ 
tical  tips  of  annulus  I  Increased  by  35.5/f'  over  that  of  a  seven  stage 
gas  film  flame  tube  (seven  stage  for  short  from  here  on).  This  is 
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Figure  Comparison  of  elliptical  and  round  tips 
1 — photograph;  2 — mm 


advantageous  for  the  improvement  of 

combustion  efficiency  to  reduce  C_. 

n 

2.  Hole  opening  forms 

Small  holes  at  the  tip  were 
used  for  annulus  I.  The  main  combus¬ 
tion  holes  were  flange  holes.  The 
tube  body  takes  in  air  vertically. 

In  the  dilution  section,-  large  holes 
(inner  ring)  were  used  (see  the 
photograph) . 


Figure  5.  Schematic  diagram 
of  the  flanging  hole  and  the 
funnel 

1 — flanging  hole;  2 — funnel 


Compared  with  flush  holes,  the  use  of  funnel  shaped  main  com¬ 
bustion  holes  would  improve  reduce  C^  by  about  0.005- 

0.01.  The  use  of  flanging  holes  as  compared  to  funnel  shaped  holes 
(Figure  5)  would  reduce  C^  by  about  another  0.01.  The  flanging  has 
a  higher  flow  coefficient  which  makes  it  easier  for  the  Jet  stream  to 
reach  the  center  of  the  flame  tube.  This  enables  the  elimination  of  the 
local  fuel  rich  vortex  nucleus  and  the  Improvement  in  combustion 
efficiency . 
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Figure  6.  Flame  tube  of  A 
type  A  engine. 

Table  1,  Temperature  field  indi¬ 
cators  using  various  apertures 


^  11  & 

%  c 

14.5 

744 

\5.5 

2tl.O 

775 

19.0 

15-^21 

l--aperture;  2 — ratio  of  heat 
area  vs.  cross-sectional  area 
of  thermal  ejection 


Figure  7.  Temperature  distribution 
along  the  blade  in  a  flame  tube 
1 — blade  height  depth;  2 — advanced; 
3 — lagging  behind;  4 — bood  distri¬ 
bution;  5 — nozzle  diameter  under 
optimal  condition 


c«)  •»•> 

Figure  8.  Improvement  of  the  rotat¬ 
ing  connection  section  in  the  tip 
of  the  annulus  I  flame  tube 
1— (a)  before;  2 — (b)  after 


The  adjustment  of  the  dilution  hole  would  directly  affect  the 
outlet  temperature  and  radial  temperature  distributions.  Because  the 
difference  of  the  radial  positions  of  the  center  line  of  the  outlet 
channel  of  the  type  A  engine  compressor  and  that  of  the  flame  tube  A 
is  too  large  (Figure  6),  the  air  mass  flow  in  the  outer  ring  II  /98 

channel  would  be  greater  than  the  one  obtained  from  the  usual  area 
distribution  law  F,/F,  —  «  •  <  i)  •  The  areas  of  the  holes 

on  the  inner  and  outer  walls  of  the  first  two  experimental  units  on 
ring  I  followed  the  —  .<«.»/ -4k..  law,  due  to  the  fact  that  the 

annular  flame  tube  could  not  maintain  the  balance  automatically  as  the 
pressure  in  the  channel  II  of  the  ring  tube.  Therefore,  three  mutual¬ 
ly  Isolated  channels  were  formed  which  caused  the  actual  external  wall 
inlet  gas  flow  to  be  slightly  larger.  The  symmetry  about  the  center 
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of  the  flame  could  not  be  assured.  Finally,  the  temperature  at  the 
root  of  the  blade  became  too  high  which  is  a  disadvantage  (Figure  7). 

Annulus  I  used  a  reduced  gas  inlet  area  on  the  external  wall 
holes  to  adjust  -^Kw/^Km  to  0.97  in  order  to  obtain  a  more  reasonable 
temperature  distribution  along  the  blade  height. 

In  the  adjustment  of  dilution  holes  it  was  found  that:  when  the 
total  bole  opening  area  was  unchanged,  the  use  of  a  smaller  number 
of  larger  boles  could  obtain  a  better  penetration  depth  than  using  a 
larger  number  of  small  boles.  The  homogeneity  of  the  temperature  field 
was  improved.  Table  1  listed  this  result. 

3-  Total  pressure  loss 

In  the  adjustment  of  annulus  I,  in  order  to  minimize  the  pressure 
loss,  the  rotating  connection  section  of  the  tip  was  changed  from 
Figure  8(a)  to  (b)  to  Improve  the  streamline.  Test  stand  results 

indicated  that:  after  modification,  R  was  raised  by  70  kg  and  more. 

z 

According  to  the  inlet  parameters  of  the  combustion  chamber  of 

type  A  engine,  t^/t^,  structural  parameters  of  annulus  I  and  hole 
d  ^  « 

opening  area,  the  magnitude  of  can  be  estimated  by  using  the  fol¬ 
lowing  equation  and  the  curves  in  [2]: 

-  1  -  K/{2ll  +  CK  -  <  0.3)  (1 ) 

where  'P’/'U"  ,  and  the  flow  resistance  coefficient 

-  </»?  -  pr)/(p,»'i/2)  - 

/  \  * 

Equation  (1)  gave  the  calculated  using  the  maximum  cross-section 
(Acr)  as  the  reference  cross-section. 

« 

Table  2  lists  the  comparison  of  the  estimated  values  of  for 
annulus  I  and  the  R  values  recorded  on  the  test  stand.  Table  3  shows 
the  value  of  P2»  and  measured  with  annulus  I.  The  numbers  4 
and  5  represent  the  measured  results  obtained  under  different  exper¬ 
imental  units  and  different  hole  opening  conditions  as  compared  to 
the  first  three.  Prom  the  estimated  and  measured  values,  we  first 
found  that: 


Table  2.  Comparison  of  the  estimated  value  of  o-i  for  the  annulus 
I  flame  tube  to  the  performance  of  the  engine  stand  test. 


.1 

1 

^  l  ♦ 

If . 

1 

S 

ft 

or 

• 

j  *1' 

0.^21 

0  921  1 

'  4 

^5 

1  .H 

i«  ‘ 

1 

1 

6^5 

tii 

6  4 

V* 

M.IhX 

0.265 

0.255 

1 

0.255 

1 

1 

22. U 

2}.l 

2(6 

«.3 

2».3 

t 

1 — parameter;  2 — experimental  unit;  3 — annulus  1-2;  4 — annulus  1-3; 
5 — annulus  1-4;  6 — annulus  1-3;  7 — annulus  1-4;  8 — remarks;  9 — 
estimated  value;  10 — experimental  value 


Table  3-  Measured  p^,  p^  and  of  the  annulus  I  flame  tube 


1^  tt 

t: 

Z 

Z 

Z 

.  1 

jr 

i 

2 

5 

4  j 

5 

pf 

6.52 

7.00 

7.0g 

7.64 

7.56 

p> 

6.00 

7.30 

7.4g 

7.04 

6.92 

a? 

0.9« 

0.946 

0.957 

0.922 

0.915 

1 — condition;  2 — number;  3 — parameter 


X  .  * 

1)  When  decreases,  is  reduced  and  R  is  also  reduced 

When  stays  the  same,  also  will  remain  the  same.  Thl 

law  has  a  very  significant  meaning  to  the  test  stand  adjustment  of 
flame  tubes . 

2)  The  measured  values  were  very  close  to  the  estimated  values 
and  the  measured  values  were  about  1!?  higher.  The  approximation  of 
the  curve  and  the  accuracy  in  the  measurements  were  the  causes  creat 
Ing  this  error. 
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« 

In  summary,  a_  was  in  the  range  of  0.925  to  0.93  when  using 
annulus  I  for  the  type  A  engine  and  the  measured  values  of  for 
the  seven  stages  were  in  the  0.905  to  0.91  range. 

Main  combustion  zone  lean  fuel  design. 

799 

In  the  adjustment  of  annulus  I,  the  gas  flow  has  been  contin¬ 
uously  increased  in  the  main  combustion  zone.  When  the  performance 
and  temperature  Indicators  were  all  good,  the  gas  flow  rate  in  the 
main  combustion  zone  was  increased  from  22.1%  to  27.8$  of  the  total 
flow  which  was  close  to  30$.  At  this  time,  a  increased  from  0.93  to 
l.l4  which  was  in  the  usual  lean  fuel  design  parametric  range  (1.1- 
1.2). 


Because  the  a  values  in  the  combustion  zone  for  annuli  1-3^2) 

and  I-4^gj  were  large  (1.14),  therefore,  the  were  lower 

than  those  of  the  seven  stage  by  0.01  or  more,  t  ,  was  also  lower 

°  max 

than  that  of  the  seven  stage  by  40°C  or  more.  Furthermore,  <S2qq 
decreased  to  l4$  which  was  22$  lower  than  the  seven  stage.  This 
indicated  that  the  lean  fuel  main  combustion  zone  had  a  higher  heat 
release  rate  than  the  rich  fuel  main  combustion  zone.  The  radial 
temperature  was  more  homogeneous  and  the  hot  zone  temperature  was 
even  lower. 

5.  The  adjustment  of  the  temperature  field  geometrical  relelvant 
parameter  K^g 

Kgg  is  defined  as: 

( 2 ) 

where  L.,  is  the  dimensionless  length  of  the  flame  tube; 

.■  >1.  is  the  ratio  of  the  total  hole  opening  area  of  the  flame 

tube  and  the  cross-section  of  the  flame  tube.  The  value  of  is 

sg 

the  relevant  parameter  for  obtaining  a  given  fuel  gas  temperature 
distribution  in  front  of  the  turbine.  The  larger  K__  is,  the  higher 

So 

the  possibility  of  obtaining  a  homogeneous  combustion  chamber  outlet 
temperature  will  be.  In  addition,  ^200  lowered.  Table  4  shows 
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Table  4.  Comparison  of  K^g  In  the  adjustment  of  annulus  I 


L'  I’i  j 

1 

If  ^ 

1 

tf  1-4 

G 

7 

0.4'i4 

0.  »7| 

0.46< 

^  j 

28.04  1 

2*i.  U 

11.28 

1 

1 

; 

16.8 

15,4 

1 — experimental  unit;  2 — parameter;  3 — seven  stage;  4 — annulus  1-3; 
5 — annulus  1-3;  6 — annulus  1-4;  7 — annulus  1-4 


Tabid  5-  Comparison  of  cooling  air  flow  rate  of  several  types  of 
engines 


4»l  m 

c,../c 

*0.25 

21.88 

7.75 

ASfsisita 

12.90 

10.06 

6.25 

^  u  flttW 

21.80 

17. 16 

8.90 

sniY 

16.40 

12.94 

10.8 

18. 00 

4«».H7 

10.0 

1— parameters ;  2 — engine;  3 — type  A  engine  equipped  with  1-4; 

4 — type  A  engine  equipped  with  seven  stage;  5 — type  B  engine 

that  when  the  values  of  annuli  1-3  and  1-4  Increased,  62^^  was 
reduced  to  5-5?  and  1.455,  respectively. 

When  A  remains  unchanged,  the  value  of  K  increases  with 

« o 

increasing  which  means  that  the  lengthening  of  the  flame  tube 

can  increase  the  duration  of  stay  of  the  fuel  gas  in  the  flame  tube 

to  more  sufficiently  mix  the  fuel  gas  with  the  diluting  air  so  that 

the  homogeneity  of  the  temperature  field  can  be  improved.  When  L 

.  sg 

is  fixed,  to  increase  K  corresponds  to  the  reduction  of  A  .  In 

— o  ®  S 

fact,  this  is  realized  by  reducing  the  flame  tube  hole  opening  area 

to  Improve  the  gas  flow  pressure  loss  which  consequently  will  lead 

to  the  reduction  of  R.  Such  as  in  the  adjustment  of  annulus  1-4, 

when  Kgg  was  changed  from  30. 80  to  31.28  ,  AR^  dropped  by  30  kg.  In 

the  adjustment  of  annulus  I,  K__  «  30  was  chosen  as  the  adjustment 

sg 

base  value. 
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6.  Plow  film  cooling  problem 


Annulus  I  underwent  several  cooling  gas  flow  rate  adjustments. 
After  successfully  accomplishing  all  the  technological  objectives 
of  the  flame  tube,  its  Cl,/C  was  about  23. 88^  (Table  5).  At  this 
time,  the  flame  tube  wall  temperature  was  uniform.  Temperature 
indicating  point  measurements  indicated  that:  70%  of  the  inner  and 
outer  annular  wall  of  the  flame  tube  was  at  550°C.  The  wall  temper¬ 
ature  of  the  inner  and  outer  tubes  was  at  700-750°C.  The  local  wall 
temperature  in  the  main  combustion  zone  reached  720°C.  Annulus  I 
did  not  show  the  local  high  temperature  zone  (850-900®C)  in  the  down¬ 
stream  downwind  of  the  flame  connection  tube  of  the  seven  stage 
flame  tube . 

IV.  CONCLUSIONS 

After  several  adjustment  and  stability  tests,  I  realized  the 
objectives  of  2%  increase  on  R  and  \.^%  decrease  in  C^  (Table  6). 

The  outlet  temperature  was  stabilized.  It  created  conditions  for 
starting  life  test  and  field  test  flights.  Through  the  experiments, 
we  realize  that: 

1)  the  adjustment  of  flame  tubes  must  control  the  gas  flow  rates. 
First,  the  flow  rate  in  the  main  combustion  zone  must  be  adjusted  to 
Improve  the  combustion  efficiency  and  then  the  dilution  zone  flow 
rate  is  adjusted  to  obtain  a  good  outlet  temperature  field 

2)  the  amount  of  air  combusted  in  the  main  combustion  zone,  the 
dimension  of  the  vortex  in  the  main  combustion  zone,  and  the  extent 
of  fuel-air  mixing  may  vary  significantly.  Annulus  I  used  a  lean 
fuel-air  mixture  main  combustion  zone  whose  flow  rate  reached  27.%%. 
Experimental  results  indicated  that  the  main  combustion  zone  gas  flow 
approached  30!?,  i.e.,  a  *■  1.2  is  appropriate. 

3)  the  use  of  larger  holes  in  the  adjustment  of  the  dilution 
zone  performance  was  not  only  applicable  to  ring  tubes  but  also  to 
annuli. 
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Table  6.  Comparison  of  the  average  performance  of  pilot  produc-  /lOl 
tion  annulus  I  and  that  of  the  seven  stage  (same  engine) 
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1 — condition;  2 — parameter;  3 — type  of  flame  tube;  ^ — annulus  I 
(average  of  five);  5 — seven  stage  (annual  average);  6 — increase 
in  design  objective;  7 — the  "fuel  consumption  rate"  with  a  "-" 
sign  represents  that  it  consumes  less  fuel  than  the  allowable  value. 
The  larger  the  absolute  value,  the  smaller  the  fuel  consumption 
rate  is;  8 — 200  point  temperature  field 


The  dilution  zone  air  flow  should  be  30$. 

H)  the  selection  of  a  proper  Ksg  is  a  more  practical  method.  In 

the  adjustment  of  aflame  tube  temperature,  K  should  be  raised  as  much 

sg 

as  possible  to  obtain  a  more  uniform  temperature  field. 


V.  EXISTING  PROBLEMS 

The  rigidity  of  the  external  tube  of  annulus  I  was  weak.  This 
was  determined  by  the  large  diameter  thin  wall  annular  shaped  struc¬ 
tural  characteristics  which  may  affect  the  useful  lifetime.  In  the 
future,  further  steps  will  be  taken  to  Improve  its-  cooling  and  to 
raise  its  rigidity.  The  hot  zone  in  the  temperature  field  carries 
out  certain  movements.  This  is  one  of  the  special  points  which  dis¬ 
tinguishes  the  annular  shape  from  the  ring  tube.  The  radial  movement 
of  the  hot  zone  of  the  tjj  temperature  field  for  annulus  I  was  statis¬ 
tically  determined  to  be  about  15-24®. 

In  the  developmental  stage,  we  have  received  guidance  and  assist¬ 
ance  from  some  comrades  at  Beijing  Aeronautical  Institute.  Those 
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who  also  participated  in  this  work  are:  Su  Xiu-shan,  Zhu  Jin-can, 
Guo  Yong-quan,  etc.  Comrade  Zhang  Qin  has  carefully  corrected  and 
plotted  the  figures  in  this  paper.  We  wish  to  express  our  sincere 
gratitude  here. 

TABLE  OF  SYMBOLS 

A  area,  cross-sectional  area  of  the  flame  tube,  m^ 

Cp^  fuel  consumption  per  unit  propulsion,  kg/kg  propulsion  hour 
E  fixed  condition  (same  as  a  subscript) 

G  combustion  chamber  inlet  flow  rate  Kg/sec 
H  cavity  height  of  the  flame  tube,  mm 

L  length  of  the  flame  tube,  mm 

n  number  of  turns,  proportionality  index  rotation/mm 

R  propulsion,  kg 

t  temperature,  °C 

Z  maximum  condition  (same  as  a  subscript) 

5  Inhomogeneity  of  the  temperature  field  ( 
a  gas  coefficient 

n  combustion  efficiency,  % 

a  pressure  restoration  coefficient,  % 


Subscripts : 

bm 

wall  surface 

Pj 

average 

cK 

reference  cross-section 

sg 

relevant 

ht 

flame  tube 

tb 

tip  of  the  flame  tube 

K 

constant 

W 

outer  ring  of  the  flame  tube 

Kt 

cylindrical  section  of 
the  flame  tube 

yg 

flame  tube 

h^L 

yj 

blade  tip 

Le 

cooling 

ya 

right  engine 

n 

inner  ring  of  the  flame 
tube 

Zu 

left  engine 

3 

combustion  chamber  outlet 

- 

depth  along  the  blade  height 

4 

outlet  of  the  turbine 

radial  average 

200 

200  point  temperature  field  • 

viscosity  coefficient 
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